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An Ultra-High Speed Image Dissecting Camera 
for Photographing Strong Shock Waves 


Kenneth B. Earnshaw and Charles M. Benedict 


(May 
An ultra-high speed image dissecting 


7 cm/microsecond. 


Waves are used to demonstrate the 
framing rate (1> 


camera patterned 
Sultanoff has produced motion pictures of shock waves 


various modes of camera 
10°/second) and the easy adaptation from streak to framing photography 


3, 1962) 


after an original design by 
having velocities greater than 


The camera is readily adaptable to either streak or framing photography. 
A discussion of the camera system design indicates the modifications of the 
which allow for greater flexibility of camera use. 


original design 
Several photographs of high speed shock 
operation. The very high 


make this camera a particularly useful tool for research in the field of very fast shock waves. 


1. Introduction 


The design of a 108 frame/sec image dissecting 
camera by Morton Sultanoff {1},! published in 1950, 
provided a means for observing extremely fast events. 
Although the Sultanoff design is invariably men- 
tioned in comprehensive surveys of high-speed photog- 
raphy [|2, 3, 4], very little information is available 
from the literature concerning the use of this type of 
camera in high-speed shock wave research. It is the 
purpose of this paper to report the results obtained 
from a Sultanoff type camera system developed in 
this laboratory. The discussion will include camera 
and film viewer design, and the techniques used in 
obtaining photographs. 


2. Camera Design 


Since the Sultanoff camera obtains very high fram- 
ing rates by using the principles of image dissec ‘tion, 
perhaps a brief discussion of image dissection is in 
order. 

Image dissection photography, as applied in the 
Sultanoff camera, consists of optically dissecting or 
sampling individual pictures at regular intervals so 
that each picture is composed of a large number of 
lines uniformly spaced over the picture area (as in a 
television picture). If picture sampling line widths 
are narrow compared to the distance between lines, 
many pictures may be recorded on one photographic 
plate by registering each picture grid pattern at a 
slightly different space to give an interlacing of all 
the pictures over a common area. Individual pic- 
tures may be seen by viewing the composite photo- 
graph through a grid composed of transparent lines 
ruled on an optically opaque plate, as in figure 1. 
The grid line widths and line spacings-must of course 
be identical to the corresponding lines and spaces of 
the dissected photographs. By registering the grid 


at different positions across a composite photograph, 
each picture can be seen in turn. 


' Figures in brackets indicate the literature references on page 310. 











As an example, if each picture is dissected into 
lines 13 » wide, with 390 » between lines, as many as 
31 pictures can be displayed on a common area. 
Figure 2 is a composite picture of 10 photographs dis- 
sected as outlined above. The 10 photographs are 
displayed separately in figure 3. Each photograph 
is obtained by indexing the grid, shown in figure 1, 
at the proper position on the composite. 

This example shows only that many individual 
pictures may be interwoven on a single photographic 
plate and that the individual photographs can be 
taken from the composite photograph. In applying 
the principles of image dissection to high-speed fram- 
ing photography, some space resolution is deliberately 
sacrificed in order to achieve time resolution between 
the photographs which are interwoven on a single 
photographic plate. In practice, the necessary time 
resolution is achieved by using a rotating mirror to 
sweep a dissected, time-varying image across a photo- 
graphic plate to obtain the interlacing of many 
images as a function of time. 

Since the basic design principles of the grid framing 
camera are covered thoroughly elsewhere [1, 2], only 
the general design features and the improvements 
incorporated in the camera construction are dis- 
cussed here. Rather complete details of the com- 
mercial components used are presented in appendix 
I. Suffice it to say that the image is dissected by a 
grid, as figure 1, placed at the focal plane of an 
objective lens and that the dissected image is re- 
imaged on film after reflection from a rotating mirror. 
The rotating mirror provides a time dependent 
indexing of successive pictures on the film. 

A pic ‘torial view of the optical system for the ultra- 
high speed camera is shown in figure 4. Although the 
system is patterned after the original Sultanoff de- 
sign, the camera uses a Sinar commercial view camera 
for holding and focusing the first lens. This permits 
rapid focusing and lens changing as well as serving 
essentially as an optical bench on which mirrors, 
filters, shutters, etc., can be mounted easily if de- 
sired. A trap door above and immediately behind 
the multislit focal plane shutter (fig. 4) provides for 
convenient focusing. The door, when opened, auto- 


297 





FIGuRE 1. 


Contact print of an image dissecting grid. 


Slit width is 13 w and the distance between slit centers is 403 ». 


Figure 2. Enlarged print of a composite 


matically swings down a focusing mirror at a 45°‘ 
angle to the focal plane, so that the image may be 
viewed from above. In addition, the camera incor- 
porates a removable grid holder so that grids can be 
interchanged in a matter of a few seconds. The grid 
holder provides for a Fresnel focal plane lens to be 


negative containing 10 dissected pictures. 


positioned just ahead of the grid so that light from 
lens 1 will be redirected to lens 2 (visible in fig. 4). 
The readily changeable grid holder permits rapid 


changeover to streak photography or to other grids 


which are designed for different shutter speeds. 
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FIGURE 3. 


Ten individual photographs taken from the composite in figure 2 by indexing a grid at 10 positions on the composite 


negative. 


ULTRA HIGH SPEED FRAMING CAMERA 








Figure 4. A pictorial view 


The mirror surface is 2 in. 2 in. and is rotated 
about an axis lying on the mirror surface. An 
electric hand grind, ‘ motor rotates the mirror at a 
maximum speed of 600 rps. A photograph of the 
camera is shown in figure 5. 


3. Synchronization 


Since this type of camera is not continuous writing, 
it is necessary to have a precise method of syachro- 
nizing the event to the framing interval of the 


CAMERA BACK 


SYNCHRONIZER —~ 


MULTI-SLIT FOCAL PLANE SHUTTER 


of the ultra-high speed camera, 


camera (the time at which the image is reflected onto 
the photographic plate). Syne -hronization is accom- 
plished by reflec i: a collimated beam of light from 
the rotating mirror (as shown in fig. 4). When the 
mirror is in proper position, the light beam enters 
the photomultiplier. The photomultiplier signal 
triggers a pulse generator which provides a 50-v 
pulse for triggering the shock generating apparatus. 
The synchronizing circuit is wired so that pulses 
cannot leave the camera until a button switch on 
the camera is pushed. Synchronizing in this fashion 
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FIGURE 5. A side view of the camera with the arid holder partiall ; removed 


PHOTOMULTIPLIER-SYNCHRONIZER CIRCUIT 


APICITANCE IN yf 
RESISTANCE IN OHMS 


Si CHARGE SWITCH 
S2 LIGHT SWITCH 
S3 FIRE SWITCH 
PULSE MONITOR OUTPUT 
J2 TRIGGER OUTPUT 


( i MIRROR 
\ 


\ COLLIMATOR 


SEC PULSE 
be 


22v 














FIGt RE 6. Schematic diagram of the photomultiplie rs ynchronize Ps 


takes advantage of an effective optical leverage so | This charges a group of condensors in parallel and 
that synchronization can be quite precise at any | then discharges them in series across the photo- 
mirror rotation speed. Figure 6 is a schematic | multiplier elements. One switching operatioa will 
diagram of the synchronizer circuit. Note that a | activate the photomultiplier for 8 min. A photo- 
simple battery-condenser power supply is used. | graph of the synchronizer assembly (removed from 
The photomultiplier is turned on by switching S-1 | the front of the camera) is shown in figure 7. 


from position 2 to position 1 and back to position 2. 
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Figure 7. Photomultiplier synchronizer (removed from front 


of camera) 


4. Film Viewer 


The film viewer for unscrambling the composite 
photographs produced by the camera is shown in 
figure 8. Correct alinement of the photograph on 


the grid is accomplished by the alinement adjustment 


which rotates the disk shaped film holder. The 
disk is anchored in the top carriage which is moved 
over the base plated by a cam. The top carriage 
is spring loaded against the cam and is held to the 


FILM HOLDER- 


ALIGNMED 


base plate with dovetail slides. The grid is held 
in a countersunk hole in the base plate. Indexing 
of the photograph across the grid is accomplished by 
turning the indexing adjustment. The worm gear 
drive (located under the cam) and the cam taper 
are such that one revolution of the knob produces 
a carriage movement of 0.001 in. Carriage position 
is easily read on the dial indicator. Photographs 
are normally viewed by placing the viewer over a 
light source. However, the entire film viewer is 
mounted on a board which can be placed in an 
enlarger for copy work. 

For moving picture photography, the film viewer 
is motor driven as shown in figure 9. The variable 
transformer motor control and the three speed pulley 
drive provide for a range of viewing speeds of about 
60:1. This range has been satisfactory for all 
photographs viewed so far. Moving picture records 
are obtained by positioning a movie camera above 
the viewer with a light source placed under the 
viewer. Figure 10 shows a typical arrangement for 
moving picture recording. 


5. Shock Wave Photography 


In any photographic observation of shockTwave 
phenomenon, it is useful first to obtain streak photo- 
graphs of the phenomena of interest. From the 
streak photographs, measurements of luminosity 
times permit one to calculate the optimum mirror 
speed for framing photography. 

Streak photographs are obtained by replacing the 
grid with an opaque plate having a single trans- 
parent slit. For the work presented here, a rather 
coarse grid was used having all but the center slit 
masked off. Good alinement was easily accomplished 
by imaging the shock tube on the center line of the 


DOVETAIL SLIDES 


/ 
INDEXING AOJUSTMENT—~ 


FIGURE 8 


Film viewer for observing high speed photographs. 
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FIGURE 9 


erid before the grid was masked off. Figure 11 
shows a series of streak photographs of shock waves 
taken at various The horizontal lines 
were made by manually alining the mirror so that a 
picture of the full grid was put on the film. This 
provides a convenient vertical time scale on the 
streak photographs which are calibrated from an 
oscilloscope record of synchronizer pulses generated 
by the camera. A horizontal distance scale is made 
by placing narrow strips of tape on the shock tube 
at suitable intervals. For convenience, a4in. x 5in. 
Polaroid film holder is used for streak photography. 
Before taking framing photographs, a calibrated 
streak photograph is used to estimate the average 
time (f,yg) for the luminous slug to pass a point on 
the shock tube, and the total time (t,) for the lumi- 
nous slug to clear the region in which framing photo- 
graph investigations are to be made. These 
estimates are diagramed in figure 12. The maximum 
rate is determined by the grid dimensions and by the 
total luminous time (t;). With a grid having a slit 
spacing )), the maximum permissible mirror rotation 
frequency without double exposure is given by: 


pressures. 


D 


4nrt 

Optical arm of camera (inches). 

Total luminous time (seconds). 

Mirror rotation frequency (cycles per second). 


Grid slit spacing (inches). 
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of luminosity which are clearly 


Film viewer with motor drive. 


Figure 13 shows a sequence of single frame photo- 
graphs without double exposure. Although no 
double exposures occur on this photograph, resolu- 
tion of the luminous front of the shock wave is rather 
poor because of the relatively low framing rate used. 
Still, the luminosity form and_ velocity 
observed quite well. 

In general it is desirable to obtain double expo- 
sures so that velocity measurements and changes 
in luminosity patterns can be observed on one photo- 
graph. In addition, by allowing double exposures 
higher framing rates can be achieved. However, 
the maximum camera speed is limited by the time 
(tave) for the luminous slug to pass a point on the 
shock tube in the region of interest, and by the grid 
spacing D. Maximum camera speed is given by: 


can be 


D 
Anrtave 
Mirror rotation speed (cycles per second), 
Grid slit spacing (inches). 
Optical arm of camera (inches). 


Maximum time (seconds) for 
slug to pass any point 
interest. 


luminous 
in region of 


Figures 14, 15, and 16 show sequences of mutiple 


exposure photographs. Note the complex patterns 
visible with the 
better resolution of detail apparent in these photo- 

















Fiagure 10. Arrangement for making high speed motion pictures. 
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Po= SOM Hg 
TIME SCALE: 0.836 »SEC/DIV 


Po=!OO0n Hg 
TIME SCALE: 0.836 “SEC/DIV 


TIME SCALE: 0840 p2SEC/DIV 


Po>= 300 Hg 
TIME SCALE: 0.840 ,uSEC/DIV 


Po = 500 Hg 
TIME SCALE: 0844 puSEC./DIV. 


Po=!000p Hg 
TIME SCALE: 0848 pSEC/DIV 


DISTANCE IN CENTIMETERS 


FIGURE 11. Streat photographs of shockwaves in electrically driven helium shocktube (taken at various ambient pressures). 


In each case energy input is 2000 joules. The shock tube is 4.15 em i.d. X 122 em. long 
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DISTANCE IN CENTIMETERS 


FIGURE 12. 


graphs. The luminous vertical lines observed in 
figures 15 and 16 are caused by the shock wave 
colliding with 3 mm o.d. glass tubes placed across the 
inside of the shocktube and spaced at 10-cm inter- 
vals. To better interpret these collisions, a mirror is 
used to obtain two simultaneous views of the shock- 
tube at right angles to the axis of the tube. Figure 
16 shows the typical overlapping of images (on the 
left) caused by the luminous slug being too long to 
be displayed separately at the camera speed 
employed. 

Figures 13 through 15 show only 6 out of a possible 
31 sequences available. Moving pictures may be 
obtained by recording all of the possible photographs 
by using the motor driven film viewer with a 16 mm 
moving picture camera. A strip of moving picture 
film taken by the arrangement in figure 9 is shown 
in figure 17. 


6. Film 


In general, the best photographs of highly lumi- 
nous shock waves have been obtained on Eastman 
Kodak Tri-X 4 in. x 5 in. glass plates. For events 
of lower luminosity, good photographs have been 
obtained on Kodak Spectroscopic 1—D 4 in. x 5 in. 
plates. Better results could no doubt be obtained 
with plates checked for flatness; however, no flat 
plates were immediately available for use in our 
work. Polaroid 4 in. x 5 in. packets (types 52 and 


653930—62 


| speed shock waves. 


Estimates of the duration of luminosity from a streak photograph. 


57) used in a Polaroid 4 in. x 5 in. film holder have 
been used for all streak photography. The Polaroid 
film is also useful in checking camera synchroniza- 
tion and for rapidly estimating the correct exposure 
for framing photography. 


7. Summary and Conclusions 


The Sultanoff type of high-speed image dissecting 
camera has proved to be a very useful tool in ob- 
taining streak and frame photographs of very high- 


Motion pictures may be ob- 
tained of events too fast to record by any other type 
of motion picture camera. In addition, if some 
sacrifice in picture detail can be tolerated, the sim- 
plicity in design of this camera allows for a much less 
expensive fabrication than is possible in other types 
of high-speed cameras. By using a variety of grids, 
streak and framing photogr: aphs | may be taken of a 
wide variety of se f-luminous events, and throughout 
a broad range of framing rates. 

Some of the disadvantages of the camera should 
be mentioned. The event to be photographed must 
be synchronized to the camera, since it is not con- 
tinuous writing. In addition, the achievable camera 
speed is limited by the time required for luminosity 
to decay appreciably at any given point of interest. 
For example, events which are strongly luminous for 
long periods of time will place serious limits on the 
framing rates that are possible. The higher framing 
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20 99 60 65 10 1$ 
DISTANCE IN CENTIMETERS 


FIGURE 13. 1 sequence of six single exposure photographs of a shock wave in helium at 175 »w Hg, using the same shocktube as 
in figure 11. 


The effective framing rate is 5.57 X 10¢/second. Energy input is 2000 joule 
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Figure 14 1 sequence of six multiple exposure photographs of a shock wave in helium at 150 p Hg, using the same shocktube as 


in figure 11. 
The effective framing rate is 48X10%/second. Shock waves are numbered in time sequence. A complex interaction of the shock wave with the shocktube walls may 
be observed in these photographs 
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FicureE 15. 1 sequence of six multiple exposure photographs of a shock wave in helium at 50 w Hg. 
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FiGuRE 16. A sequence of six multiple erposure photographs of a shock wave in helium at 50 p Hg, using the same shocktube as 
in figure 15 


The effective framing rate is 28.2 * 10*/second 





equence of 135 pictures taken from a 16 mm moving picture film of the event displayed in figure 13. 


rates in this case can be achieved only by further 
sacrifices in resolution by using a grid having a wide 
spacing between slits or by using a smaller enlarge- 
ment (or greater diminution) of the event being 
photographed. In other words, for slowly changing 
luminosity patterns, a greater sacrifice in resolution 
is necessary to obtain the higher framing rates. 

Although these disadvantages do exist for certain 
self luminous events, the Sultanoff camera is par- 
ticularly well adapted to motion picture studies of 
very high velocity shock waves. 


The authors express their thanks to C. M. Shep- 
herd for his help in designing the viewer and to W. R. 
Saxton, R. D. Dickason, and T. L. Theotokatos for 
their help in obtaining the photographs presented in 
this report. 
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9. Appendix I. Components Used in the 
Ultra-High Speed Camera 


A numbered layout of the high-speed camera is 
given in figure 18. The numbered parts are: 

1. The basic camera frame, the rotating mirror, 
and the grid holders 1(a) were manufactured by the 
Buck Instrument Company, Boulder, Colo. The 
image dissecting grids 1(b) were manufactured by the 
Bausch & Lomb Optical Company. The 8 in. focal 
length 4 in. x 5 in. Fresnel lens 1(¢) may be ob- 
tained as a standard item at almost any photographic 
supply house. 

2. The mirror drive motor is a Precise Super 40 
high-speed grinder motor. 
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\2(c)  12(a) 


FIGURE 18. 


110V, 60 ~ Variae motor speed control. 








A layout of high speed camera parts, 


10. Schneider-Kreuznach Symmar, F/5.6 5.6 in. 


4 in. x 5 in. Graphic photographic plate holder. | focal length. 


Polaroid 4 in. x 5 in. flat film holder. 
}. Tripod support; the main camera is supported 


on a movable stand and the tripod is used only for 


supporting the front of the camera at the tripod 
mount (12b). 

7. Linhof tecknika 4 in. x 5 in. adjustable and 
revolving back. This camera back is 
furnished with a bellows which has been cut off short 
and attached to the camera frame. The camera 
back is mounted on a frame which is held to the main 
camera frame by three adjustable screws. Correct 
alinement of the camera back is accomplished by 
adjusting the three screws (which are locked in 
position after alinement is correct). 

8. Kodak Aero-Ektar lens, F/2.5, 7 in. focal length. 
9. Ilex Paragon anastigmatic, F/4.5, 12 in. focal 
length. 


camera 


11. Re-imaging lens (lens 2 in fig. 4): Georz Artar, 
F/9.5, 10%4 in. focal length. 
12. Sinar commercial view camera, including: 


(a) and (b) Two bench holders with sleeves. 
Bench holder (b) is used with the tripod. 
(c) 6 in. extension bench unit. 
(d) 18 in. extension bench unit. 
(e) 12 in. extension bench unit. 
(f) Joint block. 
(¢) 4 in. x 5 in. rear carrier. 
(h) Front principal frame, with 3 lens boards 
(i) Auxiliary frame. 
(j) 4 in. x 5 in. tapered bellows. 
(k) Square bellows extension. 
(1) 4 in. x 5 in. wide angle bellows. 
(m) Front surface mirror. 
( Paper 66C4-105) 
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Biprism Method of Determining the Equivalent Focal 
Length of Flat Field Lenses 
Walter R. Darling 


(June 25, 


1962) 


A device is described that permits the rapid determination of the equivalent focal length 


of a lens. 


A transmitting biprism mounted between a collimated light source and a lens, 


divides the light incident upon the front of the lens into two parallel beams making a fixed 


angle with one another. 
plane. 


On passing through the lens, two images are formed in the focal 
The magnitude of the lateral separation of the images is determined by the angular 
separation of the two incident beams and the focal length of the lens. 


The focal length of 


the imaging lens may be determined from the measured separation of the images at the focal 
plane of the lens and the known angle of deviation of the two incident beams produced by 


the biprism. 


1. Introduction 


The National Bureau of Standards is frequently 
valled upon to determine the equivalent focal length 
and related constants of lenses with a high degree 
of accuracy. For many types of lenses these meas- 
urements are customarily performed visually on the 
precision optical bench to within +0.10 mm on 
lenses whose equivalent focal lengths range up to 
200 mm, +0.15 mm for lenses whose equivalent 
focal lengths range between 200 and 800 mm. +0.25 
mm for lenses whose equivalent focal lengths range 
from 800 to 1,200 mm, and +0.50 mm for lenses 
whose equivalent focal lengths range from 1,200 to 
1,800 mm.' 

Frequently it is desirable to determine only the 
equivalent focal length of lenses for use in evalua- 
tion of certain metrical qualities, such as distortion 
and resolution at finite and infinite distances. For 
these determinations the values of the focal lengths 
need not be as accurately known as those obtained 
by the precision optical bench method, yet it is 
desirable to know these values to a greater degree 
of accuracy than the value marked on the front of 
the lens, which is referred to as the nominal focal 
length. The tolerance for the nominal focal length 
is given by the ASA Standards? as +4 percent. 

There is small likelihood of gross error in optical 
bench methods, but to preclude the possibility of 
such error it is desirable to have a method that per- 
mits a quick check of the value of the equivalent 
focal length of a given lens obtained with the optical 
bench. It is moreover desirable that the check 
measurement be based upon a different principle of 
measurement and employing different data to mini- 
mize the danger of systematic error. In this paper 
a simple check method for measuring the equivalent 
focal length of a lens is presented. A biprism placed 
in front of the lens under test splits the incoming 
collimated beam into two parallel beams with a 


—$ ee 

1 The differences in accuracy depend more upon the relative apertures than 
upon focal lengths, as a consequence of the differences in depth of focus. 

2 ASA Standards—Focal Length Marking of Lenses. PH.3.13-1958. 











fixed angular separation. With this prism, values 
of focal length accurate to within +1 percent may 
be obtained. While the range of uncertainty is some- 


what higher than that attained by the optical bench 
method it nonetheless permits one to evaluate the 
equivalent focal length of a given lens with a degree 
of precision well within the tolerances set forth in 
the ASA Standards for nominal focal length. 


2. Description of the Instrument 


An optical wedge of borosilicate optical glass 
approximately 2% in. in diameter was ground and 
polished with an angle of approximately 2 degrees 
between the faces. The finished wedge was then 
cut on a diamond saw in the principal section across 
its diameter. One half of the wedge was inverted 
with respect to the remaining half so that the devia- 
tions were in opposite directions. These two halves 
were then cemented to a flat base of ground and 
polished optical glass of approximately 3% in. in 
diameter and % in. in thickness. For small wedge 
angles a, the deviation angle 6 for normal incidence is: 


6=(n—1)a, (1) 


n being the index of refraction of the material used 
for the biprism. For this borosilicate glass (n—1) 

0.519. For a collimated beam of light incident 
normally upon one side of the biprism, the angle 
separating the two emergent beams is 26. The 
magnitude of the angular deviation, 26, was measured 
in the refractometry laboratory with white light 
and was found to be 2.059 degrees. The probable 
error of this measurement is approximately 6 sec. 
The rest of the instrument consists of a cell to hold 
the biprism and a rod attached to the cell for adjust- 
ment of the entire unit. The component parts of 
the instrument are illustrated in figure 1. 

The wedge is not achromatic but the error re- 
sulting from the lack of achromatism is ‘negligible. 
When white light passes through a prism (fig. 2), the 
components of different wavelengths are deviated by 
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FIGURE 1. The component parts of the instrument are 


lustrated. 





»= 64.5 


4 
64.5 


FIGURE 2 When a white 
the difference 
heam 


ight beam passes through the wedge, 


in the deviation of the red beam and the blue 


forms an angle é. 


different amounts. The magnitude of the disper- 
sion ?- customarily measured by the 
between the deviations for red light, hydrogen 656, 
and blue light, hydrogen 486. Let this difference in 
deviation be called 6. It may easily be deduced 
from eq (1) that 


the Abbe 


made of 


number of 
borosilicate 


where y is 
biprism is 
Consequently 


the glass. Our 
glass, v—64.5. 


6 0.0158. 


The wavelength of maximum visibility, for which 
visual settings are made when white light is used, falls 
almost midway between the wavelengths F and C 
Since light at F and at C is definitely colored and 
also considerably less bright than that of maximum 
visibility, it is very unlikely that a setting would be 
off as far as 4% 6. Further, the deviation of the bi- 
of deviation to the left of 6 plus 


prism consists 
right of 6. Consequently the dis- 


deviation to the 
persion is also half in one image and half in the other. 


The error in setting, on one of the images, resulting | 


from chromatic aberration should not exceed 


difference | 


» where 26 is the total deviation of both prisms, 


Vv 
and the error in the measured separation of the two 
images from this then should not exceed 

1.47 
(4) (64.5) 
approximately one quarter of 1 percent of the meas- 
ured interval. This represents an upper limit of 
error that should never be exceeded, not the probable 
error, which would be one fourth as large. 


source 


» where x is the separation of the images, or 


3. Method of Measurement 


The lens under test is placed in the chuck of the 
nodal slide assembly of the optical bench. The 
nodal slide assembly is moved along the ways of 
the optical bench until the image of a collimated 
source formed by the lens under test is viewed in 
the microscope of the optical bench. The reticle 
used in the collimator consists of a vertical and 
horizontal line at right angles to each other forming 
a cross. The nodal slide asse mbly is now adjusted 
by the lead screw of the optical bench until the 
image formed by the lens of the collimated beam 
proceeding from the illuminated reticle is in sharp 
focus. The biprism is now inserted close to the 
face of the lens under test between the lens and 
the collimated beam from the illuminated reticle. 
This is illustrated in figure 3. The image of the 
reticle formed by the lens under test is split by the 
biprism into two images which are viewed in the 
ocular of the viewing microscope. Figure 4 is a 
photograph of the two split images as seen through 
the microscope. Care must be exercised in the 
orientation of the biprism to achieve maximum dis- 
placement of the split images. This orientation is 
produced by rotating the biprism within its cell 
about its optical axis until the two horizontal lines 
of the viewed images, if continued, would appear as 
one continuous line. The horizontal crosshair in 
the viewing microscope is used as a guide in making 
this adjustment. The biprism has no spherical 
refracting power so the position of best focus is not 
changed by the introduction of the biprism. The 
viewing microscope is equipped with a lateral ad- 
justment and scale that allows measurements to be 
made to within +1 uw. The vertical arm of the 
crosshair in the microscope is brought into coinci- 
dence with the vertical line in the image of one of 
the two images formed by the biprism and lens by 
the lateral adjustment of the viewing microscope. 
A reading Ff, is taken of this position on the lateral 
scale of the microscope. The microscope is now 
traversed laterally until the vertical line in the cross- 
hair of the microscope is in coincidence with the 
vertical line in the second image of the cross and a 
reading /?, is taken on the lateral scale of the viewing 
microscope. Readings Ff, and FR, are successively 
taken until an average of at least five readings of 
each position has been obtained. The focal length 
can then be determined from the following relation: 

D= Rya— Rag 


2f tan @ 
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4. Precision Optical Bench 


The operation of the precision optic al bench has 
been discussed in a previous paper.’ It may be 
mentioned here that the lens under test by the opti- 
cal bench method of measurement is subjected to a 
number of critical time-consuming adjustments be- 
fore any measurements are made. Depending upon 
the physical size of the lens and the focal length as 
much as two hours time is needed for adjustment of 
the lens on the precision optical bench. 


5. Results of Measurement 


The following table contains the results of measure- 
ments made on a number of lenses of varying focal 
lengths. 


TABLE | 


Figure 3. The briprism is inserted close to the face of the lens Webues stamnmiaioad ii 
under test, between the lens and the collimated beam from the focal length difference 
illuminated reticle. 


Lens Nominal Measured 


by optical by using 
bench biprism 
method 
f lo 


mm mm mm 
5 5&8. 3f 58. 38 
58 5S. 36 58. 37 
5 58. 5 58. 79 
58. 22 

58. 09 

58. 42 

58. 22 

149. 94 150, 29 

180. 50 179. 38 

305. 83 306. 52 

486. 58 486. 67 

610. 67 610. 59 

763. 09 763. 44 

1065. 55 1064. 99 


Average deviation 


Improper orientation of the prism may lead to 
error, but even this is unlikely to lead to serious error 
as long as the prism axis does not deviate from 
normality with the ways of the measuring micro- 
scope by an amount exceeding +2.5 degrees. A 
misorientation of this magnitude is not likely to 
occur as the splitting of the horizontal arm of the 
cross would immediately make manifest this type 
of error. 


Fiaure 4. <A photograph of the two split images as seen through 
the micros¢ ope 


6. Limitations in Measurements 


The maximum travel of the viewing microscope 
available for the measurement of ?' and FP? was 50 
mm. With the biprism whose angle of deviation 
being 2.059 degrees, the maximum equivalent focal 
where D=R,,—R,, or the difference of the averages | length that can be determined is 1390 mm. Longer 
of the lateral scale measurements of R, and R). | focal lengths can be determined by decreasing the 
The quantity 0.5 cot @ is a constant of the instru- | angle of deviation in the biprism or by increasing the 
ment. Hence the value of the equivalent focal | lateral travel of the viewing microscope. 
length may be obtained directly from the formula, | ——— 

F , : 3 F, O. Washer, W. R. Darling, Factors affecting the accuracy of distortion 
| en ee on the nodal slide optical bench, J. Opt. Soc. Am. 49, 517 


0.5 D cot 6 


Care must 


7.81 D. (2) 
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be exercised in using a prism of this nature, one must 
use a reasonably monochromatic light source, and 
the initial calibration of the prism must be made 
with near monochromatic light of approximately the 
same wavelength as that which will be used in the 
focal length determination. 

Upon examination of the results of measurements 
obtained in the experimental data the focal length 
determinations by the biprism method were found 
to be biased. The values reported in table 1 using 
the biprism method were adjusted in the amount of 
—0.07 percent. A 0.07 percent change in the focal 
length amounts to approximately 5 sec change in 
deviation of the biprism, this is well within the toler- 
ance of the measurement of the biprism. With the 
deviation from the average not greater than +0.15 
percent one can be confident that a focal length de- 
termination by means of the biprism is very un- 
likely to be error by as much as 1 percent. 

The focal length determinations in table 1 were 


made with a Wratten No. 73 filter with an effective 


wavelength of approximately 575 mu. 
7. Sources of Error 


Since {=0.5D cot 6, error can be produced by 
errors in D and by errors in @. 
without error, from the relation in eq (1) 0.5D cot @ 
is a constant equal to 27.81. Therefore the error in 


Assuming @ to be | 


the equivalent focal length as a result of error in the 
lateral scale readings is given by the expression: 
Af=27.81 AD, eq (2). 

From the above it is clear that a 2u error in D 
leads to an error of +0.055 mm in f which can 
usually be considered negligible. Errors in f arising 


from errors in the measured value of @ can be neg- 
lected as it is unlikely to exceed +0.2 percent for a 
probable error of 6 sec in an angle of 2.059 degrees. 


8. Discussion 


It is evident from the information contained in 
the foregoing sections that the biprism method is a 
satisfactory means for quickly determining the equiv- 
alent focal length of flat field lenses with an error 
that does not exceed +1 percent. It is also pointed 
out that a nodal slide assembly is not a necessity in 
making the focal length determination. Any method 
that allows the lens to be held firmly and in a reason- 
able alinement with the collimated incident beam is 
sufficient. 

—a— 
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The effect of air drag on the motion of a filament struck transversely is treated theo- 


ret ically. 


The air drag is shown to produce curvature in the transverse wave formed by the 


impact, to increase the speed of the transverse-wave front, and to increase the strain in the 


filament. 
meters/second. 


The theory is applied to the case of a nylon yarn impacted transversely at 189 
The calculated radius of curvature of the yarn was 1.1 meter, which agreed 
well with 1.2-meter radius obtained by experiment. 


The calculated effects of air drag on 


the strain and distance traveled by the transverse-wave front 601 microseconds after impact 


were small. 


When a long filament is subjected to a high-veloc- 
ity transverse impact, strain waves and transverse 
waves are generated which travel outwards along the 
filament away from the point of impact. According 
to theory [1—6]' the velocities and strain distributions 
in these waves can be calculated provided that the 
stress-strain curve of the filament under the impact 
conditions is known. Recently [7-9] the inverse 
problem has been considered for a textile yarn, i.e., 
the calculation of the stress-strain curve if the veloc- 
ities or the strain distribution resulting from an 
impact are known. 
incomplete for this practical application because it 
neglects the effect of air drag on the yarn. 
paper the theory of transverse impact is extended to 


include the effect of air drag, and the errors intro- | 


duced by its neglect are estimated for a typical case. 


1. Behavior of the Filament Neglecting Air 
Drag 


Since the effect of air drag on the motion of the 
filament is treated as a perturbation on the motion 
without air drag, the theory neglecting air drag is 
first reviewed. 
strain waves that propagate outward in each direc- 
tion from the point of impact. These strain waves 
produce a strain ¢ in the filament. In the region 
between the strain-wave fronts, material of the fila- 
ment is set into motion longitudinally toward the 
point of impact with a velocity of Wy. This material 
is tuken up by a tent-shaped wave of transverse 
motion, as shown in figure 1. The filament material 
between the projectile and the transverse-wave front 
moves vertically with the velocity V of the projectile, 
and the wave front travels in a longitudinal direction 
with velocity 1%. The section of filament between 
the projectile and transverse-wave front is a straight 
line; the solid and dotted lines in figure 1 show the 
filament at two different times after impact. 


1 Figures in brackets indicate the literature references at the end of this paper. 





Unfortunately, the theory is | 


In this | 


The impact produces longitudinal- 


The longitudinal velocity of the filament material, 
Wo, is given by 


] "ey dT F (1) 
- - ( 
yn J0 \ de ; 


where m is the mass per unit length of the unstrained 
filament, é is the strain, and 7'is the tension required 
to produce a strain e. 

The transverse-wave front velocity is 


Wo- - 


X am ae T > 
WV m(l+e,) 


with respect to the unstrained filament, i.e., with 
respect to a Lagrangian coordinate system. 7) is 
the tension in the filament resulting from the impact. 
These velocities are related to the impact velocity V 
by the equation 


(2) 


V2=2W,(1+ 6) Ao— W?. (3) 


An expression for the transverse-wave front 
velocity (’) in laboratory coordinates will now be 


\ 
fr. 


i 7 
PROJECTILE 
I~ 





Ficgurs 1. Configuration and motion of a filament impacted 
transversely at velocity V, air drag neglected. 
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derived. Suppose ink marks 1 cm apart are placed 
on the unstrained filament before the filament is 
impacted. After impact, the transverse-wave front 
velocity in cm/sec is Ao with respect to the unstrained 
filament, so the transverse-wave front will pass do 
ink marks per second. However, these marks are 
now (1+ é)) em apart since the strain in advance 
of the transverse wave is é. The transverse-wave 
front, therefore, traverses a distance dA (1+e,) 
em on the filament per second. Since the filament 
in advance of the transverse wave has a velocity W5 
opposite to the direction of motion of the transverse- 
wave front, a stationary observer detects a 
verse-wave front velocity of 


trans- 


W. (4) 


The slope, po, of the filament behind the transverse- 
wave front is seen from figure 1 to be given by 


tan y= V/U>. 


Equations (1) to (4) describe the motion of 
impacted filament in the absence of air drag. 
equations, however, must be modified slightly in 
certain For instance, the above theory is 
based on the assumption that the stress-strain curve 
is always concave downwards, whereas many stress- 
strain curves for textile yarns have concave upward 
portions. The modification required in this case 
has been discussed by Smith et al. [10]. At very 
high impact velocities, part of the region of variable 
strain in the strain wave occurs in the wake of the 
transverse-wave front instead of in advance of it as 
assumed in the discussion above. The modifications 
required in this rarely occurring case [6] are not 
considered in this paper. 


2. Effect of Air Drag 


The force on the filament between the projectile 
and transverse-wave front due to air drag may be 
expressed as a component perpendicular and a 
component parallel to the filament. However, the 
force component parallel to the filament is small, so 
only the force component perpendicular to the 
filament will be considered.’ 

This force produces curvature in the section of 
filament between the projectile and transverse-wave 
front, thereby increasing the strain and tension in 
the filament. This increase in strain produces strain 
waves of small amplitude that are propagated past 
the transverse-wave front and along the filament. 
Since the velocities of these strain waves are usually 
much greater than the velocity of the transverse 
wave, the strain, e, in the filament between the 
projectile and transverse wave is nearly uniform at 
any time, but varies with time. Thus the Lagran- 
gian velocity of the transverse-wave front is given 
by 


the 
These 


cases. 


T 
' \ m(1 


e components are 


pared later in the paper 


where 7 is the tension in the filament. Since e and 
T vary with time, \ also varies with time. 

In order to compute the length of the filament 
between the projectile and transverse-wave front, 
consider ink marks that have been placed 1 em apart 
on the filament before impact. The rate at which 
the transverse-wave front these ink marks 
is X\em/s. Therefore, the transverse-wave front will 


bd | 


passes 


have passed Adt marks at a time 


Jo 
Since the varn has strain e at time 
between the marks will be (1+e) em. 
the length of filament between the 
transverse-wave front is 


t after impact. 


the distance 
Therefore, 
projectile and 


dt. (6) 


Note that because of the extra strain caused by air 
drag, this length of filament is no longer a straight line 
but is curved. 
The velocity of 
oratory 


the transverse-wave front in lab- 
coordinates is given by \(1+e)—W and 
time. ‘Fherefore, the longitudinal dis- 
tance traveled by the transverse-wave front is 


varies with 


(7) 


where the velocity of the material of the filament in 


advance of the transverse-wave front is given by 


—{ ar de (8) 
ym. Vd 


Equations (5) to (8) will be used to calculate the 
distances S and _X after e and 7 are given as functions 
of time. 


3. Evaluation of Air Drag 
Let the strain e be expanded to give 
Cotadotat+at?+ .... (9) 


The strain is assumed to vary only a small amount 
from the strain é) that would occur without air drag, 
so a few terms of this series are expected to give a 
good approximation. For strains the 
tension may also be approximated by 


near po, 


T,+ K(e—é) (10) 
where K is the slope of 
the filament at @p. 
filament is 


the tension-strain curve of 


Therefore, the tension in the 


Ka, + Ka,t+ Kayt?-+ (11) 





With this 


approximation for the tension, eq 
becomes 


: ae 


I "0 dT 
d 
ym. ey \ de de 


jo 3a 


W,-4 \K me €o9). 


W e-+- 


ym 
(12) 


Equations (9) and (11) 
egs (5), (6), | 


as series 1n f. 


are now substituted in 


{ 
7 
Thus eq (5) becomes 
T,+Ka,+Ka,t- 
a,t- 
T / 1+kay+ka,t 4 
/ Io t 
V m(1+e5) \ poy at 
l+e, 1+6 
ka, 


1+kay 
ayt 


+-@_9+o-4 


1+ 


Teo Tilo 


where k= K/ 7). 
Dividing and taking the square root of the series 
vields 


(k t ke, i Codotit | 42 
2V1-4 ka, (1+-e9+a,)?/? 


l)yl 


(14) 
where Of? represents terms of the second and higher 
powers of ¢. The following equations are similarly 
derived: 


T kao 


(1+-e,+4a,)k ov 


t-adoy 1+kay 


tT ot yy 1 


ov 1+ €ov 1 T Cot oy | t kao 


: Kk 
Wo Vn ~ 


Cot 2ao)k Kat 
1 +e Ap —-24/- 1 +-0t?. 
T dy v1 T kao ’ \ m 4 + t 
(16) 


In order to solve for do, a), etc., a relationship must 
be found between S and X. This relationship is 
derived by considering the configuration of the 
filament between the projectile and transverse-wave 
front. 

The configuration of the filament exposed to air 
drag is shown in figure 2. Let s be the distance 
measured along the filament between the projectile 
and a segment ds of the filament. 

The air drag on the segment ds is composed of a 
frictional component due to the viscosity of the air 
and a dynamic pressure component due to uneven 


), and (12) to express A, S, X, and W | 











FicgurE 2. Effect of air drag on the shape of a filament im- 


pacted transversely at velocity V. 


it, 
aT (S+ds) 


? 


FicgurRE 3. A segment ds of the transverse wave portion of the 
impacted filament showing tensile forces required to balance 
the air drag. 


air pressure around the diameter of the yarn. The 
frictional component of the air drag is only a few 
percent of the dynamic pressure component [11] and 
so will be neglected. By the “cross-flow principle” 
[11] the dynamic pressure component of the air drag 
is directed perpendicularly to every point of the 
filament. Figure 3 shows the forces on a segment ds 
of the filament. The equilibrium equations for the 
horizontal and vertical force components of the 
segment are 

d(1 0069) sing 0 
ds 
and 

d(T = Y) 4 F cosy 0. 

as 


where F is the force per unit length on the filament 
due to air drag. These conditions simplify to 


dT 0 


and 


(18) 


Equation (17) states that the tension is constant 
along the filament. 

The force due to air drag is proportional to the 
square of the velocity of the segment, to cos? ¢ where 
¢ is the slope angle of the segment, and to the length 
of the segment. Since the case of small air drag is 


3 If the slope of the filament is close to 90°, the pressure compcnent of the drag 
becomes small. 
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being considered, the velocity may be approximated 


by the projectile velocity V and the slope ¢ by the 
slope ¢ that the filament would have if no air drag 
were present. Therefore, the force on the element 
ds is 


Fds (19) 


qV? cos® gods. 
Since F and T are independent of s and therefore 
constant along the yarn, eq (18) shows that the 
filament forms an are of a circle of radius 7/F. 
The angle subtended by the are S at the center of 
this circle is given by the length of the are divided 
by the radius SF/7, where S is the length of the 
filament between the projectile and transverse-wave 
front. The length of the corresponding chord of 
the circle is then (27/F) sin (SF/2T) and 
yX*+°V @ by figure 2. Therefore, 


is also 


(27/F) sin (SF/2T)=yX?4+V7#. 

Since F' is very small compared to 7, the sine may 
be approximated by relatively few 
infinite series. Squaring both sides of the equation 
then gives 


X24 


(20) 


(neglecting terms higher than F?), which is the de- 
sired relationship between Y and S to evaluate the 
coefficients do, a, ete. 

Equations (15), (16), and (3) for S, X, and V are 
substituted in eq (20) to yield the following identity 
in ¢. 


\ K ma, 


& l T Coy ] Tégtaoy ] + ka, W, 


+ 


Es + &% 


3+(1-4 ot 4ay)k 


dov1+kay 


\ 1 Fant +002 | +(Q#?=(0. (21) 


Since this expression is zero for all values of ¢, the 
constant part and the coefficient of ¢ must each be 
zero. Simplification of the constant part gives 


V¥1-+ep+a9 y1+ka a 
4.1 +e, 


Kay 


2m 


— kK, MA 1 T E¢ \ 1 régtaoy i+kas 


WK, m |a 1 Wdov 1 t €oly1 Tt Ey 


V1 +e9+doV1+kao| =0. 


terms of its | 


By inspection a,=0. Therefore, by eq (9), the 
strain just after impact, when ¢ is very small, is 
seen to be é@, which is the strain for no air drag. 
This is reasonable since the length of filament in 
transverse motion just after impact is too short to 
be affected by air drag. 

Substituting a@=0 in eq (21), calculating 
coefficient of ¢ and simplifying gives 


the 


[(Ao(1 + €9) — Wo) (Ao+ Ao(1 +6) k —2y K/m) 


M(1+e)(8+k+keo)]a,=0. (23) 
This equation gives a 
not vanish. 

Since d» and a, are zero, the terms in higher powers 
of t must be evaluated in order to calculate the 
effect of air drag. Thus eqs (9) and (11) become 


0 since the coefficient does 


e€=€)+ aot? +-azt 


T= T,+ Ka,.t + K,a,¢' 4 (24) 


The quantities Y and S in eqs (15) and (16) are 


recalculated to higher powers in ¢ and substituted in 
eq (20) to give 


: [(4X) +2, K/m+kW,) 


t Wi 2y K m/Xo) a,t?— F°yv6(1 


[27> 
] : ; ; 
{ 1(6A,-4 2,K m- kW) tT Eo) 
+ W,(1—2y K/m/do)]a,t?+-0t'=0. 
Putting the coefficient of ? equal to zero vields 
F?(1+-e,)*3/(4T%) 
(4, +24 K/m+kW) (1+ e 9) +Wo(1 —2y K/m/ro) 
(25) 
and putting the coefficient of ¢ equal to zero vields 
a 0, (26) 
Thus the strain in the filament is given by 
¢ + dot? 4 Of. 
Substituting this expression into eq (7) and inte- 
grating gives for the distance traveled by the trans- 


verse wave 


X=U,t+([rAki1- +-Ay—2y K/m)act®/6+-0t° (28) 


where dy is given by eq (25). 
4. Comparison With Experiment 


This theory is compared with an impact study of a 
nylon yarn by Smith et al. [12]. Their multiflash 
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Figure 4 Vultiflash photograph of a high tenacity nylon yarn struck by a projectile. 





photograph of a high tenacity nylon yarn, impacted 28 ! T 
transversely by a projectile with a velocity of 189 
m/sec is shown in figure 4. The yarn is photo- 
graphed before impact, 201 usec after impact, and 
every 50 usec thereafter. Only half of the yarn was 
photographed in order to show greater detail. The 
projectile is seen along the right-hand edge of the 
photograph. 

The curvature of the yarn in the transverse wave 
due to air drag can be seen in figure 4. The radius 
of curvature of the yarn for the first 10 exposures 
was measured and found to be approximately 
constant at 1.2 m.* 

The linear density of this nylon yarn is 92.7 tex or 
9.27 X10-° kg/m. Its tension-strain curve, measured 
at a rate of strain of 4800 percent/see by Smith et al. 
[9], is shown in figure 5. Using this curve, the strain 
é) without air drag was computed from eqs (1), (2), 
and (3). For the impact velocity V=189 m/sec, 
€o= 2.031 percent. 

In order to calculate the effect of air drag, the drag 
coefficient, /, is calculated using eq (19). The value 
of q in this equation is defined by Hoerner [11] as 


TENSION , NEWTONS 


l . 
I=5 pdC'y (29) 








where p is the density of air, d is the diameter of the 
0 3 
STRAIN ,% 





4 After the first 10 exposures, the reflected strain wave interacts with the trans- _— . 1" . . . 
verse wave and produces curvature in the yarn in addition to the curvature FIGURE ! Tension-strain curve of high tenacity nylon yarn 
produced by air drag [4] a rate of straining of 4800 percent/second. 
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varn, 0.7 mm, and (> is the dimensionless drag co- 
efficient. (Cp is a slowly varying function of the 
Reynolds number of the flow. For the case con- 
sidered, the Reynolds number is 8.8 10° and the 


| 


corresponding drag coefficient, Cp, for a cylinder is | 


1.10. 


However, the structure of a yarn affects the 
drag 


coefficient. Hoerner [11] results yield an 


average drag coefficient of 1.09 at a Reynolds number 


of 10* for two metal cables that have a structure 
similar to a textile yarn. Since the drag coefficient 
of a cylinder is 1.19 for a Reynolds number of 10*, 
the drag coefficient of the yarn is given by 


1.0$ 
i 


) 
vas 1.10=0.99. 


By eqs (29) and (19) the force acting on the yarn 
is F=11.46 newtons/m. The radius of curvature of 
the yarn, 7/F, is calculated to be 1.1 m in good agree- 
ment with the experimental value of 1.2 m. 

By eqs (25), (27), and (28) the strain, e, and dis- 
tance, .Y, transversed by the transverse wave may 
be calculated. For 601 usec after impact, the strain 
is 2.12 percent and the distance, X, traveled by the 
transverse-wave front is 19.07 cm, compared with a 
strain of 2.03 percent and an X of 19.06 cm, calcu- 
lated by neglecting air drag. Thus the effect of air 
drag on the strain in the yarn and the distance 
traveled by the transverse wave is small for these 
conditions, although it produces noticeable curvature 
in the yarn. 

The effect of air drag increases slowly with the 
velocity of impact. For example, an impact of 
300 m/sec on the same yarn produces a strain of 
4.44 percent and a distance X of 24.24 em without 
air drag, and for air drag a strain of 4.54 percent and 
distance X of 24.38 cm at 601 psec after the impact. 
Thus the effect of air drag is small for the examples 
considered. 


[1] G, 


[2] 
[3] 
[4] 
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The theory is given for 
spectral densities. The 
parator described 
and since the 


error of 
alleviates 


relative 
removes or 


a precision comparator that measures the 
a single 
many 
instrument operates under null conditions, the 


18, 1962) 


ratio of two noise 
also derived. The com- 
problems in high-speed switching, 
null position is essentially 


measurement is 


of the 


independent of amplifier noise and gain instabilities. 


1. Introduction 


The relative strength of 
signal levels is often desired. As the strength of the 
two power sources decreases the determination of 
this ratio becomes increasingly difficult. When the 
desired signals have levels that are approximately 
equal or smaller than the internal noise level of a 
very low noise amplifier, the problems become severe. 
This is due not only to the masking effect of the in- 
ternal amplifie r noise but also to the inst: ability of the 
amplifier gain. This in turn arises from the enor- 
mous amounts of gain needed. 

The usual practice is to sample the outputs of the 
two levels at a switching rate that is high compared 
to the instability drift rate of the amplifier. An 
attenuator is usually used to reduce the level of the 
stronger signal to that of the other signal. The 
early pioneer of this method was the Dicke Radiom- 
eter [Dicke, 1946]. 

The switch with its noise, and instabilities 
is often a limiting factor in such methods. The at- 
tenuator and the necessity of impedance matching 
is another problem source. These problems are 
eliminated or greatly alleviated by the method de- 
scribed herein. 


two power sources or 


losses, 


2. Principle of Operation 
The basic problem in such measurements is to 
make the amplifier gain drift essentially unimportant. 
In principle, the usual radiometer samples the two 
sources so that the variation of amplifier gain affects 
the amplified source levels of both sources equally. 
A preferable method would be to have both sources 
simultaneously present and no switching. The 
problem then, is to separate the two amplified 
sources at the output. If the two sources are 
statistically independent, they could be separated 
by correlation techniques. In principle this is 
what is done in the system described below. 

Correlation techniques are an old tool and there 
are systems described in the literature that in various 
respects resemble the one discussed here [Fink, 1959; 
Freeman, 1958, p. 274]. 
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The block diagram of the system is shown in 
figure 1. The individual components and also the 
configuration can take on different forms without 
changing the basic principles. The system described 
here is similar to one presently under construction. 





TIPLI H Hi | 
a IPLIER F-JINTEGRATOR INDICATOR 


| 





WAVEGUIDE 
BELOW CUTOFF 
ATTENUATOR 

















FiaureE 1. Block diagram of comparator. 


2.1. Sources 
The two sources may be very general. They may 
consist of random noise, pulse singals, CW, etc., or 
combinations of various types. It is assumed their 
output levels are constant during the measuring time. 


2.2. Attenuator 


For precison results the attenuator may be an 
accurate wave guide-below-cutoff (piston) attenuator. 
Small changes in attenuation may also be determined 
by measuring changes in the input level to the 
attenuator. 


2.3. Junction 


This simple but important component is shown in 
figure 2. It is a low-loss symmetrical coaxial 7. 
Opposite to where the center leg joins the 7, very 
small slits, perpendicular to the axis of the two side 
arms, are cut into the outer conductor. These are 
symmetrically placed on the 7. The junction then 
is made an integral part of the attenuator so that the 
fields (TE, mode) within the attenuator guide 
impinge upon the slits of the 7. In effect, the 7 acts 
as a center-tapped secondary of a transformer. The 
source X is applied to the center tap. 


























Fr URE 2. Junction. 


2.4. Amplifiers 


There is some value in having the two amplifiers 
similar, though this is not necessary. The amplifiers 
should be linear; at least the frequency components 
produced by the nonlinearity should be prevented 
from reaching the input to the multiplier. For some 
modes of system operation it is essential for the pass 
band characteristics of the amplifiers to remain 
constant during a measurement while the gain can 
vary. This is what primarily happens in an amplifier 
as the mutual conductance of a vacuum tube varies. 
The stability of the pass band can be enhanced by 
keeping the amplifier broadbanded except for stable 
band limiting devices. 


2.5. Multiplier, Filter, and Indicator 


The multiplier has the function of forming 
product of the input voltages. The output of 
multiplier, which is this product, is applied to 
filter. The filter permits only the d-c and very 
frequency components of this product to reach 
indicator. Of course the filter 
be a single integral unit. The indicator responds to 
both positive and negative d-c voltages and is used 
only to obtain a null. 


the 
the 
the 
low 
the 
and multiplier could 


3. General Theory 


Only the basic general theory will be presented 
This describes, under certain assumptions, what the 
system measures and the expected error (due to the 
inevitable random fluctuations in the output) in this 
measurement. Possible errors due to various com- 
ponents deviating from the assumptions will be 
reserved for a forthcoming report when a specific 
system is being described. 

The assumptions used in the following theory are: 

The system, excluding the multiplier, 
sidered linear. 

The amplifiers and filter 
considered varying in 


is con- 


* while linear 


, may be 
time. The unit 


impulse 


response and Fourier system response functions of 


these units are then time variable. The analysis 


could have also been carried through with the 
assumption that the unit impulse functions were 
sample functions from a linearly independent ergodic 
process with the essential results being the same. 

The signals present are all assumed to be sample 
functions from random processes that are bounded, 
and are singly and jointly ergodic as well as singly 
and jointly stationary. The periodic signals are 
considered to occur with a uniform random phase 
distribution over the ensemble. 

The components, other than what is stated or 
implied above, are considered to be ideal. For 
example, the multiplier performs the proper product. 
Also, the junction is symmetrical, ete. 

Sources Y and Y respectively give rise to random 
signals z(t) and y(t) that are real time varying 
functions and contain, in general, both periodic and 
nonperiodic parts. The input to amplifier A is 


x(t y(t). 
To amplifier B, it 
y(t). 


The output of amplifier A, ¢,(f), is 


a 


h(a, t){ar(t +-y(t—a) +2,(t—a)|da (3) 


a 
* 


where h,(a, t) is the unit impulse response of amplifier 
A, 2,(t) is the internal noise of amplifier A referred 
to its input, and a is the variable of integration in 
the convolution integral. Likewise the output of 
amplifier B, e,(t), is 


. 


é,(t) h,(B, t)|a(t—B) —y(t—B) +2,(¢ B)|dB (4) 


where the 6 subscripts refer to amplifier B. The 
quantities ¢,(f) and e¢,(t) are multiplied together by 
means of the multiplier. The expected value or 
statistical average of the output of the multiplier, 
indicated by E(e,é) 


BK Cold) E ( 


Or €9€, 18 


h, (a, t){r(t—a) 


t+y(t—a)+2,(t ) \de ) 


(| h,(8, t)[a(t—B) —y(t—B)+2,(t 


8) \d8 )] 


be non- 


Since h,(a,t) and h,(8,t) are considered 
random functions, 


=f . ha (a, t)h,(B, t) Ee a) 


2,(t—a)|[a2(t—B)—y(t—B)+2,(t 


ty(t—a) 


B)|)dadB (5) 


where the order of integration and averaging have 
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been interchanged. The expected value or statistical 


mean of the indicator output, /, is given by 


7 H4(0, t) | hala, t)hy(B, t) E({axct -a) 


a) +2,(t B) 


a) |[a(t 


y(t—B) + Z(t -8)|)dadB (6) 
where /7,(0, t) is the zero frequency response of the 
filter and indicator. 

The expected value of the time varying signal 
function in (6) is composed of the statistical auto- 
correlation and cross correlation functions of the 
various random signals present. Thus 

E{xr(t B)| R.(8 (7a) 


a)r(t a) 


the statistical 
Similarly 


is 


where R,(B—a) 
function of x(t). 

E\x(t a)y(t—B)| R.,(B—a) (7b) 
with similar results for other combinations, where 
R,,(8—a) is the cross correlation function of the 
functions z(t) and y(t). 


From the foregoing, (6) may be written 


hala, t)hy(B, t)[R(B—a@) 


I H,(0, t) ( 


R,,(B 


a)-+ R,.,(B—a) 
a) R, (8 —«a) 


+ Rz2,(8—a)|dadB. (8) 

If the respective signals are uncorrelated all cross 
correlation functions become zero, and (8) reduces 
to (correlation between z,(t) and 2,(t) is discussed 
later) 


I HT ,(0, t) hala, t)h,(B, t) 


[?,(8 a) R,(B 


a)|dadB. (9) 

The above correlation functions are related to the 
spectral densities of the various signals. Let S,(f) 
represent the spectral density of the signal z(t). 
Then S,(f) and RF,(r) are related (by definition) 
through the Fourier transform pair 


»@ 


S,(f) | R,(r)e~ 4dr 


« 


(10a) 


R,(1r)= S,(f)e?*df (10b) 


with similar relations for the y signal. 


| 


| 


Substituting the above values for the correlation 
functions into (9) together with the relation 7==8—a 
vields 


/ H,(0, t) | ( ( ha(a, t)e~?*/“h,(B, t) e?*# 


[S.(f)—S,(f) |ldfdadg. (11) 


Is 


Now the impulse response function, h,(a, t), 
related to the complex frequency response function, 


| H,(j22xf, t) by 


ha(a, t)e~#*da=Hq(j2xf, t)=\Ha(j2nf, t)\e% 


rs) 


(12) 


| where H,(j2z/,t) is the complex frequency gain 


autocorrelation 


function of amplifier A and @, is the associated 
phase angle. Similarly 

Ao 

| hy(, t) e*8d8—H* (j2xf, t)=|Hy(j2nf, t)|e~» 

(13) 

In general, @, and 6, are functions of time, though 
for simplicity in notation they are not explicitly 
written as such. Thus (11) becomes 


7 H,(0, t) H,(j2nf, t)||Hy(j2xf, t)\[cos (6.—6,) 





| The total power, 





+7 sin (O.—6»)|[S:(f)—S,(f)|df. (14) 
[It may be shown that H(j2z/f,t) has an imaginary 
part that is an odd function of frequency. 
Using this fact, (14) becomes 


H,( j2nf, t)||Hy(j2xf, t)| [cos (@.—6 


7=H,(0,t) q 


— 2 


[IS.(f)—S,(/))] df. (15) 


Periodic signals such as CW signals or even 


| periodic wide-sense random processes have power 


spectral densities expressed by means of the Dirac 
delta function, 6[Davenport, 1958; ch. 6]. For 
example, a CW signal, V cos (wof+6) has a power 
spectral density 


therefore, is 
RPE 
Scw( fydf==- 


vy -~@ & 


as is well known. 
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The output, /, would become zero (a null would | 
exist) when 


With such amplifiers a time-stable null condition 
expressed by (22) may be met with signals having 


| spectral densities with widely different frequency 


HT ,( j2nf, t)||Hy(j2xf, t) cos (@.—4)S,(f) df 


a 


H,( j2xf, t) H,(j2nf, t)\cos (0a—6) S,(f) df. 
(18) 


It is of considerable practical importance to 
consider those cases for which a null is obtained 
independent of the variations in the amplifiers and 
filter response functions. Since (18) contains no 
filter response function, the null is independent of 
the filter characteristics. From (18) it is also 
obvious that if the spectral densities of the x(t) and 
y(t) signals were equal over the amplifier pass bands, 
that is, if 

S,(f)=S,(/), (19) 
a null is then obtained regardless of the time varying 
aspects of the amplifiers. Thus the important | 
result is obtained that a drift or variation in the | 
gain or phase response of the amplifiers and filter 
would not influence the null condition. 

Further, if the amplifiers were constructed so that 
the frequency response determining components 
were stable, the amplifier characteristics could be 
written 
H,(j2nf, t)||Hy(j2zf, t) 6») 


COS (6g 


G,(t)G,(t) H,(j2rf) H),(j2rf) cos (8,—6,) (20) 
where G,(t) and G,(t) represent the time variable 
or unstable parts of the gains of the respective 
amplifiers and the primed H’s are new time-stable 
frequency-dependent parts of the amplifier gains. 
The @’s are now also independent of time. 

The output of the system may then be written, 
since G,(t)G@,(t) is independent of f over the bandpass, 


2 


T=H,(0,t)G,()G,(t Hi, (j2mf)||H',(j2xf) 


cos (6, 6,)[S,(f) S,(/)\df. (2] ) 


This becomes zero when 


| , j2rf) H,(j2rf) cos (6, 6,) S.(f)df 


> 


HH (j2nf) 


H,(j2rf) cos (6, 4,)S,(fdf (22) 


and is independent of the amplifier drift or variation 
in gain. Such amplifiers are not too hard to realize 
by making the unstable frequency-dependent parts 
with pass bands broad compared to the stable 
frequency-dependent parts. 


characteristics. For example y(t) may be a sinusoidal 


| signal 


y(t)=Y 9 cos (wot +8). 


Then (22) becomes 


( HT), (j2xf)||H,(j2xf)| cos (6,—8,) S( fydf 


(|, (j2nf)\|Hi(j2xf)| cos (@—8)] 


* 


ro 


athe i pate <a ; 
H1),(j2mfo)||\H,(j2mfo)| cos [O.( fo) —O(fo)] = (24) 


This is again independent of the amplifier gain 
variation. 


4. Measurement of Spectral Densities 


From the foregoing relationships, spectral densities 
of one source may be obtained in terms‘of that of 
another source. 

Assume, for example, the value S,(f).is known 
and is also known to be essentially constant over 
the pass band, and it is desired to obtain the value 
of S,(f) which is also known to be constant over the 
pass band. 

With reference to (18), 
when S,(f)=S,(f). 
thus equated. 

Let two different z(t) signals be measured. 


a null condition is obtained 
The two spectral densities are 


Then 


The ratio 
Sio(7) 
Sil J) 


may be measured by means of the precision attenua- 
tor (see fig. 1), and hence 
S(f)=AS2(f). (28) 

In this case S,,(f) and S,.(f) need not be equal 
and only relative values of S,(f) are needed. The 
relative values of S,(f) may be obtained by means 


| of the attenuator. 


A CW signal could just as well be used in place 
of S,(f). Then 


5 Sa(S) 


Ti AS;,(f) (29) 
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where 
V 
A 2 


and is again measured by the attenuator. 

If S,.(f) and S,,(/) vary differently with frequency, 
the relationship between them may be obtained by 
using (22) to relate S,,(f) with S,,(f) and S,.(f) with 
S,(f). Noting also that S,.(f)=AS,,(f), the desired 
relationship becomes 


|, (g2af)||\H,(g2xf)| cos (0,—8,)Sro(fdf 


A H!,(j2nf) 


H,(j2nf)| cos (0.—9) Sin(f)df. 


(30) 


If the spectral density of S,,(/) is constant in value 


over the pass band, (30) becomes 


| Hi, (j2af)||\H,(j2xf)| cos (0,.—4)So(fydf 


H,(j2nf)\\H( j2xf)| cos (0,.—0,)af 


AS,,(f). (31) 

In practice, this is what is usually measured. 
This gives the equivalent constant spectral density 
that would vield the same total power as the actual 
spectral density function when both are integrated 
over the actual pass band response function. This 
means the relative pass band response function of the 
actual system used should be given when stating the 
equivalent S,(/). 

A system that meets the conditions imposed by 
mad lends itself very well to a measurement of the 

‘elative pass band response function, 

K\ Hi (j2af)\|H,(j2xf)| cos (0.—4). 

Here A is a constant and represents the fact that 
only relative and not absolute values are needed as 
a function of frequency. The above funetion is 
obtained by using CW signals for both the z(t) and 
y(t) signals. The y(t) signal may be varied in level 
to maintain a null as z(t) is varied in frequency but 
kept at a constant level. The relative values of the 
attenuator setting as a function of frequency give 
the desired relationship. 

Other ways of using the system may r 
to mind. One other will be briefly mentioned. 
Another uncorrelated signal o(t) with spectral 
density S,(f) may be combined with the z(¢) signal. 
The relative amounts of this extra signal may be 
measured by an attenuator external to the z(t) signal 
in much the same way as is done with the y(t) signal. 
In this case, the y(f) signal need not be varied ‘and, 
(19) could be written 


radily come 


S,(f) =Sa(f)+Sa(/) 


So(f)+S8. Si3(f)+S8(f). (32) 





pare various random z(t) signals. 








If only relative values of S,(f) are 
relation is 


known a useful 


Si3(f) 
Sul(f) 


where 


4.1. Broad Spectral Signals Versus CW Signals 


As indicated, either broad spectral signals or 
a CW signal may be used as the y(t) signal to com- 
When CW signals 
are used, it is important that the frequency of the 
CW signal be constant. Also, practical 
dictates finite averaging times at the output. This 
means that a CW y(t) signal would produce less 
random fluctuation than a random y(¢) signal. 

Broad spectral y(t) signals, however, place less 


usage 


| stringent requirements on the constancy of the band 
| pass 


characteristics, especially if 
densities are similar to those of z(t). A broad 
spectral nonrandom signal would produce less 
random fluctuation than a random signal and still 
have those advantages of broad spectral signals. 


their spectral 


4.2. Residual Noise 


If the signals z,(t) and 2,(t) are correlated an 


| additional term 


H,(j2nf,t)||\Hy(j2af,t) ||S.2,(/) 


Za 


cos (6é,—6,-4 


would have is the 


occurred. Here, ¢, 
phase of the correlated part of the amplifier noise 
signal in (15). This represents a contribution in the 
output due to correlation between the two amplifier 


noise signals. While there is usually very little 


2b relative 


| correlation between the two amplifier noise signals 


there can be contributions by the two amplifier 
noises that have common origins and are therefore 
correlated. Whether the effect of this correlated 
residual noise can be ignored depends on the relative 
strength of signals being measured to the residual 
noise and the desired accuracy of measurement. 

This residual noise may possibly be reduced by 
proper design of the amplifier, or by using special 
networks that prevent coupling or that cause shifts 
in the phase of the correlated a of the residual 
noise such that the quantity @,—+¢,,:, in (17) 
equals 7/2, thus making the cosine factor zero. 

Also, the residual noise may be simply measured 
by the use of two sources of different known levels. 
The principle may be simply demonstrated by 
assuming (this is not a necessary assumption) all 
signals have constant spectral densities within the 
amplifier pass bands and the amplifiers have equal 
phase shifts (@,=6,). Using (15) with the 


assump- 
tion that no signals are correlated except 


z,(t) and 


327 





2,(#) and that the null condition holds, then 


S,(f) —S,(f) +|S.,2,(f)| cos $2.2, 0. 


For two different levels of S,(/), this vields 


where 
Sp(/) 
Su(/) 


5. Relative Error of a Single Measurement 


Any practical system will average the output over 
a finite time. Because this time is finite, random 
fluctuations will appear in the output giving rise to 
errors when the output indicator is read. 

For simplicity in evaluating this error, the follow- 
ing assumptions in addition to those above (with one 
exception) will be made. 

1. Each amplifier has a 
width B. 

2. The random processes giving rise to z(t), 2,(t), 
and z,(f) are considered to be Gaussian, and the 
respective signals are independent and thus uncor- 
related. Also, the various signal sources have zero 
means. 

3. The signal z(t) has a spectral density S,(f) 
that is constant over the band. 

4. The signal, y()=Y, (wot), Is now a de- 
terministic (nonrandom) CW signal centered in the 
pass band of the amplifier. The former theory still 
applies. The present system merely represents 
sample functions of the ensemble with time origins 
fixed relative to the CW signal. 

5. The two amplifier noise signals, z,(t) and z,(f), 
are noncorrelated. ‘Their respective spectral densi- 
ties S. (f) and S.,,(f) are constant over the band. 


square band pass of 


cos 


6. The exception to the former assumptions is 
that the response functions of the amplifiers and 
filter are constant with time. 

It is quite common to represent a narrow band of 
noise of band width B around a central frequency by 
[Davenport, p. 158] 


X (t)=X.(t 
Rit 


X ,(t) sin wot 


COS Wol 


g(t) (36) 


COS | wWol 


where ,(f), Y.(#), R(t). and (tf) are random fluc- 
tuations that are slowly varying (due to the narrow 
band) with respect to the central frequency, fo 
w,/2mr. Similarly 


t)=Z,(t 
t)=Z,.(t 


Zautt 
Zyx(t) SIN wot. 


COS Wot SIN Wy 


COs Wot 


It can be shown |[Davenport, p. 158] that 


> 


S,(f)df=2S,(0) B 


Z:, 28. (O)B 





where S,(0) denotes the constancy of S, with 
quency. 
Also, it can be shown [van der Ziel, 1954] that 





X?2(t).X2(t+7) (Xx? 


2h? (T) 7 


where 


X.()AAt+1) Sipe wd} 
sin rBr 


2S8.(0)B B 
TOT 


(44) 


since S,(0) is constant. Similar relationships hold 
for X,(t), Z.(t), and Z,(#). 
[It will also be assumed that the two amplifier 
noise signals are of equal strength. Hence 
S, (0)=S,,(0) =S,(0). (45) 
The inputs to the two amplifiers will be 
Input a 


+-Z,.(t)| Sin wet (46) 


Input 
ep(t) =[.X- “ot+Zre(t 


| COS Wot 


X(t) +Zp,(t)] sin wot. (47) 

For mathematical convenience, the filter is con- 
sidered to consist of two parts; (1) a zonal filter 
(i.e., a filter that function that is 
unity over the pass band and zero elsewhere) that 
effectively integrates [Freeman, p. 226] the output 
of the multiplier over a period long compared to 
that of the central frequency fy, and (2) a filter that 
is effective at much lower frequencies. 

The output of the multiplier, M, is effectively 
time averaged over a cycle of the center frequency, 
wy, by the zonal filter. Thus 


has a system 


M(t) é,(t)e,(t) 


)+-Z,.(t)| |X,(t)+-Z,,(t)]} (48) 
where <e,(t)e,(t) > represents the temporal average. 
To obtain the spectral density of M, its auto- 
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Now 


correlation function is needed. 


Ry (7) M(t) M(t T). (49) 


If the value of M(#) and M(t+-7) are used as given 


in (48) together with other relationships expressed 
above, it may be shown that 


| i? wai Br [1682(0) 
4 (rBr)? 


16S8,(0)S.(0) 


Riyy(r) 


; sl 3 
t8§S?(0)|+-4Y2S.(0)B sin wt id 
; rBr 


+ 1682(0) B?—sY2S,(0) BA rs}. (50) 


Since PR,,(7r) is an even function of 7 
oe 
Salt) =2 


J0 


Ry (r) cos (2xfr)dr. 


(51) 


Substituting the value of Ry(r) as given by (50) 
and integrating 


l ; 
Su() {145.(0) B YP) 


f Y-S.(0) +2/282(0) +2S8,(0)S2(0) 


.$2(0) ][ B— |] 


> 


for 0< | fi<s5 


2(2.82 (0) +2S8,(0)S8,(0) +S3(0)|[ B—|f\] 
ae 
lors 





L 0) for fi >B. (52) 


The mean value of the indicator deflection, /, is | 


given by 


ii : H,(0)|4S,(0) B—Y?| (53) 


which may be obtained from (48) by taking the aver- 
age over the ensemble and multiplying by the d-c 
response function /7,(0). The same value can be 
obtained by integrating (52) over the vicinity of 
zero frequency, taking the square root, and multi- 
plying bv /7,(0). 


occurs when 


At null, of course, / is zero which 


y?—48S,(0)B. (54) 

In adjusting the output to a null, the fluctuations 
in the output cause an uncertainty in the value of 
Y, for a null condition. This means that the 
measured value of ¥) and consequently the meas- 
ured value of S,(0) are random variables. These 
are designated Y>,, and S,,,(0) respectively and are 
related through (54); that is 


245. (0) 8: (55) 
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The relative error, ¢, in a single measurement is 
now defined as 


Vv (Sm(0)— Srm(0) )* 
Srm(0) 


> 


VAS tm (0) 


S,(0) 
where 
S-m(0) — Sam (0) 


“rm 


AS,(0) 


and 


S-m(0)=S,(0). (57) 


It is assumed the operator (or servomechanism) in 


| adjusting the attenuator makes an error, AY, (which 
| is a random variable) that just compensates for the 
| fluctuation, AJ, in J(t), where 


This means, using (53), that 
sa las ace ee 
Al=—; H,(0)(48,(0) B—(Yo+ AY)’}. 


Simplifying by the use of (54) 


AI=H,(0)Y AY (60) 


| where AY? has been neglected compared to AY Y5. 


Using (58) it can be shown that 


AP | H,( j2nf)|*Su( fydf—(1)?. (61) 


| When /7,(j2zf) is nonzero except at zero and near- 


zero frequencies, it follows that essentially 
Al’: S(0) ( H,(j2nf) "df. 
Now, from (52), and using (54) 


S4,(0) 481(0)B| 142747 | 


where 


S,(0) 4 
Y=<sa (604) 
S,(0) 
and the prime denotes S,,(0) less the term involving 
the delta function. 

The quantity 1+¥ is equivalent to the operating 
noise figure. This would reduce to the standard 
noise figure when z(t) is equivalent to a signal that 
originates from a resistor at a temperature of 290 °K. 

Using various relations as given above, the follow- 
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ing chain of relations may be obtained 


4A]? 
Hi(0)¥5 


4 F De f)\20 \ Af T)2 
avert |. Lax f)Su( PO 1] 


ae oh 14.24 +¥) ( 


VS 
H,( j2nf)\2 
H,,(0) 


AS?2(0) 4AY? 
S2(0) r: 


Hy q2nf) : 


H,(0) | 


df. (65) 


For a simple RC filter of time constant 7, the 


integral is equal to For a critically damped 


l 
2T 
system the integral would be ' that given above. 

Using the simple RC filter, the equation for « 
becomes 


(66) 
For example, if 
B=5x<«10° ¢ 
EA 
i] 3.5 see 
then 
10~* or 1 percent. 
If y=10, then 7 must be increased to SO 


obtain the same value of «. If y=0, then 
only be 1.0 see for the same value of «. 


sec to 


T need 


6. Conclusions 


The S\ stem described has several desirable char- 
acteristics. 

1. The normally used high-speed switch has been 
eliminated. This removes the problems of switch 
noise, insertion loss, and instability. 

2. The effect of amplifier gain drift 
essentially eliminated. 

3. The system is operated under null conditions 
and the indicator null position is therefore inde- 
pendent of amplifier gain settings. 

4. A precision, continuously variable piston at- 
tenuator may be used. 


has been 


5. Signals of widely varying strengths may be 
compared. 

6. The system is capable of measuring its own 
effective pass band response to a high degree of 
precision. 

7. Problems of impedance variations in a dissi- 
pative attenuator have been removed. 

8. The system is quite insensitive to “hum” and 
other such signals that might possibly modulate the 
normally present signals. This follows since the null 
condition is independent of a time variable system 
response. This time variation if produced by “hum” 
or other signals is what often gives troubles in 
systems, especially those using linear or envelope 
detectors and switching frequencies commensurate 
with the undesirable modulating frequencies. 

Some undesirable aspects are: 

1. One must have two low noise amplifiers. 

2. The signal power is now divided into 
channels. 

3. A determination of the effect of correlation be- 
tween the internal amplifiers may be necessary in 
some situations. This may be done by comparing 
two standard sources of different levels. 

4. An appropriate multiplier is needed. Only the 
output at zero frequency is of interest and this fact 
removes many problems. Also, an ideal multiplier 
is not necessarily needed. The main consideration 
is that no output at zero frequency should exist 
when only one input is energized 


two 
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distributions in 
central portions of two opposite edges are 


The s ress 


rectangular plates 


evaluated theoretically. 


and Wilhelmina D. Kroll 


1962) 


with loads over 


shown to be 


uniform 
Results are 


symmetrical 


applicable to the stress analysis of a uniaxially loaded plate which is elastically constrained 


by an unloaded plate perpe ndicular to it. 
welded assemblies. 


A comparison of experimental results with theoretical results for a simple 


structure shows reasonable agreement. 


1. Introduction 


Following the occurrence of serious structural 
failures in welded ships during World War II, 
investigations were undertaken to determine the 
cause and possible remedies for these failures. It 
was apparent to early investigators that at least part 
of the failures in large monolithic structures could be 
attributed to the brittle behavior of material com- 
monly considered ductile. Experimental studies 
{1, 2}' were conducted to confirm this assumption. 
The specimens used in these studies, as shown in 
figure 1 which is taken from reference [1], represented 
interrupted bulkhead intersections in welded tankers. 


Information from these and other studies led to 


modifications of existing ships and improved design 


of new vessels, 
number o 

Strain these specimens also 
revealed the existence of a complex stress distribution 
near the intersections of two plates in the test 
specimens. No theoretical analysis was found to 
explain the distribution of stresses in the test 
specimens. This may be due to the fact that a 
rigorous investig even a simple monolithic 


and resulted in a substantially lower 
f failures [3]. 
measurements in 


ration of 
structure, such as shown in figure 2, in which the base 
plate under uniaxial tension is elastically constrained 
by two flanges, is a three-dimensional problem with 
mixed boundary conditions and is difficult to solve. 

However, it is shown herein that the problem of 
determining the stress distribution in the structure of 
figure 2, consisting of plates whose thickness is small 
compared to other dimensions, can be reduced to the 
stress problem of a plate loaded over portions of its 
boundaries. The latter problem has been adequately 
treated in standard texts but numerical results are 
available only for the cases of narrow strips under 
concentrated loads [{4, 5]. This paper gives the nu- 
merical results of stresses due to a uniform load over 
central portions of two opposite edges of plates of 


*Present addres lhe Catholic University of America, Washington, D.C, 
1 Figures in brackets indicate the literature references on page 346. 
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This structural configuration is often found in 


welded 


various dimensions and indicates how the results may 
be used to estimate the stress distribution near the 
intersection of two perpendicular thin plates of 
arbitrary lengths. 


2. Analysis 


The simplified welded structure shown in figure 2 
can be considered to be equivalent to the plate shown 
in figure 3 since, in these figures, the shaded portions 
which correspond to each other have the same re- 
straining effect on the respective base plates. The 


| equivalence of these two problems can be verified as 


follows. If the plate shown in figure 3 is thin, the 
shaded portions at the ends can be bent through 90 
degrees without changing the stress picture. Then 
the plate can be cut at the middle along the y-axis. 
The continuity effect is replaced by boundary forces 
which are equal to p. The two halves of the plate 
may then be joined at the bent ends with the shaded 
portions in the up and down positions shown in figure 
2. Note that the deformations along the bent edges 
are compatible because of symmetry and that the 
loads p cancel each other when the halves of the plate 
are jomed. The boundary forces produced by cutting 
then become the end forces shown in figure 2; 
these two figures are equivalent to each other. 

The analysis of the plate in figure 3, therefore, will 
form the basis of evaluating the stresses in the welded 
structure of figure 2. 


hence 


A direct solution of the problem shown in figure 3 
is rather difficult. However, the problem becomes 
almost trivial through the use of the superposition 
principle. It will be shown in the next few para- 
graphs that the problem shown in figure 3 is statically 
equivalent to the combination of cases shown in 
figure 4b and 4e. 

First, it is assumed that the joints between the 
shaded and unshaded portions of the plate shown in 
figure 3 are cut. A uniform tensile load would pro- 
duce rigid body displacements of the shaded portions 
of the plate without producing stresses in them. 
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For the central part of the plate of figure 3, which 3. Theory 
is under the unit tensile load p, the stresses and 
strains in the transverse direction of the plate | In order to determine the stresses in a plate such 
would be: as shown in figure 3, use has been made of existing 
and theoretical work. The expressions for the stresses 
in a plate, figure 5, with a symmetrical load distrib- 
ution on two of the boundaries representable by a 
Fourier cosine series have been developed in texts 
where on elasticity. See, for example, [6]. 
o,=stress in y-direction 
€,—strain in y-direction 

load per unit area 
Young’s modulus 
Poisson’s ratio. 


The deformed configuration of the plate is shown WY 


in figure 4a. fj VA 


To prevent the contraction corresponding to e¢,, a 
uniformly distributed load g must act along the 


edges of the central portion of the plate. See figure (a) Contraction of bose plate under tensile forces 
4b. The magnitude of q is then, 


q 
qd = — Key=vp. Ltt tee e ttt ttt ttt 


Due to the uniform load g, the deformations in 
the shaded and unshaded portions of the plate are 
compatible. The portions of the plate can now be 


joined to restore continuity. However, to be con- RERRERIREAE ALL 


sistent with the condition of a free boundary, a (b) Uniform load added to prevent contraction of 
load —q must be applied to the edges of the central port of plate under tensile forces 

portion of the plate to neutralize the edge load gq 
as shown in figure 4c. It is clear that the effect of “qd y 

loads p on the structure shown in figure 3 is equiva- Sod ddd ddd ddd dds 
lent to the effect of a system of forces shown in 7777 

figure 4b superimposed on another system of forces 
shown in figure 4c. 
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(c) Uniform load added to obtain system of forces 
statically equivalent fo original one 


= Figure 4. Method of analyzing welded joint. 
ad i. = 











FiaurE 2, Simple welded structure. 
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Fiaure 3. Structure considered equivalent to welded structure. Fiaure 5. 








Dimensions of plate under uniform load q. 
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The equations for the stresses are: 


mre 
(7 


mre . mare mary 
cosh j —sinh j -) cosh v. 


a 
° MTree m7re 
sinh 2 ——+-2 
alle l 
Mine 


] cosh 


mre, . mare mary 
+sinh -) cosh - 
l l l 
Mme . 9 mre 
—+.2 
l l 


sinh 2 


mare mre mary mry 
cosh —— : 


j j sinh j j 


Wiry 
CC sh 


° Mine mre 
sinh 2 +2 
l l 
where 


half length of loaded portion of plate; 
half width of plate; 
half length of plate; 
',y=coordinates of points in plate, meas- 
ured from the mid-length 
mid-width of the plate, 
tively; 
m=O, I, : 
etc. 


The coefficient 


and 
respec- 


) 


Ao, Aj, 


Fourier coefficients. 


If a normal load g is uniformly distributed over the 
central portion 2a of the plate, the Fourier coeffi- 
cients for m=1, 2, 3 © are [7]: 


4.2 
j 


1 . 
q COS 


For o,, 
Wind 


da 
l 


mra 


l 


Sin 


Substituting the values of the coefficients and the equivalent expression 


the equations for the stresses become: 


g mare mare ° mre 
| ( cosh ———sinh 
I l 


l l 


° Mm7TCe mr7re 
sinh - cosh 


l l 


mre mane 
( “ cosh 
Mirada MmaTL / / 


“4 n7re 

| T sinh : - ) cosh 
l 1 

\ 


sin COS 


¢ mre mre . mary 
—— cosh sinh -— 
mrr l / 


lets 


mary 
— Sin 


° mare mare mane 
sinh / cosh T 


q el 


unless modified as indicated in the following sections. 


mary . 


mry 


j sinh 


; ) mary 
cosh — 
l 


° Mm7re 
sinh cosh 
l 


h —= 


mary 
| 


mre 
sinn “ 


l 


sinh 


mary 


ma7Ac 
sinn S s 
l 


inh j 


mare 


“ mre WTC » 
2 sinh —— cosh - for sinh 2 


l l 


mal 


y. mre . may) 

j sinh j sinh j 

mae 

l 

mary mary mre. mary 
cre sinh >— sinh ~ 


ih so ma7re 


+ 


[a 
mary) 
l 


cosh 
(4) 


J 


The above equations were used in computing the stresses in plates for various dimensions and loadings 
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4. Computations 


The computations for the stresses in plates with 
discontinuous loadings were carried out in an IBM 
704 electronic calculator. The dimensions of the 
plates and the parts of the plates under uniform 
load are.given in table 1. 

The values of the stresses using eq (4) converge 
rather slowly; a large number of terms are usually 
necessary to assure sufficient accuracy in the results. 
Based on hand computation for the stresses at 
particular points in the plate, the following procedure 
was used to obtain stresses to five decimal places. 
At first, the stresses were computed for values of m 


mrad 


for which 0 2r; then the stresses were com- 


mra 
l 


increased by 27, a compari- 


puted for values of m for which 0 4r andsoon. 


» . mrad 
Each time the angle 


l 
son was made between the value of o, for this cycle and 
the preceding one. If the value of o, changed by 


0.00001 or less, the computation was stopped. It 


f 
mane | 
Os 
' 1 
\ 


mre 
2g . mra l 
sin c 


mr l 


TABLE 1, 


Loaded Nonloaded 

length length of 

of plate plate, 
2a 2(l-a 


Loaded length 
Width of 
plate Total length 
2c a 


l 


This reduces to 


; Gis 2q 
l m=1 tw 


However, from the loading shown on figure 5 


’ 


, Mire . mare ° i 
cosh? j +sinh > cosh — sinh? 


Loaded length 


was thought that this criterion was adequate but 
in several instances the computation was stopped 
before the series had converged. Therefore except 
on the boundary, 100 terms in the series were used 
. , mrm(e—y) . 
for all points, as long as : 16. 


l 


The computa- 


mr\c—y) 
l 

greater than 16 since the increment in the stresses 

Mme _ mary 

T - O! ry 
greater than 8, the hyperbolic sine was considered 
to be equal to the hyperbolic cosine with no de- 
tectable error within the limits of accuracy desired. 
The terms under the summation sign in eq (4) were 
accordingly modified as the summing proceeded to 
avoid overflow in the computer when computing 
the hyperbolic functions of very large angles. 

Along the boundaries, y= +c, the values of the 
stresses showed no signs of convergence, even when 
as many as 300 terms were taken, but became 
oscillatory. The equation for o, at y=e is from 
eq (4): 


tion was stopped if the quantity became 


would be essentially zero. For angles 


mare mane mane 


= 

l l | * 

: mre mane mre a (9) 
sinh cosh +—- | 
J 


l l l 


Dimensions of plates investigated and computed maximum and minimum stresses 


Nonloaded length Maximum and minimum 


stresses 


Width Width 


l-a 
c 
0. 040 
—1, 025 


0. 099 
1. 023 


0. 176 
—1,023 


0. 126 
—1. 038 


0. 023 
—1.025 


0. 022 


—1. 044 


0. 022 
—1.024 


0.018 
—1. 001 


-—g for 


:() forasrsil. 





When the stresses in the y-direction are known, for 
example, and those in the z-direction are expressed | 
by a slowly convergent series, the stresses o, at the | 
boundary can be computed by making use of the 
following relationship [s]. 


0; : | — fro : 


The equations for ¢, and a,, eq (4) 
eq (8) give 


(S) 


, substituted into 


mra Marr 


. / cos jj 


sin 


. mre mre, mre 
sinh —— cosh ——+-—— 


l l l 


yore 
/ 


(9) 


The last term is given by eq (7). The series part of 
eq (9) is rapidly convergent and only 20 terms were 
needed in the summation as long as the value of 
mare , ~ mre 
7 Was less than 9. When the value of j 

became greater than 9, the computation was stopped 


(es cosh es toh we 
c c Cc 


, 2 
sinh 


2A 


‘ 2r 
sinh 
- 


al rl . wd I 
cosh — sinh —— 
“a2 - - 
2. = = 
sinh a PE 


( c 


where A=—0.017q and 


—(.055q for cases 1 and 2, 
respectively. 


The stress residuals for case 3 
parabola, figure 6c. To satisfy the boundary con- 
dition that the stresses in the z-direction be zero 
at the boundaries r= -+/, the approximate solution 
for the stresses when the forces on the ends are 
distributed according to a parabolic law [10] was 
used. If, as it is pointed out, the stress function 
is in the form of 


approximated a 


] 


9 


> ] y > \9 2 9\9 
o=5Py(1—-4 )+ (2?—a?)*(y’—e)? 


6c" 


(ay +-aox?+a,y") (11) 
and the plate is twice as long as it is wide, the stress 
distribution across the plate at z=0 is almost uni- 


form. The equations derived for the stresses with 


Td ae 
) cosh —— sinh 
Cc c 
rl rl e rl TI TI ° rl 
( cosh ——sinh ) cosh —— sinh 
? c c Cc c “4 


a rl Td 
sinh - cosh 
Cc . 


as any error in not taking more terms in the series 
would be insignificant. 


5. Correction of Boundary Values 


An inspection of the computed values for cases 
1 to 4 showed that eq (4) satisfied the boundary 
condition that the shearing stresses at the boundary 
z= +l be zero but did not satisfy the condition that 
the stresses normal to the boundary be zero. The 
stress residuals in the z-direction were small for case 
1 but became progressively larger as the extent of 
the loading approached the boundary of the plate. 

For cases 1 and 2, table 1, the distribution of the 
stress residuals followed approximately a cosine 
curve (see the solid curve on fig. 6a and 6b). 
For case 3, the distribution of the stress residuals 
approximated a parabola, figure 6c, and for case 4, 
a combination of a parabola and cosine curve, 
figure 6d. For cases 5 to 8 where the loading is on 
a relatively small area and distant from the plate 
boundary, there were no stress residuals in the 
z-direction. 

Te +-/ for cases 1 
and 2 and so satisfy the boundary conditions, stresses 
were computed by use of the following equations [9], 
and added to those previously computed 


eliminate the stresses o, at z 


wr . rl 
sinh 
Fe 


tls rl 
c c 


° TI 
sinh 
= 


as 


rl 
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Cc . Ty 
sin ’ 
c 


this stress function were used, therefore, to compute 
the corrections to the stresses. The equations are: 


0.37 (i—4)+4(; 1) (34-1)(0.02986 


a 
2 rp 
+0.04675 %+-0.00683 4 )—0.2455 
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Stress residual 


Correction 
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(ce) Cose3 (d) Case 4 
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Ficure 6. Stress residuals on transverse boundary of plates fo 
cases 1 to 4. 
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FicuRE 7a. Distribution of oy/q in plate for case 4. 
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FiGuRE 7b. Distribution of o,/q in plate for case 4. 
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FiGurRE 7c. Distribution of rxy/q in plate for case 4. 
The last term in the equation for o,/q is a constant 
term corresponding to a uniform load. 

The equation for the distribution of stresses in 
the z-direction at the boundary of the plate r=+1 
for case 4, figure 4d, is 
| f Da 


_977(1—¥ , 37Y_ 9 range 
m0 (1 r) 0.08 cos 20 0.53048, 
(13) 


oz/q 
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FIGURE 8a. stribution of oy/q an plate for case 8. 
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FIGURE 8b. Dist 
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FiGURE 8c. Distribution of Txy/q tn plate for case 5. 

which is a combination of a parabolic distribution 
and a cosine distribution. The corresponding 
stresses in other regions of the plate were obtained 
by superposing the stresses due to the cosine dis- 
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FIGURE 9b. Distribution of o,/q in plate for case 8. 


tribution on those corresponding to the parabolic 
distribution as obtained from solutions similar to 
those for cases 1 and 2 and case 3, respectively. 
These stresses and those corresponding to a uniform 
load, o,/g=—0.53048, were added to the computed 
stresses to give essentially zero stress on the bound- 
ary at gz tl, 





Rotio of distance from 6. Results 
longitudinal ¢ to width, = 

Values of the ratios of stress to applied load 
per unit area, o;/q, o,/q, and 7,,/q, were computed 
so that the results would be applicable for any 
value of load g. In the remainder of the paper, 
these ratios will be referred to as stresses. 


‘4 


Typical plots of the stresses oe 2, ond 2 ee 
q q ( 
corrected are shown in figures 7, 8, and 9. The 
maximum and minimum stresses are listed in 
table 1. 

The computed values of stresses are given in 
tables 2 to 9. These tables are presented in terms 
of nondimensional parameters for generality. Since 
the stress distribution in the plate is affected only 
in the region of load discontinuity, one of the 
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parameters used is - where (rx—a) is the distance 
P 


measured from the point of discontinuity of the 
thution of rxy/q in plate for case 8. load p in firure 4c or figure 3. 


TABLE 2. Stresses in plate for case 1—Continued 


704 Resul Corrects 


Corrected 


0, O83 
O80 
O64 
020 
OR7 


OR1 
074 
016 57 .012 | — 050 
030 : 2 0 
O38 0. O82 
031 180 0 


064 012 0 
O56 013 —(). O16 
033 013 030 
OO8 .010 . 038 
004 031 

0 


—(), 029 026 0 
025 024 —(). 010 
oll 019 —. 017 
007 O11 019 
026 004 014 
036 0 


017 023 0 
O14 021 0 
005 016 0 
005 009 0 
014 003 0 
O15 0 0 
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TABLE 3. Stresses in plate for case TaBLE 4. Stresses in plate for case? 


Unloaded length l—a Unloaded length l—a 
a =15:l>> ¢ ’ =1.0; 1>>¢ 
Width oo Width ; 7 


No. of tesults Corrected | 704 Results Corrected 


terms, 
m 


y/9 Tzy/Q 


0. 002 - | —2.5 0. 002 
001 . : : 001 
U | f _ 0 
—(. 001 7 R —), 001 
001 | - 001 
0 . g 001 
001 
00S 
007 
003 
—. 003 
007 
OOS 


. 028 
024 
011 
006 
025 


036 


0. 065 

016 057 
030 033 
038 —. 008 
031 —. 061 
-114 


O82 
075 
050 

0 
—(). 083 


0 


0 
0.010 
022 
036 
039 
037 
024 
0 


0 


OS] 
037 
001 
. 005 


029 
024 
010 
004 
O09 
007 


016 
017 
O18 
O17 
011 
007 
- 0 
0. 044 099 | i 2 x at 0. 018 
041 OR6 5 4 76 026 
032 -- 049 ; 6 7 . 014 
.019 002 | : 0 
006 5 038 2 ; 3 | —0.013 
O55 . 205 017 
0 





TABLE 4. Stresses in plate for case 3—Continued TABLE 5. Stresses in plate for case 4—Continued 


704 Results Corrected No. of 704 Results Correctec 


o:/q 


0. 012 
—. 012 
009 
003 
013 
. 025 
. 033 
031 


0 

0 

0. 003 
009 
012 
O11 
003 
O19 


Stresses in plate for case 


nloaded length 7 
Width 


704 Results 
oy/q 


1. 004 
1. 004 
1. 003 
1. 002 
1. 001 
1. 000 
1. 000 
1.0 


1. O11 
1. 010 
1. 007 
1. 004 
1. 003 
1. OO1 
1. 000 
1.0 


1. O18 
1. O17 
1.013 
1. OO8 
1. 005 
1. 003 
1. 001 
1.0 


0. 999 
1. 001 
1. 004 
1. 005 
1. 005 
1. 003 
1. OO1 
1.0 


0. 954 
960 
975 
992 
QOS 

1. 000 

1. 001 

1.0 


0. 857 
869 
902 
950 
974 
YYl 
gay 

1.0 


0. 695 


a 


0.5;l>>¢ 


Fy/q 


0. 047 
. 046 
. 043 
. 032 
. 024 
.014 
. 005 


176 
. 165 

133 
. O83 
. 055 
. 029 
. 009 


Corrected 


0 
0. 010 
O17 
019 
017 
014 
008 


O18 
038 
056 
. 060 
055 
037 
001 


006 
019 
042 
O54 
059 
. 047 
. 003 


016 
. 029 
. 031 
. 020 
. 003 
031 
. 008 


| 
| 
I 
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o2/q 


064 

057 

~. 033 
008 
033 

061 

. O89 

f —. 387 
{ 613 
-0. 138 


—. 131 
—. 105 


JE 6. 


Loaded length a 
Width 


No. of 
terms, m 


—1.0 
0 


—0. 311 


, 


0 
0 
0 
0 
0 
0 


0 
0. 045 
OSY 
. 116 
110 
. 079 
025 


0 
0 
0 
0 
0 
0 
0 
0 


=2.0; 1 


f—0. 500 


500 


Stresses in plate for c 


ase 
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TABLE 8. Stresses plate for casé Continued 


3 


No. of 
terms, ™ 


016 
030 
~. 038 
031 
O19 


7. Discussion 


7.1. Plate With Distributed Load 


The results given herein are directly applicable to a 
finite plate loaded by uniformly distributed equal 


loads along portions of two longitudinal edges. The 
ratio of the length of the loaded portion to the total 
length of the plate is given in column 5 of table 1. 
The ratios of the loaded length of plate to the width 
and of the unloaded length of plate to the width are 
given in columns 6 and 7 of table 1. 

A comparison of figures 7 to 9 shows that, in ac- 
cordance with St. Venant’s Principle, the discoa- 
tinuity in the loading affects the stress distribution 
in the plate for a distance to either side of the dis- 
continuity equal to the width of the plate. The 
stresses for case 1, table 2, have the same distribution 
and are of about the same magnitude as those for 
case 5, table 6. It is inferred, therefore, that if the 
ratios of the loaded length of plate to the width of 
the plate and of the unloaded length to the width 





TABLE 9. Stresses in plate for case 8 


Loaded length a 
Width ’ ec 


No of 
terms, ™ 


are greater than 2.0, the stress distributions will be 
similar to and the stresses will be of the some magnitude 
as those shown in figure 8 and listed in tables 2 and 6. 

For those points in the plate at a distance twice 
the width of the plate from the boundary of the load, 
the stresses in the direction of loading, ¢,/q, will be 
either 1.0 or 0, depending upon whether the points 
are in the loaded or unloaded parts of the plate, 
respectively. The stresses transverse to the direc- 
tion of loading, o,/q, and the shear stresses, r,,/q, 
will be zero. 
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The stresses in the direction of loading are, at 


the most, only 4 percent higher than the applied | 


load. However when the loading is over only a small 
portion of the plate length, case 8, the stresses de- 
crease from values equal to the applied load at the 
boundary to values equal to only about 76 percent 
of the applied load at the center, figure 9a. 

The stresses normal to the direction of loading, 
o,/q, are of the same order of magnitude and have the 
same distribution but are of opposite sign on either 
side of the discontinuity for cases 1 and 5 (see figure 
8b). The distribution of the stresses changes some- 
what as the ratio of loaded length of plate to width 
of plate is decreased from 2.0 to 1.0, cases 5 to 7, 
but the maximum stress does not. For case 8 where 
the ratio of loaded portion of the plate to width of the 
plate is 0.5, there is a buildup of stresses in the 
loaded portion of the plate near the discontinuity 
in the loading, figure 9b. This same increase in ¢,/q, 
is found in the loaded portion of the plate as the ratio 
of nonloaded length to width of the plate is decreased 
from 2.0 to 0.5, being most significant when the ratio 
is smallest, case 4, figure 7b. The distribution of 
the axial stresses across the width of the plate for 
case § is shown in figure 10. 

The shearing stresses increase to almost 30 percent 
of the applied load near the boundary of the plate 
in the region of the discontinuity in the loading, 
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Ficure 10. Distribution of axial stresses across width of plate 


for case 8. 





figures 7c to 9c. These stresses should approach 
a theoretical value of + 1/ras the edge of the plate is 
approached and then go to zero on the boundary. 
In other regions of the plate, the shearing stresses 
are insignificant. 


7.2. Comparison With Previously Published Results 


Barton [8] investigated the problem of a long cir- 
cular cylinder loaded with a band of uniform pressure 
over a finite length. His figures 19 and 20 for a load- 
ing length of one diameter and figures 22 and 23 
for a loading length of one-half diameter are of partic- 
ular interest. These load lengths are similar to 
cases 7 and 8 of this paper, with a loading length 
equal to the width of plate and one-half of the width, 
respectively. Barton’s intermediate values of 
stresses differ considerably from the values shown in 
figures 8 and 9, but the general patterns and extreme 
values are strikingly similar. Considering the dif- 
ference in geometries between these problems, this 
general similarity between the two sets of results is 
remarkable. 
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FicuRE 11. Stresses in longitudinal direction along center line 


of welded joint. 
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Fiaure 12. Shearing stresses in welded joint. 


(a) Along center line; (b) 5 in. from center line; (c) 8 in. from center line. 


8. Application of Results to a Welded 
Structure 


Having computed the stresses in a plate loaded 
along portions of two opposite boundaries, an esti- 
mate of the stresses in a simple welded structure can 
be made. The unloaded flange is considered cut 
and attached to the end of the loaded base plate, 
figure 3. For a Poisson’s ratio of 0.29, the stresses 
corresponding to figure 4b are 


Or P 


oy=—q=0.29p. 

The stresses in a plate corresponding to figure 4c 
can be taken from tables 2 through 9 for appropriate 
geometrical parameters, noting that (r—a) is the 
distance measured from the point of discontinuity 
of the load g. The stresses corresponding to those 
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Figure 13. Stresses in transverse direction in welded joint. 


(a) Along center line; (b) 5 in. from center line; (c) 8 in. from center line. 


in figure 4c are then added to the stresses correspond- 


ing to those in figure 4b. 


Using the tables as described above, the stress 
distribution in a welded structure similar to figure 2 
was calculated. These theoretical results are given 
in figures 11 to 13 where they are compared with 
the experimental results described in the appendix. 


| The results are found to be in reasonable agreement. 
| 


9. Summary 


Values of the stresses in plates having uniform 
symmetrical loading along portions of two opposite 


| longitudinal edges are presented in tables 2 to 9. 
Lan) tn) 


The discontinuity of loading affects the stress dis- 
tribution in the plates for a distance to either side 


| of the discontinuity equal to the width of the plate. 


By suitable transformations, these results can be 
used to evaluate the stresses near the joint of two 
perpendicular plates subjected to uniaxial tension. 
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The computed results show that a partial loading 
of the edges of the plate produces a biaxial stress 
condition. In any welded structure subjected to 
the loading considered, the state of stress would be 
biaxial tension. 

Theoretical results obtained by the use of tables 
2 to 9 show reasonable agreement with experimental 
results. 


The authors are indebted to Ruth M. Woolley for 
her competent work on the preliminary phases of 
this problem so that it could be expeditiously pro- 
gramed for the computer, and to L. K. Irwin, A. F. 
Kirstein, and others in the Engineering Mechanics 
Section for their valuable suggestions and help, 
especially in the experimental phase of the problem. 
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11. Appendix 


In order to verify the theoretical results obtained, 
a test of a simple welded structure was conducted 
in a 1,150,000-lb capacity horizontal testing machine. 

















(6) Stations ot which stresses were determined. 


FiGurRE 14. Welded structure tested. 


The specimen was made of 40.8-lb ship-plate, 1?in. 
in thickness. The base plate was 8 ft. 4 in. long 
and 20 in. wide. ‘Two flanges, each 20 in. wide and 
10 in. long, were welded to the base plate (see fig. 
14a). The overall shape of the specimen was similar 
to that shown in figure 2. Pulling heads having the 


| shape of eyebars were welded to the ends of the base 


plate. The whole assembly was completely stress- 
relieved to a stress-free state. The pulling heads 
were then connected to the testing machine by pins. 
Bonded wire-resistance strain gi and three- 
component rosettes of bonded wire-resistance strain 
gages were installed at strategic locations on the 
specimen, figure 14b. * With this setup it was possible 
to determine the Young’s modulus and Poisson’s 
ratio, as well as the magnitude of strains in different 
directions in the specimen. Poisson’s ratio was 
determined to be 0.29. The applied load was 500,000 
lb or p=25,000 Ib/in. in the longitudinal direction. 
The strains in the specimen were read directly to 10 
uin./in. and were estimated to 2 uin./in. Due to the 
high Young’s modulus of steel, these values are 
equivalent to approximately 300 psi and 60 psi, 
respectively, in stresses. Therefore, the results can- 
not be claimed to have a high degree of accuracy. 

The test results were reduced and indicated on the 
graphs of theoretical values (figs. 11 to 13), which 
are computed as described in section 8. It can be 
seen that the theoretical and experimental results 
are in reasonable agreement. 
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Low-Carbon Sieels 
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The mechanism of stress-corrosion cracking in low-carbon steels was investigated using 
notched specimens, stressed in tension to produce a slight amount of plastic deformation, and 


exposed in a boiling 20 percent aqueous solution of NH,NOs. 


Initiation and propagation of 


cracks were studied by removing specimens for metallographic examination after stated 
periods of exposure and by simultaneously recording extensions and electrochemical solution 


potentials. 
occurred in 200 minutes or more. 


Cracks were initiated after about 20 minutes’ exposure and complete failures 
Both crack initiation and propagation were postulated 


to result from an electrochemical process in which the anodic areas were submicroscopic. 
Tensile fractures were initiated in the corroded specimens several times before complete 
failure occurred but were prevented from going to completion because of the energy involved 
in forming the surfaces of the cracks and the strain hardening and the strain aging of the 


steel at the tips of the cracks. 


1. Introduction 


Failures of riveted steel pressure vessels have 
frequently been attributed to stress-corrosion crack- 
ing, called in these instances ‘‘caustic embrittlement.’’ 
The construction of pressure vessels by welding has 
prevented the leaks that were a factor in causing the 
cracking of riveted structures but has by no means 
eliminated stress-corrosion cracking of low-carbon 
steel vessels. The mechanism of ‘the initiation and 
propagation of stress-corrosion cracking in low-carbon 
steels is not well understood. Several theories have 
been advanced to explain the mechanism; none are 
universally accepted. The data reported in this 
paper are used for the reexamination and refinement 
of some of the current theories of stress-corrosion 
cracking. 


A brief review of conditions under which stress- | 


corrosion cracking has been reported in low-carbon 
steels was given in an earlier paper [1].'!_ Corrodents 
have included hydroxides, nitrates, cyanides, com- 
plex organics, and even humid weather conditions.? 

In our earlier laboratory study of the stress- 
corrosion cracking of low-carbon steels [1] extensive 
exposure periods were required to produce failures in 
conventional (but subsize) tensile specimens in a 
boiling aqueous solution of NH,NO,. Failures of 
specimens in less than 300 hr were rare and some 
specimens were removed from the corrodent after 
1,800 hr exposure with no evidence of the develop- 
ment of stress-corrosion cracks. Failures, when they 
occurred, were in regions of stress concentrations 
resulting from the combined action of stress and the 
general attack of the corrodent. <A typical failed 
specimen is shown in figure 1. 

Anodic polarization of specimens produced prefer- 
ential grain boundary attack at a small number of 
boundaries, whether or not the specimens were 


1 Figures in brackets indicate the literature references at the end of this paper. 

? Failure of steel hooks, designed to hold in place the concrete block facing on 
the Bureau of Engraving and Printing Annex, Washington, D.C., was attributed 
to stress-corrosion cracking [2]. 





subjected to externally applied stresses. In the 
stressed specimens the penetration of the corrodent 
along the grain boundaries as a result of the applied 
current (average current density, 0.06 ma/cm?) 
eventually caused such stress concentrations that 
failures occurred in regions having fairly uniform 
cross-sectional areas. Again, however, the exposure 
periods to failure were more than 200 hr. 

These data suggested that specimens designed to 
contain definite and reproducible regions of stress 


concentration, e.g., notches, would most probably 


fail after relatively short exposure periods. Accord- 
ingly, a notched specimen was designed that sub- 
sequently failed after a few hours’ exposure under 
stress in the corrodent. 

In earlier studies time lapse photographic records 
of the extension of a conventional specimen, that 
failed in 315 hr, showed that most of the strain indica- 
tive of the development of stress-corrosion cracks 
occurred during the last 3 hr of the life of the speci- 
men. With aspecimen that would fail in a few hours 
it was feasible to determine the incubation period to 
the initiation of cracking and to determine the exten- 
sion-time and electrochemical solution potential-time 
relationships during the life of the specimen. 

In the early work [1] with fine-grained materials 
the exposure periods to failure were much longer for a 
“ship plate steel’ * having a high transition temper- 
ature from brittle to ductile fracture than for a 
similar steel having a lower transition temperature. 
Although the temperature of the corrodent was well 
above the higher transition temperature, it seemed 
desirable to compare these materials as notched 
specimens where the inherent brittleness of the 
steels, if present, would be expected to affect the 
results. ‘Two steels, one silicon killed, the other 


3 Plates had been removed from welded ships. A crack had run completely 
across one of these plates (high transition temperature from brittle to ductile 
fracture) while the crack had penetrated into and stopped in the second plate (low 
transition temperature from brittle to ductile fracture). Fora report of the NBS 
work on “ship plate,” see ‘Correlation of metallurgical properties and service 
performance of steel plates from fractured ships’ by Morgan L. Williams, 
Welding J. Research Suppl., p. 445, (Oct. 1958). 
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FIGURE 1. 


Conventional stress-corrosion specimen that failed afte r 300 hr in the aqueous bowling NH,NO solution. 


Cracking was limited to reduced section containing regions of stress concentration 


aluminum killed, with about the same carbon content 
as the ship plate steels were also available. The 
susceptibility of these steels to  stress-corrosion 
cracking was compared with the ship plate steels and 
data obtained used in an attempt to correlate 
mechanical properties, chemical composition, and 
susceptibility to stress-corrosion cracking. 


2. Materials and Methods 


2.1. Materials 


Four of the steels used in the previous work [1] 
were also used in the current investigation. Their 
chemical compositions, tensile properties, transition 
temperatures (if determined), hardnesses, and grain 
sizes are given in table 1. A and B were ship plate, 
C was the silicon-killed, and D the aluminum-killed 
steel. 


2.2. Experimental Methods 


The silicon- and aluminum-killed steels (C and D 
of table 1) were supplied as %-in. diam rods. These 
steels were normalized before use by holding them 
for 1 hr at 1,700 °F and air cooling. Notched speci- 
mens, having the dimensions shown in figure 2, were 
machined from these normalized rods. The ship 
plate steels were % and % in. thick, respectively, for 
the high transition and low transition steels. In 
order to avoid differences in notch geometry, coupons 
from these steels were machined to %-in.-diam rods 
and the notches were then machined as indicated for 
the killed steels. Neither ship plate steel was heat 
treated in any manner at NBS. 

In order to investigate the relative susceptibility 
of these various steels to stress-corrosion cracking, 
all specimens were insulated with plastic tape except 
over the notched areas, and placed in Pyrex cells 


TABLE 1 


A. Chemical composition of steels used in investigation 


Manganese Phosphorus Sulfur 


0.01 0. 029 
OO 025 
008 022 

. OO . 029 


Silicon Aluminum Chromium Copper 


0. 03 0. 070 0.05 
061 O16 O15 
17 ol ND 
. 018 04 ND 


. Mechanical properties of steels used in investigation 


Yield 
strength 


Ultimate 
tensile 
strength 


pst pst ¢ 
45, 400 60, 200 63 
32, 800 62, 500 5g 
49, 700 
46, 200 63, 200 68 


NoTtE: ND—Not determined. 


Reduction 
in area 


66, 800 6614 


Trans ‘tion 
temperature 


Grain size 
ASTM No. 


Hardness 
Rockwell B 


scale 
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Ficure 2. Drawing of notched specimen used in the investi- 


gation. 

that were wrapped externally with resistance wire 
and heated electrically. The corrodent was a boiling 
aqueous 20 percent solution of NH,NO,. Loss of 
liquid was prevented by reflux condensers connected 
to the tops of the cells. The cell used in much 
of the work had a side arm that could be connected 
through a salt bridge to a calomel electrode. The 
electrochemical solution potential of the specimen 
(in terms of a calomel electrode of the saturated 
KCl type) was thus readily determined manually or 
recorded on one channel of an oscillograph. 

Previous work in this laboratory [3] and studies 
of the data of other investigators have indicated 
that the order of procedure may be important in 
obtaining data, particularly those including periods 
to failure and electrochemical solution potentials. 
Most of the data to be reported were obtained using 
the following procedure. The corrodent was heated 
to boiling (103.5 °C) and the electrochemical solution 
potential of the specimen was noted. When the 
potential reached a value of —0.40 v, the specimen 
was stressed in tension by dead loading using a 
lever system. A few specimens, however, were 
stressed prior to adding the corrodent and bringing 
it to its boiling point. 

In laboratory investigations, stress-corrosion crack- 
ing of low-carbon steels occurs only in specimens 
subjected to stresses high enough to produce some 
plastic extension. In order to determine the stresses 
to be used with the notched specimens, true stress- 
logarithmic strain curves * were obtained for two 
notched specimens from each lot of steel. The 
stresses selected for the various lots of steel were 
those that would produce approximately the same 
deviation (logarithmic strain of 0.01) in each case, 
from the modulus line in the room temperature true 
stress-logarithmic strain diagram. The logarithmic 
strain of one specimen removed from the corrodent 
after 10 min under stress, at temperature, was 0.033. 

Most of the deformation of the specimen occurred 
on the initial loading. However, notched specimens 
continue to elongate but at a diminishing rate for : 
considerable period after the stress is applied. This 
elongation was followed by reading (or recording, 
using time-lapse photography) dial gages actuated 


~ . load , areao 

4 True stress is defined as , and logarithmic strain as In Areao is 
area area, 

the original cross-sectional area at the root of the notch and area is the area at 


any time ¢. 


| steels studied. 





by the motion of the ends of the loading levers. 
Inasmuch as the lever ratio was 25 to 1 and the dial 
gage could be read to 1 * 107‘ in., a change in 
length of the specimen-specimen support system of 
4 < 10°° in. could be observed. 

In a few instances the elongation of the specimen® 
alone and in many instances the elongation of the 
specimen-specimen support system was recorded 
simultaneously with the electrochemical solution 
potential of the specimen. 


3. Results and Discussion 


3.1. Comparison of the Relative Susceptibility of 
Four Lots of Steel to Stress-Corrosion Cracking 


Four specimens of each of the four lots of steel 
whose compositions are given in table 1 were exposed 
under stress in the boiling NH,NO, solution until 
failure occurred. As was indicated earlier, while the 
applied stresses were different for each lot of material 
they were determined from true stress-logarithmic 


| strain curves to give approximately equal plastic 
| strains. 


The average exposure periods to 
and the standard deviations are given in 


failure 
table 


steel 
strain, in a 


TaBLE 2. Exposure periods to failu-e for the lots of 
exposed, subjected to approximately the same 
boiling NH,NQs solution 


Average 
exposure period 
to failure 


standard 
deviation 


min. 


85 


9 « 


» 


17,6 


Specimens from the ship plate steel (lot B), with 
the highest transition temperature, were the most 
resistant to stress-corrosion cracking of any of the 
The fact that this steel had brittle 
characteristics below 98 °F (37 °C) apparently had 
no effect on its properties in the NH,NO, solution 
at 103.5 °C. 

There were slight differences in the chemistry of 
the four steels. The carbon content of all was 
approximately 0.20 percent and according to Parkins 
[5] they should be resistant to stress-corrosion 
cracking. He reported that the resistance of steels 
to stress-corrosion cracking in a Ca(NO,).+NH,NO, 
solution increased with increased carbon content up 
to 0.18 percent and that steels containing more of 
carbon were relatively resistant to stress-corrosion 
cracking. 

The principal differences in the 
the four steels are in their silicon 
content. The silicon contents of A and D are 
relatively low; B, however, has a silicon content 
between that of A and C. The aluminum ene 
of A and D are higher than that of either B or 
Whether the difference in the resistances of es 


5 In this case a specimen 9 to 10 n. long was used. 


composition of 
and aluminum 
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steels is due to the differences in silicon and aluminum 
contents or other factors cannot be determined from 
these data. 

A fractured surface, characteristic of that found 
on specimens from lots B, C, and D, is shown in 
figure 3. In all of these specimens, the area of 
fracture due to stress-corrosion cracking was clearly 
delineated from that produced in the final tensile 
failure of the specimen. The fractured surface 
from a typical specimen of lot A steel is shown in 
figure 4. With careful study the area of fracture 
due to stress-corrosion cracking could be determined. 
The unusual appearance of the fractured surfaces 
of these specimens was due to shearing of the speci- 
men, in the final tensile fracture, along planes 
approximately parallel to the tension axis. A 
longitudinal section of the specimen, including the 
fractured surface, is shown in figure 5. The steel 





contained an unusual amount of stringers of non- | 


metallic inclusions and the shearing action occurred 
along these stringers. 


FIGURE 3. su 


Fractured face of specimen characteristic 
that found on specimens from lots B, C and D. 
K 4h. 


FIGuRE 4. 
Not 


A steel. 


and D 


Surface of fractured specimen of lot 


difference in appearance from fractures from lots B, ¢ 


0, f 


Figure 5. Longitudinal section through fracture in specimen 

shown in figure 

Void in the center of the figure resulted from 
Appearance of failed surface at the left of the 

corrosion crack probably penetrated to that point 
right of the void occurred during tensile fracture 


; 

4. 

shearing action at failure. 

void indicated that a stress- 
Distortion of grains to the 

Etched 1 percent nital. X< 50. 


3.2. Initiation of Stress-Corrosion Cracks in a Notched 
Specimen 


The use of notched specimens has made possible 
the study of the initiation and propagation of stress- 
corrosion cracks in low-carbon steel exposed in the 
ammonium nitrate solution. 


The room temperature (25 °C) electrochemical 
solution potential of a typical steel from this group 
(~0.20% carbon) was —0.6 to —0.7 v with respect 
to a calomel electrode of the saturated KCl type. 
As the corrodent and specimen were heated to the 
boiling point the electrochemical solution potential 
became more noble. The change was generally 
rather slow until the solution began to boil at about 
103.5 °C. With contirued boiling, the potential 
drifted slowly in the positive direction and could 
reach —0.30 v after several hours exposure in the 
boiling corrodent. The change in potential from 
—0.6 to —0.380 v on heating the corrodent from 25 
to 103.5 °C indicated that a protective film (or 
polarized condition) had formed on the exposed 
surface of the specimen. 

Stressing of a specimen that was at the boiling 
temperature of the corrodent and had an electro- 
chemical solution potential of —0.38 v or less changed 
this potential almost instantly as much as 70 mv in 
the active direction. A typical oscillographic trace 
for this change is shown in figure 6. The stress 
applied to the specimen was calculated to be suffi- 
ciently large to produce plastic deformation at the 
root of the notch. Thus, the change in potential 
was due to the rupturing of the protective film (or 
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FIGuRE 6. 


Markings or wer scale indicate successive minutes. 


disruption of the polarized condition) at the root of 
the notch. 

Some minutes were required for the specimen po- 
tential to change from the maximum value (see fig. 
6) to its value prior to the application of the stress. 
Extension measurements made directly on the speci- 
men have indicated that it continued to extend, but 
at a diminishing rate, after the initial extension 
accompanying the stressing of the specimen.® 


The measured potential of the stressed specimen 
depends on the relative rates of film rupture and re- 


pair. Measurements of the rate of change in poten- 
tial of a specimen, taken immediately following its 
complete failure, indicated that 10 min or more may 
be required for the freshly exposed metal to be polar- 
ized to the prefractured condition of the surface. 
This is not unexpected considering the relatively slow 
rate ’ at which the potential of the specimen changes 
during the heating period. Hence, the data indicate 
that the rate of film repair rather than the continued 
extension of the specimen was the primary factor 
that delayed the return of the specimen potential 
to the prestressed value. The electrochemical solu- 
tion potential continued to change slowly in the 
noble direction, as was indicated above, after the 
prestressed potential had been reached. 

There was no abrupt change in the extension rate 
of the specimen, such as was found for stainless 
steel [6], indicating the initiation of stress-corrosion 
cracking.” Specimens were removed from the cor- 


—_—_—_—_—___.. 


6 The stress was applied by releasing a hydraulic jack supporting the end of 
the lever with the load in place. Hence, the stress was applied rapidly but 
without impact 

? This slow rate of film repair contrasts sharply with the very rapid film repair 
found in the AZ31 magnesium alloy [3]. In this latter case film repair takes place 
in less than 1 sec 

§ Changes in the extension rate for stainless steels were noted using conventional 
tensile specimens with the extension resulting from the formation of many cracks. 
Normally on a section through the root of the notch only two cracks were found. 
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Ose illographic record of change in pote nlial accompanying loading of a specimen. 


Irregularities in potential curve at equilibrium positions are caused by boiling of the solution 


roding medium after various exposure periods, sec- 
tioned longitudinally, and examined metallograph- 
ically for evidence of stress-corrosion cracks. In 
this manner, it was determined that cracks may 
start to develop within 20 min after the specimen 
was stressed. Figure 7 shows a _ stress-corrosion 
crack in a notched specimen removed after 30 min 
exposure under Examinations of newly 
initiated stress-corrosion cracks have usually re- 
vealed corrosion products on the surface of the 
metal at the mouths of the cracks in these specimens. 
Cracks followed intercrystalline paths (see fig. 8) 
and usually developed in pairs in a notched specimen 
as is shown in figures 8, 9, 10, and 11. Hetenyi 
and Liu |7] have shown that the maximum shearing 
stress is not at the root of the notch but is displaced 
about % the root radius from that point. This 
may explain the location of the stress-corrosion 
cracks in figures 8, 9, 10, and 11. Kies, Smith, 
and Irwin [8] have postulated that the direction 
of the maximum tensile stress near the root of the 
notch is not parallel to the tension axis but makes 
an angle of about 60 degrees with it. This may 
explain the direction taken by the cracks. 

The results obtained by stressing a specimen in 
the corrodent prior to heating it were different in 
some respects than those reported above. For 
example, applving the stress to specimens at room 
temperature did not produce any appreciable change 
in the electrochemical solution potential. The 
potential of the stressed specimen, however, be- 
came more positive with increasing temperature as 
was noted above for the unstressed specimen heated 
in the corrodent. Specimens remained in the solu- 
tion at temperature for 60 min after rapid boiling 
had begun and during this period there were no 
changes in potential, in the active direction, indic- 
ative of film rupture. Specimens were then removed 


stress. 





Figure 9. Stress-corrosion cracks that had developed in a D 
specimen exposed under stress for 60 min without any 


9 & evidence of sudden eatension. 


FicurRE 7. Stress-corrosion crack developing in specimen 
exposed under stress fo O min in the boiling 20 percent 
NH,NOz solution. 


Root of the not to the t netched 


FIGURE 10. Stress-corrosion cracks that de veloped in a D 
specimen with sudden extension occurring during the last 
minute of exposure. 


Compare the extent of crack in this specimen with figure 9. Unetched. X 100. 


and sectioned longitudinally to determine the 

extent of stress-corrosion cracking. The depth of 

penetration of the cracks was not appreciably 

different from that found in specimens of the same 
lot of steel (C) that had also been stressed for a 

in a C specimen after 1 hr.exposure in the corroding medium. — of 60 min but after the solution had come to 
Unetched x 250 a DOL. 
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FIGURE S. Inte rerustaliine stress-corrosion cracks developing 





Figure 11. Cracks that developed in specimen from lot D steel 
after 2 hr exposure in the corrodent. 


Unetched x 90. 


Most investigators in the field agree that stress- 
corrosion cracking is initiated by an electrochemical 
process. The fact that no cracks developed in a 
specimen to which cathodic protection was applied 
during the entire exposure period under stress sub- 
stantiates this idea for low carbon steels. 


If a continuous anodic region exists at the grain 
boundaries, and if this region is attacked in the 
corrodent in the absence of stress, there is no prob- 
lem in explaining the initiation of intererystalline 


stress-corrosion cracking. Material at grain bound- 
aries may be in a higher energy state (and hence 
more active) than that on the exposed faces of 
crystals. There was no evidence in this investi- 
gation, however, that the grain boundaries were 
sufficiently energetic, in the absence of stress, to 
provide a continuous anodic path for either the 
initiation or propagation of stress-corrosion cracks. 
Hence, the appearance of active anodic areas where 
certain grain boundaries contact the surface of the 
metal must be associated with the application of 
the stress. Two explanations of this mechanism 
have been proposed. 

The first of these postulates that the polarized 
condition on the surface film developed during the 
heating of the corrodent protects the specimen from 
corrosion generally as long as it is intact. Applying 
stress to the specimen immersed in the boiling 
corrodent destroys the polarized condition (or 
ruptures the film) in the highly stressed region at 
the root of the notch. Electrochemical solution 
potential measurements obtained at the final failure 
of specimens indicated that this exposed surface 
may actually be 0.20 to 0.25 anodic to the material 
on the sides of the notch. Thus, an ideal condition 
exists for corrosion by an electrochemical process | 
with all of the required elements present: a rela- | 


tively small anode, a large cathode, and an elec- 
trolyte. 

The data, however, indicate that cracks do not 
begin to develop until the film rupture associated 
with the initial loading has been repaired. While 
several minutes were usually required for the repair 
of the protective film (or polarized condition) at the 
root of the notch metallographic data given above 
indicate that approximately 20 min of exposure was 
required for the initiation of cracking. Furthermore, 
stress-corrosion cracking was approximately as 
extensive in those specimens stressed in a cold 
solution, which never showed an appreciable change 
in potential that could be attributed to film rupture, 
as in specimens showing a change in potential of 
0.07 v on initial stressing. If rupturing of the 
protective film on the surface of the specimen is 
responsible for the initiation of stress-corrosion 
cracking the active areas must be small compared 
to those exposed on the initial loading of the specimen, 
in fact, probably little wider than grain or crystal 
boundaries.’ These active areas must also result 
from the extension of the specimen after the initial 
loading. Transducer data obtained in one experi- 
ment indicated that 20 min after the initial extension 
the specimen was continuing to extend at the rate 
of 1 to 2 10-° em/see. Theory, metallographic 
evidence, and the data listed above, all indicate that 
this extension occurred in a region of the steel less 
than one root radius in length, i.e., less than 10? 
em and hence with the strain rate greater than 
10-‘/see. If the extension was limited to the width 
of the stress-corrosion cracks, the strain rate could 
by as large as 10~?/sec. A conservative estimate 
would be 107%/see and the strain would be expected 
to be confined largely to those regions at the tips 
of the cracks. The strain rate calculated above is 
postulated to be high enough to produce an anodic 
region at or on a submicroscopic area adjacent to a 
grain boundary by rupturing the protective film [9], 
destroying the polarized condition [10], or inducing 
a condition favorable for anodic dissolution {11]. 
Such anodic areas would have so small an effect 
on the average recorded electrochemical solution 
potential of the specimen as to be masked by the 
background (see fig. 6) but could be highly effective 
in the initiation of stress-corrosion cracks by an 
electrochemical process. 

A second possibility that must be considered for 
the initiation of stress-corrosion cracking is a strain- 
induced precipitate at the grain boundaries. Waber 
and McDonald [12] suggested that an Fe.N 
precipitate is the responsible anodic material at the 
grain boundaries. More recently, in explaining 
the mechanism of propagation of stress-corrosion 
cracks in stainless steel. Uhlig and White [13] 
had postulated the formation of a precipitate 
cathodic to the matrix, with the immediately 
adjacent material anodic to both the precipitate 
and the matrix, ie., an “impoverished region.”’ 
Such a theory might also apply to low-carbon steel. 

* Initial loading ruptures an area at the root of the notch that is relatively 
large, probably containing both crystal boundaries and crystal interiors. Inter- 


crystalline corrosion can take place, however. only if the exposed crystal faces 
are cathodic to the material at or adjacent to the crystal boundaries. 


353 





3.3. Propagation of Stress-Corrosion Cracks 


Two general mechanisms 
for the propagation of stress-corrosion cracks. 
These are designated the electrochemical and the 
electrochemical-mechanical theories. The propaga- 
tion of cracking in the electro- 
chemical postulated as resulting from 
an electrochemical process in which the advancing 
tip of the crack is the anode and the sides of the 
crack the cathode. This provides a small anode- 
large cathode condition, ideal for rapid corrosion, 
but nevertheless, may require current densities of 
the order of several amperes/cm? [8, 10] at the tip 
of the crack. 

In the electrochemical-mechanical theory, cracks 
are postulated to both start and penetrate into the 
metal by an electrochemical process until sufficiently 
high-stress concentrations have been built up for 
them to propagate as a mechanical fractures. 
Keating [14] suggested that such cracks would 
travel until they were stopped by an inclusion, 
unfavorably oriented grain boundary, or other 
stacle. More recently the cracks started by elec- 
trochemical processes have been considered “Griffith 
eracks,’’ and further cracking has been treated as 
“brittle fracture.” [15] According to this theory 
the decrease in elastic strain energy of the specimen 
(due to the reduction in cross section of the support- 
ing metal) must be greater than the energy associ- 
ated with forming the surfaces of the new cracks 
plus the energy used in the plastic deformation of 
the metal at the apex of the crack, that is 


have been proposed 


stress-corrosion 
theory is 


ob- 


YP 


where 
elastic, 


y is the energy and the subscripts indicate 

surface, and plastic energies respectively. 
Perryman [i6] has suggested that “‘it is 

that the spread of cracks along boundaries is assisted 


possible 


by the reduction in surface energy associated with 
the wetting of the crack ” This idea has 
been developed further by ( ‘oleman, Weinstein, and 
Rostoker [17]. They state, however, that at 
stage an electrochemical reaction enters 
process. 


some 
into the 


In the present investigation stress-corrosion cracks 
have penetrated as deeply as 1 to 1.5107? em in 
1 hr under stress in the boiling corrosive medium. 
If cracking occurred only after the first 20 min of 
exposure the subsequent penetration rate was ap- 
proximately 5x10 em/sec. If the crack tip is 
assumed to be continuously anodic to the crack sides 
(as a result of one of the mechanisms discussed 
earlier) and if it is further assumed that the resulting 
electrochemical dissolution 100 percent 
efficient, a current density of about 0.15 amp/em? 
would be sufficient to produce crack penetration at 
the rate indicated above. Such a current density is 
considered entirely reasonable. 

Engell and Baumel [18] have determined the 
electrochemical solution potentials of grain bound- 
aries and grain interiors of ‘coarse-grained samples 
of mild steel.’’ From their curves one can estimate 
that, with the steel in the unstressed condition, the 


process Is 


grain boundaries were approximately 75 mv anodic 
to the grain interiors in a boiling 55 percent Ca(NQs). 
solution. These authors state, without explanation, 
that the application of stress appreciably increased 
the difference in potential between the grains and 
grain boundaries. 

In the present work if the exposure periods were 
continued beyond 60 min there were sudden 
marked increases in extension '’ of the specimens 
that were accompanied by changes in the electro- 
chemical solution potential of as much as 70 mv in 
the active direction. Typical oscillographic curves 
showing the changes in extension and potential are 
shown in figure 12. Such changes have: been re- 
corded as many as five times before the specimen 
failed completely. Engell and Baumel [18] have 
reported similar behavior in specimens of a very 
low-carbon steel. Figures 9 and 10 show the develop- 
ment of stress-corrosion cracks at the roots of notches 
of two specimens of lot D steel exposed for 60 min 
in the boiling corrodent. The specimen shown in 
figure 9 extended in a normal manner. During the 
last minute of exposure of the specimen shown in 
figure 10, a sudden extension and change in potential 
occurred, as described above. The photomicro- 
graphs and the a data indicate that 
there was a sudden the metal at the tip 
of the stress-corrosion crack. This produced a sharp 
increase in elongation and exposed unpolarized 
film-free metal to the corrodent with a resulting 
change in potential in the negative direction. <A 
considerable amount of energy is required to form 
the new surfaces. Furthermore, the energy required 
to produce the plastic deformation most probably 
increases rapidly with the crack propagation because 
of the resulting strain hardening and strain aging of 
the material. Hence, the decrease in elastic strain 
energy as the result of fracture . very quickly 
balanced by the sum of energies associated with the 
formation of the new surfaces and the energy of 
plastic deformation at the tip of the crack and rapid 
cracking stops. 

Tensile tests on notched specimens of lot C steel 
made in boiling water indicated that this steel was 
subject to strain aging at 100 °C. Geil [19] reported 
that a 0.02 percent carbon ingot iron, similar in 
composition to that used by Engell and Baumel 
was also subject to strain aging in boiling water. 

The final failure of any material resulting from 
stress-corrosion cracking is by tensile fracture of 
material whose supporting section has been reduced 
by the development and propagation of 
corrosion cracks, see figure 3. The sudden extension 
of the stress-corrosion cracks, just discussed, may be 
considered as the initiation of the final tensile 
fractures of the specimens. ‘These fractures are 
prevented from going to completion by the strain 
hardening and the strain aging of the metal ahead of 
the tips of the cracks. Rapid crack propagation 
then ceases until cracks developing further by the 
relatively slow electrochemical process again produce 
stress concentrations and a condition of tensile 
overload. 


earing of 


stress- 


One D specimen extended abruptly during the 60th min, see figure 10. 
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L——_—_ 5 min 


Oscillographic charts showing increases in elongation of specimen 


FIGURE 12. 





(upper curve) and accompanying change in potential 


(lower curve) 


3.4. Cathodic Protection of Low-Carbon Steels 


The application of cathodic protection can prevent 
stress-corrosion cracks from Starting in low-carbon 
steels exposed under conditions described in the 


earlier sections of this paper. Cathodic protection 
can generally stop further propagation of cracks 
even if they have progressed a considerable distance 
into the steel, and one or two rapid crack extensions, 
such as those described in the preceding section, 
have occurred. 

Cathodic protection is effective in that it polarizes 
the whole specimen (including the tips of the cracks) 
to the potential of the film-free metal. Hence, 
there are neither anodes nor cathodes but a uniform 
potential on the surface of the specimen. 

Cathodic protection stops the electrochemical re- 
action that accelerates the rate of extension of the 
specimen. In this investigation, however, it did 
not stop the extension at once. This is to be ex- 
pected. The tips of stress-corrosion cracks are 
extremely sharp notches that produce very high 
stress concentrations in the metal. Hence, the spec- 
imen will extend until the metal at the tips of the 
cracks is sufficiently strengthened by straia aging or 
cold work to resist further plastic deformation. 
This will have the effect of shifting the stress pattern 
and reducing the stress concentration. Thus, the 
net effect of stopping the electrochemical reaction 
would be a continued extension of the specimen but 
at a steadily diminishing rate. This has been 


found to occur after the application of cathodic 
protection, not only in low carbon steel but also in 
stainless steel and in the AZ31 magnesium alloy. 

Some specimens containing very severe stress- 
corrosion cracks, before cathodic protection was ap- 
plied, would probably fail eventually in spite of 
cathodic protection as has happened with one stain- 
less steel specimen. The important point is that the 
final failure, that could be expected in seconds or 
minutes, may be delayed for minutes, hours, or 
possibly days, by cathodic protection. 

The fact that the rate of extension is slowed by 
cathodic protection is a very strong argument for 
the importance of electrochemical processes in the 
propagation of  stress-corrosion cracking. Once 
stress-corrosion cracks have been initiated, sharp 
notches, effective as regions of stress concentration, 
are present whether or not cathodic protection is ap- 
plied. The essential difference is that with cathodic 
protection the electrochemical solution potential of 
the specimen is the same over the entire surface and 
if further crack development occurs it must be en- 
tirely mechanical and generally very slow. In the 
absence of cathodic protection there are differences 
in potential at different areas on the surface of the 
specimen and crack propagation is rapid. The 
logical assumption is that electrochemical attack, 
occurring at anodic regions at the tips of the cracks, 
accelerates the cracking process by the continued de- 


| velopment of regions of high stress concentration. 
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4. Summary 


1. The mechanism of stress-corrosion cracking 
and the relative susceptibility to cracking of four 
low-carbon steels (carbon content 0.20%) have been 


studied using notched specimens and boiling 20 | 


percent NH,NO, as the corrodent. 
2. Stress-corrosion cracking failures of these steels, 


stressed to produce a small amount of plastic ex- | 


tension at the roots of the notches, occurred in 4 to 


8 hr compared to several hundred hours for un- | 


notched specimens. 

3. At a temperature well above 
temperature (for either steel) a steel with a high 
transition temperature from ductile to brittle frac- 
ture, in the Charpy V notch test, was more resistant 
to stress-corrosion cracking than a steel with a much 
lower transition temperature. There was no corre- 
lation between the resistance of the four steels to 
stress-corrosion cracking and their chemistry. 

4. Cracks at the root of the notch were most prob- 
ably initiated after major damage to polarization or 
the protective film (that resulted from initial stress- 
ing) had been repaired. Cracks generally develop 
in pairs on opposite sides of the root of the notch 
and in the areas where the shear stresses are postu- 
lated to be the greatest. 

5. Crack initiation is at grain boundaries and is 
postulated to be by an electrochemical reaction with 
submicroscopic anodes formed as the result of plastic 
extension at the strain rate of the order of magni- 
tude of 0.001/sec. 

6. In its early stages crack propagation is inter- 
crystalline. It is postulated to be an electrochemical 
process requiring a current density of about 0.15 
amp/cm? at the tip of the crack. The anode at the 
tip of the crack may result from plastic extension at 
a strain rate of the order of magnitude of 0.001/see or 
from a stress-induced cathodic precipitate at the 
grain boundary. 

7. Tensile fracture, resulting from overload, was 
initiated as many as five times before it went to 
completion. Earlier failures were prevented by the 
strain hardening and strain aging of the metal ahead 
of the tips of the stress-corrosion cracks. 

8. Cathodic protection can prevent the initiation 
of stress-corrosion cracking in these steels and can 
stop or delay propagation of cracks in material 
already severely damaged. This fact is considered 
to be a very strong argument for an electrochemical 
mechanism for the initiation and propagation of 
stress-corrosion cracks. 


the transition | 


Much of the work reported above was carried 
forward under the joint sponsorship of the National 
Bureau of Standards and the Corrosion Research 
Council. 

The “ship plate’ steels, lots A and B, were supplied 
by Morgan L. Williams of this Bureau. Lots C and 
D of the steel were supplied through the courtesy of 
Samuel Epstein of the Bethlehem Steel Co. 

Joan P. Calvert, Michael J. McBee, John M. 
Wehrung, H. Thomas Yolken, James G. Early, Jr., 
Mrs. June Chung Fugh, and Joseph E. Himes 
assisted in the experimental work in this investi- 
gation. 
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A Dual Centrifuge for Generating Low-Frequency 
Sinusoidal Accelerations 
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This paper describes an experimental 


constrained to have nonrotational motion, 


“dual centrifuge’’ 
being calibrated is carried around a circular path in a horizontal! plane. 


accelerometer 
instrument is 


which an 
If the 


in 


such as is provided by a parallel link device, 
sinusoidal excitation along its sensitive (horizontal) axis is obtained. 


The excitation obtain- 


able is equivalent to linear excitation but at unusually low frequency and large displacement. 


For example, : 


a machine has been built which has a frequency range from 0.5 to 30 cycles per 


second, a displacement (zero to peak) up to 12 inches, and develops an acceleration amplitude, 
useful for calibration, up to 100 gravity at 10 cycles per second and above. 


1. Introduction 


Because of the wide use of low-frequency acceler- 
ometers there is a need for a calibration system 
tailored to the requirements of these instruments. 
This necessity has led to the consideration of several 
methods of generating the required low-frequency 
accelerations. At low frequencies most systems 
suffer from one or more limitations, such as poor 
waveform, restricted amplitude, or excessive super- 


imposed vibration. 

The several systems for obtaining low-frequency 
sinusoidal calibration have been reviewed, which may 
(1) electrodynamic, 


be classified as 
(3) mechanical iain 

The electrodynamic shaker is undoubtedly the 
most useful calibrator for laboratory use. However, 
it usually suffers from insufficient amplitude at low 
frequencies. A typical example of a low-frequency 
electrodynamic shaker has a frequency range of 
0.35 to 500 c/s but provides a maximum displacement 
of only 2 in. (double amplitude), whereas about 20 
n. displacement is required to develop 1 g at 1 ¢/s. 

Transient shakers usually take the form of a 
cantilever spring upon whose free end the instrument 
under test is fixed. They suffer from the fact that 
the excitation is transient, so that it is difficult to 
correlate excitation amplitude with response. In 
addition each shaker is usually restricted to a single 
frequency and several shakers are required to cover 
any appreciable range. 

The mechanical shakers reviewed 
divided into three classifications: 
ages, (2) inertial shakers, and (3) rotary calibrators. 
Four-bar linkages, , scotch yokes and _ slider- 
crank mechanisms, cannot generate a truly sinusoidal 
test motion because of the varying effective moment 
of inertia. In order to get the requisite stiffness in 
the reciprocating parts it is required that a relatively 
large mass be subjected to the test motion so that 
the sinusoidally varying inertial forces are neces- 
sarily large. If the test motion is to be nearly 
sinusoidal it is required that the angular velocities 
of the rotating parts be nearly uniform and it will 
therefore be required that the rotating members have 


(2) transient, ¢ 


may be sub- 
(1) four-bar link- 


e. g. 


| zontal plane. 
| on 





a high moment of and a massive 
device results. 

An inertial shaker suitable for calibration purposes 
has been designed by the Engineering Mechanics 

Section of the Bureau [1] ? whie h consists of a spring- 
mass driven by an eccentric weight. Its frequency 
range is 20 to 110 e/s. Inertial shakers of lower 
frequency are rarely used for calibration because of 
poor waveform. 

A rotary dynamic calibrator which uses the earth’s 
field for excitation has been devised by Wildhack 
and Smith [2]. This device is limited in acceleration 
amplitude to a maximum of | g at all frequencies. 
It is presently in use by the Mechanical Instruments 
Section of the Bureau over the frequency range 1 to 
30 ¢/s. 

The dual centrifuge, which uses its own centrif- 
ugal field for excitation, was described by Woolard 
in 1939 [3]. It consists essentially of a small turn- 
table mounted on a large one, each turning in a hori- 
The instrument under test is mounted 
the small table. When the turntables have 
constant angular velocities the seismic mass of the 
accelerometer responds to a sinusoidally varying 
component of the centrifugal force field generated by 
the rotation of the large table. Woolard published 
an equation showing the acceleration amplitude in 
terms of angular velocities, original displacement, and 
relative displacement of the seismic mass, but did 
not publish a complete mathematical analysis of 
the response. However, he did point out that when 
the small table turns backwards so that it has zero 
absolute rotation, the accelerometer is submitted to 
sinusoidal acceleration which may be of large ampli- 
tude. This is the practically important feature of 
the dual centrifuge and one that has been generally 


inertia, bulky, 


| overlooked. 


The following discussion, which describes a 
mechanical configuration similar to Woolard’s design, 
and develops the equation of motion of the seismic 
mass, shows the dual centrifuge to be a useful tool 
for dynamic calibration, particularly when a large 
amplitude is required at low frequencies. 


1 Figures in brackets indicate the literature references at the end of this paper. 
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2. Analysis of the Response of a Linear’ 


1-Degree of Freedom Seismic Acceler- 


ometer to General Plane Motion Exci- 
tation 


Although the conclusion of this paper is that the 
useful form of the dual centrifuge is that one in 
which the small table has zero absolute rotation, a 
general analysis of the motion without this limitation 
is instructive. 

In the general case, the housing of the instrument 
under test will have motion of 3 degrees of freedom, 
1.e., two of translation and one of rotation. The 
physical situation is shown in figure 1. The large 
turntable turns with an absolute angular velocity 
w) about its center 0. The small table rotates about 
point p with constant absolute angular velocity w,, 
and an angular velocity Q=w,—wp relative to the 
large table. The point 0 is chosen as the origin of a 
plane nonrotating coordinate system, with positive 
directions chosen according to usual convention. Let 
the radius from 0 to p be Ro. In general, the test 
position p’ of the seismic mass m will not coincide 
with p. If a-b represents the direction of the sensi- 
tive axis, then p’ is located an amount d from p, 
perpendicular to a-b, and an amount e parallel to 
a-b. Response of the instrument consists of dis- 
placement S along a-b and the equation of motion is 
therefore in terms of S. 

The following assumptions apply to the analysis: 

(1) Plane horizontal motion, 

(2) Point 0 is fixed, 

(3) Ro, wo, and w, are constant,- 

(4) The instrument has no transverse response, 

(5) The instrument is a linear 1-degree of freedom 
system, 

(6) Accelerometer damping is viscous. 

Referring to figure 1, it will be seen that if the 
seismic Mass m is constrained to motion along the 
axis a-b the motion of the mass may be described 
by the following coordinate equations in which s, 


AY 


FicurRE 1. Schematic diagram of the dual centrifuge. 


| centrifuge in which w,=0, 
takes a 


, and § are relative to the instrument case: 


Ry COS wot —d sin w,t+-(e+-s)cos wot (1) 


y=Ro sin wot+d cos w,t+-(e+s)sin w,ft- (2 


Double differentiation shows that the coordinate 
equations which describe the acceleration of m are: 


> 4 2 > 2 " 
Ro w COS wot +d w) sin wyt—(e+s)w? cos wot 


28, SIN wt 


$ COS wf (3) 


Rows sin wet —dw), cos w,t—(e+s)w} sin wet 
+28w, COS w,f 4 


8 sin W yt - (4) 
The acceleration along the axis a-b to which the 

mass m is subjected is the sum of the projections 

of these coordinate components along a-b, iLe., 

acceleration=Z cos w,f+ i sin w,f or acceleration 

s—(e+s) wy Row cos Qt. 

So that the equation of motion may be derived: 


ms +es+ks=mRous cos Ot+m(e+s)w%. 


and p=c/e critical when ec critical 


k 
Ifo, \ m 


2ymk, this equation can be reduced to: 


§ +2 pa,8 + (w,— ws) s= Rows cos Ot+- ew}: (6) 


The solution of eq (6) is: 


Ag NPV E 4) ap 


wt (p-y P—-14+25) 


r¢ 


> 2 2 2 
Rowol (w,— Q 


+ 


ws) cos Qt- 2pw, Q sin Qt] 
(w,— wp? — 2*)?+-4p*wi 0? 


The steady state portion of eq (7) reduces to: 


Rw, cos (Qt—o) 


> 
3= ; “ s+ = = 


[ (wa—w3 — 27)? + 4p? w20?)|"/? © wi—w? 


in which 


2pw»y Q 
tan @ : 
Wy wy 


It should be noted that in the form of the dual 
l.e., Q=—wo, eq (7) 
identical with the familiar 


form almost 


equation for a spring-mass accelerometer subject to 
sinusoidal linear motion 


> Be 


wnt (P— yp 1) 


> 217.2 H\ cae 
Rowe | (w2— 2?) cos Qt 


+ 


+-2nw,2 sin Ot 
Pon | (9) 


(w2— 0°)?-+-4 pw2Q? 


In which Q 
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The steady state portion may be written: 


Rw cos (wot —d) (10) 


V (w? — ap)? +4 p?w2a5 

From the foregoing it is evident that for the special 
case in which the instrument has only translational 
motion, the motion of m is simple sinusoidal and 
hence may be directly compared with data using an 
electrodynamic shaker. In fact, for a given instru- 
ment, data from the double centrifuge and the 
electrodynamic shaker may be used together to 
construct a single response curve (see figs. 7, 8 
and 9). 


’ 


3. Limitations of the Dual Centrifuge 


Consideration of the response of a linear, single 
degree of freedom system to the general case of the 
dual centrifuge in which the angular velocity of the 
little table is independent of that of the large table 
indicates that this form of the equipment is not 
attractive for accurate calibration. Examination of 


the response equation (7) shows that the magnitude 
of the transient at any instant depends upon w, as 
well as upon the instrument natural frequency w, 


When por. P 
Wy 


and the instrument damping p. 


: ; : ., op 
the transient is oscillatory. At p?+ P L... the 
W,, 


2 
, , a >, wp 
transient is critically damped and when pce 1, 
This 


large values of the 


n 


it is over-damped. that at 


latter 


. TOES = 
suggests 


Wp 


very ratio 


the decay of the 
Wp 

transient may be so slow as to prevent the practical 
attainment of the steady state condition. Moreover, 
examination of the steady bias portion of eq (8) 
shows that, even after steady state conditions how 
been attained, there is a range of W, OVer which the 
seismic mass will be frozen against the instrument 
stops, the frequency range over which this condition 
will exist depending on the magnitude of the offset e 
in figure 1.- 

Without further evaluation of either of the above 
conditions, a decisive consideration is found by 
examination of the periodic portion of eq (8) which 
discloses that the condition of resonance wil be 
realized whenever w,+Q=o,, so that a frequency- 
response curve of the instrument can be drawn only 
if the condition ,=constant, 1.e., w)+2=constant, 
can be maintained. Since the frequency-response 
curve is the objective of sinusoidal input testing, 
this limitation appears to make use of the general 
vase of the dual centrifuge unprofitable. Accord- 
ingly the test machine which will be described in 
the next section was made to operate with w, 
that eq (10) is applicable. 

It is to be noted that, since a transverse axis 
undergoes the same motion as the sensitive axis, 
the instrument under test is subjected to a sinosoidal 


0 so 


transverse excitation of the same amplitude and 
frequency as the excitation along the sensitive axis. 
If the transverse response of the instrument is linear 
and independent of the principal (axial) response, 
the transverse response will be sinusoidal and will 
lag, or lead, the axial response by 90°. Under this 
circumstance the error in calibration can easily be 
estimated by assuming the electrical output of the 
instrument to be proportional to the vector sum of 
the axial and transverse responses. At a transverse 
response of 5 percent of the axial response, the error 
is thus found to be about 0.12 percent. In some 
instruments, however, the response along one trans- 
verse axis is dependent on the axial response, in- 
creasing with the axial response. For this case the 
error is not easily computed as the phase between 
the axial and transverse responses will depend upon 
this rate of increase of transverse response. More- 
over, in this case, the transverse response is most 
likely to be nonsinusoidal. However, it appears 
that instruments with an appreciable dependence of 
transverse on axial response along more than one 
transverse axis are rare. 


4. Description of the Dual Centrifuge 


In view of the obvious advantage in having avail- 
able for instrument test purposes the large accelera- 
tion amplitude at low frequency that can be gener- 
ated with the dual centrifuge and in view of the 
great simplicity of the dual centrifuge when it oper- 
ates with w,=0, a dual centrifuge of this type has 
been built. 

A pilot model was first designed, constructed, and 
tested in order to investigate the potential perform- 
ance of this type of apparatus. This was subse- 
quently replaced with an improved version suitable 
for continuous laboratory use. The pilot 
dual centrifuge shown in figures 2 and 3 consists 
essentially of a small turntable (2) mounted on a 
large turnts able (1) by means of a roller bearing cart- 
ridge (6) and a radius arm (4). The instrument (3) 
is mounted on the small turntable (2). Its leads are 
brought to the overhead arm (15), no slip rings 
being required. Since the pulleys which drive the 
small table are all the same size, the absolute angular 
velocity of the small turntable is controlled only by 
the coaxial control pulley (14). In this case the 
control pulley was fixed and hence the small turn- 
table has zero absolute angular velocity. Motion 
is transmitted through a system of open timing 
belts and pulleys (7, 8). The main table is turned 
by a motor (11) which is a shunt wound d-c motor 
with armature rheostat speed control. Power is 
transmitted via belt drive (12). The roller bearing 
cartridge of the countershaft (5) serves to center 
radius arm (4). The radius arm assembly is made 
symmetrical by the addition of a small balance 
weight (9). Thus the effective radius from the cen- 
ter of the large turntable to the center of the small 
turntable may be varied without altering the dy- 
namic balance of the machine. Counterweight (10) 
serves to balance the combined effect of the radius 
arm assembly and the countershaft assembly. In 


model 
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FicgurE 2. Experimental dual centrifuge. 
preparation for a calibration, the accelerometer is 
balanced by adding an equal weight to the small 
counterbalance (9) and twice that to the large 
counterbalance (10). The whole assembly is sup- 
ported by the machine frame (13). 

The machine shown in figure 4 conforms with the 
basic design illustrated in figures 2 and 3, but has an 
extended amplitude and frequency range. It has a 
frequency range from 0.5 to 30 c/s, a displacement 
amplitude range up to 12 in., and an acceleration 
amplitude range up to 100 g at 10 ¢/s and higher 
suitable for calibration purposes. The zero to peak 
amplitude of the displacement is measured to within 
+0.02 in. by means of two measurements of the dis- 
tance from a fiducial line on the frame to the edge 
of the small table. The angular velocity w)=- 
measured by counting on an electronic events per 
unit time meter the number of pulses from a small 
generator driven by the main shaft which develops 
600 pulses foreach revolution. If the count is over a 
period of 10 sec the accuracy with which the average 
velocity can be determined is 1 part in 6,000 at 1 
rps and with better accuracies at higher velocities. 
Errors in velocity measurement due to a progressive 
change in velocity can be minimized by observing 
the velocity over a period of time both before and 
after as well as during the calibration. Errors in 
calibration due to wow and to flutter at any fre- 
quency other than the frequency of test can be 


Q is 





Sri = 


FicurE 3. Experimental dual centrifuge. 


View of lower control pulleys 


Figure 4. Dual centrifuge. 


detected as a deviation from sinusoidal wave shape 
of the response of the instrument under test. It is 
estimated that in practice the total error in velocity 
measurement does not exceed +0.02 percent of the 
velocity at any velocity above 1 rps. 

An additional source of error is found in the fact 
that even the best rotating devices are not in perfect 
dynamic balance. The machine pictured in figure 4 
has a critical speed around 22 rps, varying with the 
load, near which no calibrations have been made. 
There is some unwanted sinusoidal vibration due to 
dynamic unbalance at any frequency of test, but an 
independent measurement of it can be made with a 
piezo-electric accelerometer, which has no zero 
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Fiaure 5. Dual centrifuge calibration accuracy of sinusoidal 


linear accelerations. 


frequency response, mounted to the big table with 
its axis along a radius. The response of this instru- 
ment is then a measure of the unwanted vibration. 
Its phase relation to the response of the instrument 
under test is unknown, so that it can only be said 
that the sum of the error due to this vibration plus 
any other periodic errors due to either vibration, 
flutter or wow, or similar source does not exceed the 
response of the monitoring piezo-electric instrument 
to within the accuracy of the monitoring instrument. 
In the case of the machine pictured in figure 4 no 
such error greater than +0.02 g has been observed 
during any test. 

The estimated calibration accuracy attained with 
this device over the acceleration range +1/4 to 
+100 g and the frequency range 0.5 to 30 ¢/s is 
shown graphically on figure 5. 


The phase of the response can be measured to | 
| 


within an estimated + 5° by recording the output of a 
photocell excited by a fixed light and a mirror at- 
tached to the large table on a dual trace oscillograph 
with the response of the instrument under test. 


5. Experimental Results 


Figures 6, 7, 8, and 9 show the results of calibra- 
tions on the dual centrifuge. They appear to be 
typical frequency-response curves. Figure 6 shows 


Figure 6. Output of potentiometer-ltype accelerometer under 
dual centrifuge excitation at 0.86 c/s and 0.5 g (zero peak). 
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the response of a potentiometer type acceleration 
transducer when tested on the experimental equip- 
ment pictured in figures 2 and 3, as recorded by a pen 
writing oscillograph. Figure 7 shows the amplitude 
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response of the same instrument obtained by cali- 
brations on the experimental equipment compared 
with the results of calibrations on an electromagnetic 
shaker. Figure 8 shows the frequency response, both 
amplitude and phase, of a bonded strain gage ac- 
celerometer determined by calibrations on the equip- 
ment pictured in figure 4 and by calibrations on an 
electromagnetic shaker. 

Figure 9 shows the amplitude response of an un- 
bonded strain gage accelerometer from calibrations 
on the equipment pictured on figure 2 and on an 
electromagnetic shaker. 
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A rotational manometer may be 


defined 


as an instrument which measures 


7, 1962) 


a pressure 


difference by balancing against it a known pressure difference generated by a rotating element. 


A single fluid is used, 
one type, 
ular viscosity and is predictable. 
and tested in air, 
repeatable to within about 
about 


unsatisfactory. 


and balance is determined by the 
the centrifugal manometer, the rotor-generated pressure is independent of molec- 
This instrument, 
is here adapted for use in water. 
2 percent, for pressure differences above 1 dyne ecm~?, 
10 percent for pressure differences down to 1/10 dyne em~?. 
which the rotor-generated pressures depend upon viscosity 


absence of flow in a detector. In 


which has been previously investigated 
Under the conditions of test, it proved 
and to 
An alternate type, in 
, was found, when tested, to be 


A simple general analysis, supported by observations, shows the dependence 


of manometer sensitivity upon the resistance of the system and the design of the flow detector. 


1. Introduction 


Kemp ' has described a micromanometer in which 
a known pressure difference is generated by centrifu- 
gal accelerations in a rotor to balance an unknown 
pressure difference. Balance is determined by the 
absence of flow in a detector containing suspended 
particles. He has constructed such an instrument, 
the “centrifugal manometer,” for use in air and has 
found it to perform very satisfactorily. 

A manometer if this type has several advantages 
for measurement of a pressure difference between 
two points in a given fluid, providing only that the 
difference is nearly steady. It contains but a single 
fluid, the same as that given. Consequently, with 
no fluid in contact with another, there is no meniscus 
or bubble and no column of liquid to be measured. 
Any change in the unknown pressure difference 
accelerates the flow in the detector almost immedi- 
ately, and there is no lag associated with the filling 
and emptying of reservoirs. This advantage is 
much greater for measurements in water than in air, 
as such a lag increases with the inertia and viscosity 
of the fluid. Since the operation of the manometer 
is independent of the action of gravity, it does not 
need to be level and does not require a steady plat- 
form. ‘The principle of operation is simple and, 
with careful design, the performance can be predicted 
with accuracy sufficient to make a calibration un- 
necessary. Finally, since the known pressure dif- 
ference created within the instrument does not 
require thermal expansion or compressibility of the 
fluid, such a manometer is suitable for use in a liquid. 

The primary effort described here is the design 
and test of a centrifugal manometer adapted for 
use in water. This adaptation involves, on the one 
hand, such practical problems as the choice of 
particles to be suspended and the control of bearing 
leakage, bubble formation, and temperature varia- 

J. F, Kemp, ¢ 


entrifugal mancmeter, J, Basic Eng., pp. 241-347 (Sept, 1959), 


tions. On the other 


hand is the question of to 
what extent the 


predictability of performance, 
found by Kemp, will be repeated in a_ pressure 
generator of quite different size and design, in 
particular, with a rotor of much smaller diameter. 

Secondly, as will be seen, the concept of the 
centrifugal manometer is easily generalized to that 
of the “rotational”? manometer, a class which 
includes many additional possible designs. Since 
these designs differ only in their pressure generators, 
retaining the same system of rotor drive and flow 
detection, it was found convenient to test one and 
compare its performance with that of the centrifugal 
manometer. 

Finally, the sensitivity of manometers of a rota- 
tional type is investigated to determine upon what 
factors it depends and how it can be increased. 


2. General Description 


2.1. Principle of Rotational Manometers 


The general principle of the rotational manometer 
is easily derived. It will be found most convenient 
to treat the centrifugal manometer and the alterna- 
tive ‘‘viscous’’ manometer later as special cases. 

If a solid of revolution is rotated about its axis 
in a cavity filled with a fluid, the pressure varies 
from point to point throughout the fluid. Points 
within the solid may be included if they are accessible 
through suitable passageways. It is assumed that 
the rate of rotation is constant and that the pressure 
at all points is steady. Denoting the pressure 
difference between any two selected fixed points 
by Ap, a characteristic length of the system by L, 
the density and kinematic viscosity of the (incom- 
pressible) fluid by p and », and the rate of rotation 
by Q, a simple dimensional analysis shows that 


bars) a) 
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or alternately, 
Ap _ 
TL? pQ? ~ 


The function f, or F, 
ment and, in 
analysis. 

If now, in a separate fluid of the same properties, 
a pressure difference Ap’ exists between two given 
points (as say on a pitot-static tube) and these 
points are connected by suitable conduits to the 
previously mentioned points in the cavity, a current 
will, in general, flow through the conduits. But 
if the rate of rotation is adjusted such that Ap 
equals Ap’, the current will stop. Thus, if f is 
known and Ap’ is unknown, the cavity and rotating 
solid, with a motor and drive to rotate it, and with 
a sensing element to detect conduit flow, constitute 
a rotational manometer which determines Ap’ by 
eq (1) whenever the flow is stopped. Conversely, 
if Ap’ can be given known values, the same apparatus 
can be used to determine f. 


can be determined by experi- 
exceptional cases, by theoretical 


3. Water Manometer System 


The components and general scheme of the water 
rotational manometer are shown in the photograph 
in figure 1. The motor and drive at the right turn 
the rotor within the cavity housing seen at the left 
center. The pressure generator shown is that of the 
centrifugal manometer. The flow detector is at the 
far left. In operation, the conduit at the left of the 
three open nipples shown in front of the cavity-rotor 
unit is connected to the external source of lower 
pressure. It then passes through the detector and 
into the rotor or cavity. The circuit is completed by 
the conduit in the center which connects the cavity 
to the external source of higher pressure. 

It is necessary to prevent any leak or bubble for- 
mation in the system as these cause spurious move- 
ment through ~ flow detector. The conduit at the 
right in figure which is not part of the circuit, 
connects an exter ele water source to a point near the 
outside end of the drive shaft bearing to make the 
pressure at this point nearly equal to that in the 
cavity. Any leak out of the bearing, when this 


FicurE 1. The water rotational manometer: 





| pressure is greater than atmospheric, is then drawn 
| from the external source only. If the cavity pressure 
is less than atmospheric, the bearing conduit can be 
used to keep air from entering the system. This 
method of leak prevention, rather than the use of 
some form of packing gland, was chosen to keep the 
drive shaft bearing friction negligible compared with 
the friction between the cones in the drive system 
described below. To prevent bubble formation, 
water used to fill the system was previously degassed 
by boiling. Since the volume flow during an obser- 
vation was extremely small, this water remained in 
the system practically indefinitely. 

In any rotational manometer it is necessary that 
the rate of rotation can be given any value within 
the operating range and there be held constant. For 
these purposes, a synchronous motor was used with 
a drive which included a system of gear selections and 
a pair of slightly tapered cones for continuous speed 
control. Torque was transmitted from one cone to 
the other by a weighted wheel with a narrow rubber 
(O-ring) perimeter. These components can be seen 
in figure 1. 

In the calibrations, the period 7=27/Q, rather than 
Q, was measured directly. This was done by timing, 
with a stop watch, a convenient number of cycles of 
the rotor shaft, the number chosen so that the timing 
lasted for around 1 min. The readings were found to 
be repeatable to within about 0.3 sec, that is, to 
about % percent, except at the highest speeds when 


| some slip between the cones became noticeable. In 


the calibrations, values of 7 
down to 0.35 sec. 

The flow detector was similar to that described in 
footnote [1], except that there the suspended particles 
were oil droplets. A narrow shaft of light from a 
vertical filament passed through a glass- walled section 
of the conduit and illuminated small particles in the 
water. Light rays emanating from the particles 
were then deflected upward by a prism and observed 
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Figure 2. The design of the flow detector in the vicinity of the 


suspended particles. 


(The glass walls are 1/16-in, apart.) 





through a low power (about 70) microscope. The 
angle between the shaft of light and the normal to 
the walls was sufficient (about 45°) to provide a dark 
background for the illuminated particles. These 
were of fine clay which, it was found, usually con- 
tained miscellaneous organic material. The particles 
could be replenished through a hole in the top of the 
section from a syringe serving as a reservoir. The 
smaller particles, whether clay or organic (it was 
seldom obvious), showed no tendency to settle and 
the reservoir was needed only occasionally. 

The design of the flow detector in the vicinity of 
the suspended particles is shown schematically in 
figure 2. The circle, drawn on the side view, shows 
the approximate field of vision, seen through the 
microscope, and the two vertical lines within it 
represent thin wires which were placed in the focal 
plane within the barrel of the microscope. 


4. Calibration Systems 


Known pressure differences were provided by two 
calibration systems. The first consisted of two 
bottles of known diameters and partly filled with 
water, with a syringe, connected to a micrometer 
which injected a known volume of water into one 
of the bottles, thereby increasing the depth of water 
by a known amount. The useful range in Ap’ for 
the bottles was, roughly, between 100 dynes em~? and 
1 dyne cm~*. In a second system, the ‘‘Poiseuille 
apparatus,” the pressure difference between two 
taps on the side of a long circular pipe, with carefully 
measured diameter, was calculated from the known 
rate of flow. This flow, from one vessel to another, 
passing through the pipe and any one of several 
capillaries, was previously determined in terms of 
the kinematic viscosity v and the differences in water 
level in the These differences were sub- 
stantial, between 5 and 50 em. This apparatus 
could provide any desired value of Ap’ less than 
about 20 dynes cm~*. In figures 4 and 5, the 
calibrations with each system are differentiated. 

The probably error in the pressure difference 
calculated with the Poiseuille apparatus can be 
estimated in terms of errors in the separate measure- 
ments upon which the calculation depends. With 
the diameter of the pipe and the distance between 
pressure taps given by d and /, and with the molecular 
viscosity of the liquid and the volume rate of flow 
given by uw and gq, the expression for the pressure 
difference 


vessels. 


a\4 
Ap > (5) luq 
shows that 
dA p 
Ap 


=a ey le (3) 
et. 2 a 


where 6 denotes an error in 
diameter of the section of 
taps was determined with 
Of over a dozen readings, 
mean (0.3031 in.) by more 


measurement. The 
the pipe containing the 
a slightly tapered plug. 
none deviated from the 
than 0.10 percent which 
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is taken as 


es om - ee 

7" The error 6 is taken to be the 
a | 

diameter of the pressure taps (0.74 mm) so that, 

| dl 

with 7=10.0 em, |7|\=0.74 percent. The error in 
viscosity is proportional to an error in temperature 
taken to be 0.1 °C, the smallest division on the 
thermometer. At room temperature (25 °C) this 


6 
r| =0.23 percent. 


Finally, in the calibration 
of the Poiseuille apparatus, the difference in water 
level in the vessels was plotted as a function of time. 
The derivatives of these curves, taken at numerous 
points, were than plotted to give qg, the volume rate 
of flow. On the basis of the scatter in these plots, 
ldq| . ss 

“t) is taken as 0.5 percent. 
q 

' 


Taking as the probable 


error in Ap the square root of the sum of the squares 
of the terms on the right side of eq (3), it is found 


that 
(=F) 
Ap 

If the fixed error due to the first two terms is disre- 
garded, the variable error due to the last two terms 
is found to be 0.55 percent, which is a measure of the 
degree of scatter in the manometer calibrations 
attributable to the Poiseuille apparatus. 

The accuracy of the bottle apparatus is difficult 
to estimate since it is limited by the regularity of the 


=1.10%. 


prob 


| shape of the menisci in the bottles and by the steadi- 


ness of the supporting platform. However, that the 
probable percentage error in the bottle apparatus in 
its useful range is at least no larger than that of the 
Poiseuille apparatus may be inferred from a com- 
parison of the degree of scatter associated with each 
instrument in the region where they overlap, as 
found in the calibration in figure 4. 


5. Centrifugal Manometer 


The principal features of the pressure generator 
of the centrifugal manometer tested here are shown 
in figure 3. The cavity and rotating solid were 
concentric circular cylinders, separated by a constant 
gap, small compared with their length. The differ- 


ence in pressure was taken between the inner wall 


of the cavity and the axis which were connected by 
16 passages arranged in two rings. ‘The axial pres- 
sure was transmitted outward through the hollow 
bearing to the left, while the pressure at the cavity 
wall was transmitted through four passages to a 
collecting ring. The hollow bearing at the left of 
the diagram and the drive shaft, which was soldered 
into the rotor, were of stainless steel; all other parts 
were of brass. 

The diameters of the cavity and rotor were 1.746 
and 1.581 em (compared with 22.8 cm in Kemp’s 

















FiGURE 3. Design of the cavity and rotor of the centrifugal 


manometer, 
Che diameter of the cavity is 1.75 


manometer) leaving a gap of 0.825 mm, about one- 
tenth of the radius of the rotor. Their lengths were 
7.91 em and 7.80 cm, almost one hundred times the 
gap between them, and leaving at each end a gap ol 
0.55mm. Orifices in the rotor were of 0.74 mm diam. 
for a length of 1.20 mm thereafter increasing to 1.60 
mm diam. Those in the cavity walls were of 
approximately the same design. The diameter of 
the bearings was 0.476 cm (%j¢ in.), and the inner 
diameter of the hollow bearing was 0.325 em. 

As shown by the analysis in footnote 1, for a 
centrifugal manometer long enough to ignore end 
effects, and with laminar flow in the annular cavity, 


eq (2) takes the simple form: 


<1 ==K, a constant. 
L? pQ? 

In this case the pressure difference is unaffected by 
the viscosity. Taking as Z the radius of the rotor, 
the analysis gives K as a function of the ratio of the 
gap between cylinders to the radius. As the gap is 
reduced to zero, K becomes equal to 1/2. For the 
manometer here described, A takes the value 0.533. 
However the pressure along the rotor axis may differ 
from that within the connecting bearing conduit 
where rotation of the fluid does not penetrate. 
Assuming that the average pressure within the inner 
radius of the bearing conduit is the same in the rotor 
as in the conduit and substituting this average for 
the axial pressure, A is reduced to 0.5224. Intro- 
ducing 7=27/2 and the numerical value of Z in cm, 
and with A=0.5224, the predicted calibration is 


Ap 
p 


T?=12.89 cm’. (4) 
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FIGURE 4. Calibration of the centrifugal manometer, pressure 


COE fic tent against the pe riod of rotation. 


refer to calibration with the bottles and the Poiseuille 
pparatus respectively 


Figure 4 shows the calibration of the centrifugal 
manometer obtained by test and plotted according 
to eq (2). The abscissa, 10°v7’, is a substitute for 
[?Q/v. A constant value of the ordinate, of 12.55 
cm’, provides about the best fit for the data. This 
value is about 2}; percent less than that predicted by 
eq (4). To see more readily what pressure differ- 
ences the observations represent, the data have been 
replotted in figure 5 according to eq (1). Since for 
water at room temperature, 10°»~0.9 or, roughly, 
unity, the ordinate in figure 5 approximates the 
measured pressure difference in dynes em-*, Sim- 
larly, the abscissa in figures 4 and 5 approximates 
the period in seconds. The figures show that, under 
the conditions of the calibration, the manometer 
readings had a scatter of less than about +2 percent 
for Ap >1 dyne cm~ and of about 10 percent for 
Ap around 1/10 dyne em~*. As Ap is further 
reduced, the scatter increases rapidly. 

The deviation of the experimental calibration from 
the theoretical prediction is roughly of the same 
magnitude as the correction in A from 0.533. to 
0.5224 which was based upon an assumption. <A 
reduction in the inner radius of the bearing conduit 
reduces both the correction and the theoretical 
uncertainty and, it may be assumed, would also 
reduce the discrepancy between the observed and 
predicted calibrations. In Kemp’s design, where 
the ratio of the radius of the bearing conduit to the 
radius of the rotor was about \ of that here, no such 
discrepancy was apparent. 

An investigation made by G. I. Taylor ? shows 
that the simple laminar flow in the annulus of the 
present manometer breaks down when 10%v7' is 
reduced to about 0.166. This value, corresponding 
to 10~* S 


pv 


p cle sg : . 
5 ~457, is beyond the range of the calibra- 


2H. Schlichting, Boundary Layer Theory, 


Ist English ed., p. 359 
Hill Book Co., New York, 1954 
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Figure 5. Calibration of the centrifugal manometer, pressure difference against the period of rotation 


rhe symbols £ 


and © refer to calibration with the bottles and the Poiseuille apparatus respectively. 
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tions. Presumably, however, the change of regime 
would have little effect on the calibration curve, since 
Kemp’s data, which contain the point of inst: ibility 
within their range, show no such effect. 

There is no apparent upper limit to the straight 
line calibration curve in figure 5. However, at the 
higher rates of rotation a slipping of the drive 
between the cones began to appear. Thus an 
upward extension of range would require an improve- 
ment of the drive system. 

The flow through the test section was found to 
respond to any change in rotor speed without 
appreciable lag. This is because the liquid in the 
radial passages of the rotor moves almost as a rigid 
body, and its angular acceleration does not depend 
upon viscous action. 


6. Viscous Manometer 


In types of rotational manometers other than the 
centrifugal manometer the pressure difference de- 
pends in general, upon the viscosity. This is a 
practical disadvantage. Conceivably, however, some 
designs may be of practical interest because of the 
possibility of extreme simplicity in the configuration 
of cavity and rotor. 

The design of the pressure generator of the 
“‘viscous’? manometer tested here is shown in figure 
6. The cavity and rotor are circular cylinders 
with axes parallel but eccentric producing a varying 
gap and pressure along the perimeter. The points 
of pressure difference are both on the inner surface 
of the cavity at opposite ends of a diameter. Two 
configurations were examined. In one, type A, 
te, and in the other, 


the pressure taps are at @ 


— 


2 


ROTOR 
AXIS 


CAVITY 
AXIS 
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FicureE 6, Design of the cavity and rotor of the viscous 


manometer. 


ivity i 


| 





type B the taps are at 6=0,7. The design is very 
simple with no conduits in the rotor. All parts 
were made of brass except the shaft which was of 
stainless steel. The temperature was taken in a 
brass well fitted along the outside of the cavity 
casing. 

The dixmeters of the cavity and rotor were 1% 
in. and 1% in. The rotor axis was displaced ' 
in. from that of the cavity so that the gap between 
cylinders varied from %4 to 1%, in. The lengths of 
the cavity and rotor were 1'%¢ and 1% in. leaving 
a gap of M4 in. at each end. The design of the 
orifices in the cavity wall was similar to that already 
described. The shaft and bearing were \ in. in 
diameter. 

The principle of operation is simple. 
the distance and velocity along a streamline 
s and u, the equation of motion for flow in 
annulus can be written, approximately, 


ona Pte) 


The term on the left is the resultant force per unit 
volume along s due to viscous action and is positive 
when the velocity profile is concave in the direction 
of flow. Because of continuity of flow through the 


Denoting 
by 
the 


Oru 
Y Oy 


Op 
Os 


annulus, the profiles are concave for 


and convex otherwise. 


averages between the rotor and the 


approximately, Denoting 


by p and u?* the 


u? 


cavity wall for given 6, P+ p-~ is @ maximum near 


and a minimum near @ And since at 


these two points the values of uw? may be assumed 


equal, there difference in pressure between 
them which is proportional to the molecular vis- 


is a 


a ’ has about the 


cosity. On the other hand, 
0 and @ 
u? is not the same, there is a difference in pressure 
between these two points approximately equal to 


same value at 0 x and since, by continuity, 


the difference in p = 


The gap and eccentricity were made sufficiently 
large to insure that imperfections of construction 
would be of no consequence, but otherwise theic 
dimensions were selected to produce a range of 
pressure about the same as that of the centrifugal 
manometer for the same range of rotational speeds. 
This selection was based upon an analysis of the 
pressure variation around the cavity wall. Because 
of its approximate nature, this analysis could hardly 
be very accurate and it is not included here. How- 
ever the calibrations thus predicted for types A 
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Fiaure 7. Calibration of the viscous manometer, 


types 


The symbols A and e refer to 


and B are shown in figure 7. case the 


while in 


In the first 
pressure difference is proportional to yQ 
the second it is hig seems to pQ?, 

The outstanding feature of the calibration data 
plotted in figure 7 is the high degree of 
compared with that of figures 4 and 5. The proven 
reliability of the sensing apparatus in calibrating 
the centrifugal manometer shows that, primarily, 
this scatter is not a random error in the observation 


scatter 


O° VT 


A and B, pressure difference against the pe riod of rotation. 


types A and B respectively. 


of Ap, but rather is due to variations in Ap itself. 
That is, the initial assumption of steady pressure is 


violated. Indeed, for a given v7’, an observed value 
of Ap/pv? was frequently repeated several times 
consecutively, with small variation, and then ob- 
served to shift to a new value. <A nonrandom 
variability in Ap is apparent for case A in the lower 
portion of figure 7 where the data fall on two dis- 
tinct curves. In a first calibration only values on 
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the lower curve were found, but in a second, some 
weeks later, values on both curves appeared. A 
general increase in the scatter of case B as 102r7' 
is reduced below about 0.9, and a tendency for the 
slope of case A to steepen below this same value, 
suggest some change in regime at this point, pos- 
sibly a transition to turbulent flow. At least some 
of the variation for higher values of 10°77 might 
derive from cellular types of flow related to those 
between concentric cylinders and investigated by 
G. 1. Taylor. But however interpreted, the scatter 
in figure 7 demonstrates the failure of a 
manometer of the present design. 

It must also be noted, that because any change 
in the pattern of flow in the cavity requires the 
action of viscous forces, the rotor response of this 
tvpe of manometer is relatively slow. In the range 
of the smaller pressure differences, the establish- 
ment of a near steady state required from one to 
several minutes. 

There may _ be designs for a viscous 
manometer which would prove more successful. In 
the present case a reduction of the gap might reduce 
the scatter to an acceptable degree. Alternative 
designs of simple geometry construction are 
readily conceived but not investigated 
here. 


viscous 


other 


and 


these are 


7. Manometer Sensitivity 


7.1. Effect of Circuit Resistance and Test Section 
Design 


If the pressure Ap is not precisely adjusted to 
Ap’, the imbalance, say 6p, is accompanied by a 
volume rate of flow through the system, say 6g, 
and in the illuminated cross section, say of area A, 
there is an average velocity of flow év equal to 
5q/A. 

Taking 6v to be inversely proportional to the time, 
t, that a particle remains within the two wire lines 
shown in figure 2, and determining 6p by either of 
the two calibration systems, it was found that 6v 
increased, approximately, in proportion to 6p for 
any given combination of manometer and calibration 
system. Further investigation, with the centrifugal 
manometer, confirmed that 6v/ép remained un- 
affected by Q. 

Thus the flow through the system due to the 
pressure imbalance can be taken as laminar, and 


1 6 
P AR (5) 


u Ov 


where Ff defines the resistance of the system to the 
flow 6g, through it. Since the smallest pressure 
which can be measured depends upon the smallest 
velocity which can be detected, the sensitivity of a 
manometer system with 6v/dp and, to 
maximize sensitivity, AR must be minimized. 
The design of the conduit in the vicinity of the 
illuminated particles, the test section, determines A 
and greatly affects the sensitivity of a manometer 
system. Denoting by r that part of the total 


increases 


resistance 2? which is affected by the test section 
design, and by 7 the remainder, which may be 
assumed irreducible, there is an optimum design 
which minimizes A (7+79). Any change in 7, as by 
replacing one pitot-static tube with another, re- 
quires a new optimum design. In these investiga- 
tions, different values of 7) correspond to different 
combinations of manometer and calibration system. 

One practical approach to an optimum variable 
detector design is to vary the depth of the test 
section. In the case tested here, the glass walls of 
the test section were plane and parallel and fixed 
at Ms in. apart. Between the walls were fitted a 
fixed roof and a moveable floor of brass. Both were 
curved to give a minimum cross section at the plane 
of illumination. These features are indicated in 
figure 2. Moving the floor up or down changed 
both A and r and, consequently, the sensitivity. 

Figure 8 shows the observed variation of the 
sensitivity with depth of test section, for three 
values of ro. The ordinate, tép/u, is approximately 
proportional to 6p/yév, that is, inversely proportional 
to the sensitivity. The abscissa, is the depth of 
the test section. The solid lines (1, 2, 3) represent 
A (r+ro) with values of ro in the ratios 1:6:13 and 
with r a funetion of 2 chosen best to fit the data. 
While this fit is not very close, the data show that, 
with increasing 7, the sensitivity decreases and the 
value of z for ‘optimum design’ also decreases. 

The test section depth was not made adjustable 
until after the manometer calibrations were com- 
pleted. During the calibrations it was constant at 
0.0625 in. (equal to the width of the section), a 
value selected from a preliminary calculation. As 
seen in figure 8, varying degrees of sensitivity were 
thereby sacrificed. 





Figure 8. Variation of sensitivity with depth of test section z, 


and resistance to flow, To. 


On the curves (1, 2, 3) values of ro are in the ratio 1:6:13. 
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It is noted that the resistance of any element in 
the circuit can be determined easily by connecting 
the ends to open vessels of water of differing head and 
measuring the volume flow in a given time. 


7.2. Drift Errors and Corrections 


Experience showed that the factor limiting the ac- 
curacy of the calibrations was not, in general, an in- 
ability to detect év but rather drifts through the 
system caused by variations in temperature. Pre- 
sumably these drifts arise primarily from temperature 
induced variations in the density of the water from 
place to place in the system, and perhaps in part from 
differing rates of expansion of the solid members of 
the system. In any case, before calibration, with Q 
and Ap equal to zero, any drift observed through the 
microscope could be eliminated. This was done by 
adjusting the position of the end of a warm copper 
rod between the vertical elements of a U-shaped 
section of the conduit. The rod was heated by the 
light source which was near its other end. However, 
immediately after a reading with Ap again zero, a 
drift was often observed. This drift implied an error 
6p to be added to or subtracted from Ap, and whose 
magnitude was given by the known value of tip/u. 
For values of Ap above some limit the drift correc- 
tions were inconsequential but became critical as Ap 
was reduced to the average magnitude of Op. In the 
case of the centrifugal manometer these limits were 
roughly 1 dyne em—* and 1/10 dyne em—*. The ma- 
nometer system was partly insulated from the light 
source and was ventilated by a fan. No attempt was 
made to measure the temperature field around the 
rotor housing except to note that temperature differ- 
ences were too small to be reliably determined by a 
thermometer with increments of 0.1° C. However 
the addition of the fan was found to reduce the drift 
corrections markedly and further attempts at tem- 
perature control may hold promise. 

The probable significance of a drift correction in- 

creases with the time required for the reading. Un- 
fortunately, for the manometers. tested, there 
occurred a slight pulsing of the observed particles 
with each turn of the rotor, presumably due to ir- 
regularities in bearing friction combined with the 
possibility of slip between the cones. As a result, 
verification of no flow required an observation lasting 
through several cycles, the period of which increased 
with decreasing Ap. This pulsing, requiring longer 
eadings which allow the temperature drift time to 
develop, was found to be the single factor most 
damaging to the convenience and accuracy of the 
centrifugal manometer. 


8. Conclusions 


There are many possible designs for a rotational 
manometer, that is, a meniscus free instrument which 
balances an unknown pressure difference against the 
known action of a rotating element. One design, the 
centrifugal manometer, which had been previously 


] 


| of 


investigated and tested successfully in air, has been 
found here to be also suitable for use in water. This 
type of rotational manometer has several distinct 
advantages. The calibration is of simple form and, 
with careful design and construction, can be pre- 
dicted, that is, no calibration is necessary. The rotor 
response of the instrument is rapid, and the pressure 
difference produced by the rotation is independent of 
the molecular viscosity, so the measurement of tem- 
perature is unnecessary. The single alternate design 
here investigated, the ‘“‘viscous’’ manometer, had 
none of these advantages and, because of variability 
in the calibrations, proved unsuitable. The flow pat- 
terns in the cavity upon which the somewhat variable 
pressure differences depend are probably complex 
and dominated by viscous forces. While improve- 
ments are possible, it seems apparent that the centrif- 
ugal manometer is the most satisfactory of rotational 
manometers. 

Under the conditions of calibration, the centrifugal 
manometer gave readings repeatable to within about 
+2 percent for Ap dyne cm~? and to about +10 
percent for Ap around 1/10 dyne em~*. In nearly all 
the readings for Ap<10 dynes cm~? the Poiseuille 
calibrating apparatus was used. Because the two 
orifices in the side of the Poiseuille flow tube were 
small, the resistance of this apparatus was probably 
larger than that of many pitot-static tubes of likely 
dimensions and design. Therefore the performance 
under calibration was reasonably close to that which 
might be expected under normal operation. 

The sensitivity of the centrifugal manometer 
tested in water compares favorably with that of 
micromanometers involving a fluid surface and is 
similar to that attained by most meniscus free 
manometers used in air. The smallest detectable 
pressure differences were roughly twice those de- 
tectable by Kemp’s instrument in air. Comparison 
manometers for use in different media by the 
pressure differences they can detect is somewhat 
arbitrary. An alternate criterion is the percentage 
error in measuring, with a pitot-static tube, the 
velocity V at a given point in a conduit for a given 
Reynolds number. This error is proportional to the 
error in the pressure difference 1/2pV? which, at 
room temperature, is about three times larger in 
water than in air. 

The sensitivity of a rotational manometer, or the 
error to which it is liable, depends upon the resistance 
of the circuit in which it operates. An increase of 
the resistance of any element decreases the sensi- 
tivity. Therefore any anticipation of probable 
error in actual operation must rest not only upon 
the results of calibration but must also compare the 
resistance of the calibrating apparatus with that of 
the elements which are to replac eit. A manometer’s 
sensitivity is also dependent upon the design and 
dimensions of the test section in the flow detector. 
Ideally, the optimum design differs for each value of 
the circuit resistance. One approach to a variable 
design has been found practical. 


(Paper 66C4-112) 
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Study of Gypsum Plasters Exposed to Fire 


J. V. Ryan 
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1962) 


Experimental results are reported from a study in which small gypsum plaster specimens 
were exposed to controlled fires similar to those to which large building elements have been 


subjected in tests by a recognized standard method. 


either structural load or restraint. 
perature rise criterion, 


The small specimens were tested without 


Their fire performances were judged on time to a tem- 
The results were analyzed to determine the effects on fire performance 


of variation of mix ratio, aggregate type and density, duration of aging, and relative humidity 


of the ambient in which the specimens were aged. 
ranges normally used, had little effect, if any, on fire performance. 


Mix ratio and aggregate density, over the 
The times to temperature 


rise for perlite and vermiculite aggregate plasters were essentially equal, but shorter times 


were observed for sanded plasters. 


Duration of aging and relative humidity of the ambient 


had a significant effect only for short aging periods and very high relative humidity conditions. 
Kstimates of thermal properties of gypsum plasters at elevated temperatures were 


derived from the data. 
building elements involving gypsum plaster. 


1. Introduction 


Prehistoric man probably plastered the interior of 
his stick-and-grass hut with mud. Gypsum plaster,’ 
similar to that of today, was used in Egyptian 
basape a and tombs at least 4,000 years ago; white 
ime plaster was used in Greece 3,000 years ago, and 
by the Romans later. It is not clear that any of the 
advanced civilizations of Egypt, Greece, or Rome 
recognized the value of plaster as a fire protective 
material. In more recent times, official recognition 
included an Ordinance by King John, following the 
London fire of 1212, in which plastered walls were 
required in certain buildings [1].2~ Although plaster 
has been used extensively for fire protection through 
the years, the particular application was based more 
on accumulated experience than on the results of 
planned investigations, until comparatively recent 
times. Even today much of the data available on 
the fire protection of structures with plaster was 
derived from isolated tests conducted for “approval” 
purposes, rather than from systematic study. 

A number of cementitious materials are employed 
in the plastering of building elements. Among these 
are gypsum, portland cement, lime, and various 
combinations. Gypsum neat plaster consists of not 
less than 66 percent of hemihydrate calcium sulfate 
(CaSO,-%H,O) plus other materials added at the 
factory to control time of set and working qualities. 
When mixed with water, the gypsum changes to 
CaSQO,-2H,O, making the hard plaster of wide use. 
Under fire exposure, the fully hydrated gypsum is 
calcined, or changed back to hemihydrate, at about 
260 °F and to anhydrous CaSO, at about 350 °F. 

1 In accordance with established technical and practical usage, the term plaster 
is used in this paper for both the dry, powdery material sold on the market and the 
hardened end product resulting from chemical reaction with water, with or with- 
out aggregate. 

2 Figures in brackets indicate the literature references at the end of this paper. 
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These estimates should be useful for predicting fire endurances of 


Both these phase changes are endothermic. Also, 
the chemically combined water released, representing 
21 percent by weight of the original material, plus 
additional free water ordinarily present is vaporized. 
Consequently, great quantities of heat are absorbed 
without causing temperature increases during the 
calcination of gypsum. This fact makes gypsum 
plaster an important material for the fire protection 
of buildings. 

The present study on gypsum plasters consisted of 
an experimental program on the fire endurances of 
small plaster specimens plus determinations of 
thermal conductivity and other physical properties. 
The plasters represented a variety of aggregate and 
mix combinations. The data obtained were ana- 
lyzed to make comparisons among the plasters in 
terms of fire endurance and to make estimates of 
their thermal properties at elevated temperatures, 
both useful for predicting the fire endurance of 
building elements involving gypsum plaster. 


2. Investigation 


Several factors affect the behavior of a plastered 
building element when exposed to fire. Some, such 
as method of attachment between plaster and struc- 
tural elements, are characteristic of the assembly. 
Others, such as mix proportions and aggregate type, 
are characteristic of the plaster itself. The former 
were eliminated as far as possible for the experiments 
reported here. 


2.1. Factors Varied 


The investigation was designed to provide com- 
parative data on the effects of four factors, as well 
as the interactions among them. The factors were: 

A. Aggregate. Three different aggregates were 
used: perlite, vermiculite, and sand. In addition, 
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the perlites used were in three density ranges, 
designated as high, intermediate, and low. 

B. Mir. The mixes chosen included those in com- 
mon use plus some leaner and some richer in gypsum 
content. 

co Length of aging. The length of aging period 
prior to the fire test was varied from 28 to 103 days. 

D. Humidity. Each specimen was aged at room 
temperature in air at either high or intermediate 
relative humidity. 


2.2. Specimens 


Three sizes of specimens were prepared: those for 
fire tests, those for thermal conductivity determina- 
tions, and those for compressive strength tests. 
Each fire-test specimen consisted of approximately 
1 in. of plaster applied to a 24-in. square base of 
metal lath. Each thermal-conductivity specimen 


was a 1- by 8- by S8-in. slab of plaster; each com- 


pressive-strength specimen, a 2-in. cube of plaster. 
The specimens were made in sets, each set consisting 


of two fire-test specimens, two thermal-conductivity 
specimens, and either three or six cubes for com- 
pressive strength determinations. All the specimens 
except two sets were prepared by the same experi- 
enced craftsman, whose work was excellent. The 
forms for specimens are shown in figure 1. 


a. Materials 


The plasters were all regular commercial gypsum 
plasters and complied with the requirements of 
ASTM (C-—2s; that for all but one set of the sanded 
mixes was unfibered; those for the lightweight aggre- 
gate mixes were extra fibered (sisal) and formulated 
for use with lightweight aggregates. The program 
was started with a supply of the latter type all from 
a single production lot at one plant. This proved 


Figure 1. Specimen 
with two lath-and-channe 
P-1n. cubes: 


forms: upper, for fire test specimens, 
lower left, for 


8- by 8- by 1-in. specimens. 


assemblies in place; 
lower right, fo 
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TABLE 1. Data on aggregates 


Perlite 
Vermic 
ulite 


Inter Hleavy 


mediate 


DENSITY, lb/ft 


Maximum 12. 90 
Mean 5 12. 20 
Minimum 7. 8 11. 50 


ANALYSIS 


cumulated percent retained 
By 
weight 
By volume 


> 


19.2 21.6 
10.6 56 67.3 


83.9 4 92. 7 
95. § 4 OS. 5 


100.0 100.0 100, 0 


FINENI 


» Dry weight of 1 ft? of damp sand, average value 

to be an inadequate supply and a supplementary 
supply was obtained. The plaster for lightweight 
aggregates proved unsatisfactory with sand aggre- 
gate, so the unfibered gypsum plaster was obtained. 
The hemihydrate content of the initial supply was 
88.4 percent, that of the supplementary supply 85.7 
percent. The unfibered gypsum plaster for use with 
sand was not analyzed. 

The aggregates were sand, perlite, and vermicu- 
lite. The sand was from the White Marsh pits. 
This sand, known to be highly siliceous, was received 
with a minor amount of small gravel. It was sieved 
to remove this gravel, placed in fiber drums, and 
stored indoors. The moisture coatent was 2.1 per- 
cent upon receipt. One sei of specimens was made 
with the sand at the original moisture content using 
extra fibered plaster. The moisture content of the 
sand was less than 0.2 percent by the time the un- 
fibered gypsum plaster had been obtained for the 
remainder of the specimens. 

The perlite aggregates from an ore deposit near 
Grants, N. Mex., were provided in three density 
ranges: low, intermediate, and high, varying from 8 
to 13 Ib/ft®. The intermediate density perlite was 
reported as representative of normal production. 
The perlites were received in fiber drums and stored 
indoors. The vermiculite was obtained in paper bags 
from a local building supply firm and was stored 
indoors. Densities, sieve analyses, and fineness mod- 
uli for the various aggregates are given in table 1. 

The water was drawn from the municipal supply 
into clean buckets for weighing and dumped into 
the mixer. The metal lath and furring channels were 
standard commercial products. 





b. Mixes 


The proportions of the plasters with perlite or 
vermiculite are expressed as bags (nominally 99.25 
lb, about 1.5 ft *) of gypsum to volume of aggregate 
in cubic feet; those of the sanded plasters are 
expressed as weight of gypsum to the dry weight of 
the sand. The amount of water used was determined | 
by the ‘workability’ as judged by the plasterer. 
The various mixes prepared are listed in table 2. 

All constituents of each mix were weighed. The 
usual charge for a batch of lightweight aggregate 
plaster consisted of one cubic foot of aggregate and 
the proper weight of gypsum based on nominal 
99.25 lb per bag. Each batch was mixed 24% min 
after the last water was added. Data on the 
mean mix proportions are given in table 2. 


c. Preparation of Specimens 
The plaster bases, details of which are shown in 
figure 2, for two fire-test specimens were prepared 
with temporary nonabsorbent backing behind and 
braced against the lath. The backing produced 
smooth rear surface, thereby facilitating reliable 
average thickness measurements that could not have 
been obtained with the usual ‘‘keyed”’ rear surface. 
The two plaster bases were placed in a frame, shown 
with the forms for other specimens in figure 1, 
which was mounted in a nearly vertical position. 
The plaster, aggregate, and water were mixed in 
a 4 ft® capacity commercial plaster mixer. The 
resulting wet mixture was applied to the metal lath 
base by trowel, the surface struck off with a metal 
‘‘featheredge” and, after partial set had developed, 


TABLE 2, 


Based on mean values for m 


Perlite 


1:0,10 
1:0.20 
1:0.24 
1:0,21 
1:0,32 
1:0,.21, 
1;0.16, 
1:0.24, 
1:0.38 


Vermiculite 


1:0.18 


1:2 
1:4 


« Mixed units 


» L=low-density perlite 


iterials used in preparing 


bags of gypsum to cu ft of lightweight aggregates; by weight, for sand aggregate. 
, 1 =high-density perlite, 


was again troweled to produce as smooth a surface 
as practicable. Exc ept for specimens of 1:2, 1:3 
mixes (E, E—L, E—H in tables 2 and 3) that required 
two-coat application, the full thickness of the 
lightweight aggregate plaster specimens was _ built 
up in a single continuous application. The sanded 
clan being considerably heavier, tended to slump 
so that two-coat applications, both of the same 
mix, were necessary. In all two-coat applications, 
the first coat was scratched with a wire brush. The 
second coat was not applied until the following day, 
although the first usually set within a few hours. 

The thermal-conductivity specimens were prepared 
by placing in each form an amount of plaster esti- 
mated at about half that required to fill the form, 
with an attached vibrator in operation. More 
plaster was added, care being taken to eliminate, so 
far as possible, air pockets. The top surface was 
struck off and troweled. 

The compressive-strength specimens were made by 
partially filling the forms with plaster, attempting 
to eliminate air pockets, filling, and striking off the 
top surface. 


d. Aging 


As soon as the troweling of the specimens was 
finished, they were covered with wet burlap. The 
next day they were removed from the forms and 
placed in the conditioning rooms, at which time the 
temporary backing was removed from the fire-test 
specimens. From each set, generally, one fire-test 
specimen and one thermal-conductivity specimen 
were placed in a high relative humidity box. The 
other fire-test and thermal-conductivity specimens 


Miz information 


from one to eight batches of the various mixes 


Water-to-gyp- 
sum ratio Dy 
weight 


Compressive 
strength 


By volume 


iggregates 


0.45 
59 
bb 
64 
.69 
.59, 0.80 
.59, 0.80 
60, 0.80 960, 490 
87 460 


320 
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THESE CHANNELS (2) 
REMOVED WHEN SPECIMEN WAS 
PLACED ON FURNACE 


+4 FURRING CHANNELS AND LATH 
FASTENED WITH SPOT WELDS OR 


SHEET METAL SCREWS 


SPECIMEN AND FRAME ASSEMBLY 





ASBESTOS MILLBOARD 


MINERAL WOOL 
INSULATION 


A 7 7 7 
"x 12"x24" Fl 
PLASTER 


x 


LOOSE FIRE BRICK 
BACKING 


BRICK SLAB 


23°X23"° FURNACE OPENING 


PLACEMENT OF SPECIMEN ON FURNACE 





X — THERMOCOUPLE LOCATIONS 


FIGuRE 2. 


The 


were stored in intermediate relative humidity. 
groups of three cubes were kept together, all in either 


high or intermediate relative humidity storage. The 
high-humidity box was operated at relative humidi- 
ties from 100 to 67 percent; the intermediate humi- 
dity ranged from 58 to 26 percent. All specimens 
were at room temperatures (generally 72+4 °F). 
The length of the aging period prior to fire tests was 
varied from 28 to 103 days. 





Details of fire lest specimen and placement on furnace. 


2.3. Fire Tests 


The fire tests were conducted using a small, gas- 
fired furnace on which the specimen was laid hori- 
zontally (see fig. 3). A fire brick slab was placed 
over the specimen and the edges of the specimen- 
slab assembly were insulated against lateral heat flow 
(see fig. 2). The arrangement of plaster specimen 
with fire brick slab above and a %-in. airspace be- 
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TaBLE 3. Thermal properties at 117 °F (47 °C) mean temperature 


Thermal conductivities and densities measured on 8-in. square specimens; specific heats computed by formulas in reference 


[3]; 
diffusivities a=k/pc. The computed ¢c and @ values for mean only. 


Conductivity, k, Density, p 
Plaster des Average (No. of spec Btu/hr ft? (°F /in. 
ignation * | weight loss imens 
in drying * 


Specific Diffusiv- 
ee m3 . heat, c ity, @ 
Mean Max, , “a , Min, 


1b/ft8 3tullb 4 ft2/hr 
93 27 0. 0098 
68. : , 2 0106 
54. 25 . 0108 


49. : . 256 0110 
52 257 0110 
50. 5 . 252 0109 
51 55 0105 
51 5 . 0107 


OLOS 


0110 
0113 
. 0163 
. 0152 


* See table 2, 
Drying at 140 °F prior to thermal conductivity determinations, 


tween was designed to facilitate experimental 
measurements rather than to simulate a construction 
type. 

The test procedure differed from the standard 
method in several respects. The plaster specimens’ 
area exposed to fire was slightly less than 4 ft?, and 
the specimens were ne ither structurally loaded nor 
restrained. The Standard Test Method, ASTM 
K-119 [2], requires that a minimum ceiling area of 
180 ft? be exposed and that the floor specimen be 
loaded. 

Also, the procedure differed from the standard in 
the way that the temperature of the upper surface 
of the fire brick slab was measured, as described in a 
later section. 

a. Furnace 


The furnace consisted of a cubical combustion 
chamber about 23 in. on a side and open at the top. 
The fuel was that from the public gas supply of 
Washington, D.C., and had an average heating 
value of 1025 to 1035 Btu/ft®. The gas was fed 
from six burners into the combustion chamber, pro- 
viding diffusion flames. The flow of gas was regu- 
lated by a remotely controlled valve. 

The furnace temperatures were measured by four 
chromel-alumel thermocouples encased in porcelain 
insulators and closed iron pipes. The pipes were 
placed so that the thermocouple junctions were 
approximately 2% in. below the specimen. The gas 
flow was regulated automatically so that the average 
furnace temperatures approximated those of the 
standard time-temperature curve given in ASTM 
E-119 [2]. A metal pan that could be raised and 
lowered was suspended above the furnace. The pan 
had water supply and drain designed to permit a flow 
of water while maintaining an average depth of 
about in. This pan, with the water flowing, pro- 
vided a nearly constant temperature receiver for the a <a : oe 
heat passing through the specimen. The furnace SORE ¥ rmape Sone Joe UNE TY « Bee tout. 
and water pan are shown in figure 3 





377 





b. Instrumentation 


In addition to those in the furnace, 24 gage (B&S) 
chromel-alumel thermocouples were cemented on 
the centers of the upper surface of the plaster and 
upper and lower surfaces of the fire brick slab placed 
above the plaster, as shown in figure 2. This slab 
was 24 in. long, 12 in. wide, and 1 in. thick, with a 
density of 50.2 lb/ft. A 6-in. wide by 24-in.-long 
piece of 's-in. asbestos millboard was placed on each 
side of the fire brick slab to complete the cover over 
the specimen. Thermocouples were attached to the 
upper and lower surfaces of one of the pieces of mill- 
board to permit comparison of temperatures for 
detection of lateral differences, if any. 

The thermocouple leads were connected to self- 
balancing potentiometers, and those from the fur- 
nace to automatic programing equipment for regu- 
lating the fuel supply. The program was intended 
to be the time-temperature curve defined in ASTM 
E-119. Very close agreement with this curve was 
achieved in the tests. 

The thermocouple bead on the upper surface of the 
fire brick slab was partially embedded in the slab and 
partially exposed to air. No pad, such as used in 
standard floor tests, was placed over this thermo- 
couple. 


2.4. Supplementary Tests 


Steady-state thermal conductivity determinations 
on the 8-in. square specimens were made in the 
National Bureau of Standards pentane hot-plate 
apparatus, at a mean plaster temperature in the range 
116 to 120 °F. The specimens were dried to constant 
weight at 140 °F immediately prior to measurement 
(maximum weight 2.0 percent; average 0.7 
percent). The data are given in table 3. 

The compressive strengths of 2-in. cubes of the 
plasters were determined in a_ universal testing 
machine and are given in table 2. The data on the 
aggregates are given in table 1. 


loss 
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3. Results 


The results obtained in the supplementary tests 
are given in tables 1, 2, and 3. Therefore, the fol- 
lowing comments are limited to the fire tests. 


3.1. General Behavior 


In most of the specimens, one or two transverse 
cracks were observed during the fire test but were not 
more than }42 in. wide and therefore were judged not 
of great importance. These cracks became much 
wider as a result of cooling after the tests; however, 
most of the specimens were intact upon removal from 
the furnace. The gypsum-sand plaster specimens 
were exceptions. ‘These specimens did not appear to 
crack more than the others during the fire exposure 
but did upon cooling, to the extent that sizable por- 
tions of plaster fell off some specimens during removal 
from the furnace. 


3.2. Data Reduction 


The performances of the various plasters in the 
fire tests were compared on the basis of the elapsed 
times at which a temperature rise of 250 °F above 
the initial temperature was observed on the top 
surface of the fire brick slab, although the tests 
were continued beyond these times. Typical time- 
temperature curves are shown in figure 4. The 
250 °F rise is the same as that taken as a limiting 
criterion in standard fire-endurance tests of floors 
and walls, but was measured without the felted 
asbestos pad required in the standard test method. 
These temperature-rise times were plotted against 
plaster thickness figs. 5 and 6). Previous 
experience had indicated the fire endurances of solid 
specimens were a power function of thickness [4]. 
Therefore, the data for the plasters were subjected 
to least-squares analyses for fits to power functions. 

The temperature-rise times for the individual speci- 
mens were then adjusted, or normalized, to a 
thickness of 1.0 in. by the power function. The 
resulting times, adjusted for thickness variations 
among the individual specimens, were taken as 
best estimates upon which to evaluate the effects of 
and interactions among the variables of mix, aggre- 
gate, and conditioning. 

The observed temperature data and known ther- 
mal properties of the fire brick slab were used to 
estimate values of the apparent thermal conduc- 
tivities of the plasters at elevated temperatures. 
Values of density were computed, knowing the densi- 
ties at room temperature and the temperatures at 
which gypsum phase changes occur. Values of 
specific heat were obtained from formulas and data 
given by Kelley [3]. From these values of k, p, 
and ¢, values of the thermal diffusivity, a=k/pe, 
were calculated. 


(see 


4. Discussion 


The results were analyzed and are discussed from 
the viewpoint of estimating the effects of the vari- 
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ables in the program: mix, aggregate, and condi- 
tioning. There was some variation in thickness of 
the plaster, and in the times to the limiting 250 °F 
rise condition; the discussion will compare the results 
in terms of the mentioned variables with considera- 
tions of the dimensional and thermal characteristics 
that influenced the results. 


4.1. Thickness 


The techniques employed for applying plaster in 
this study, representative of normal building prac- 
tice, lead to normal variation of thickness. Among 
the perlite plaster specimens for fire tests, the 
standard deviation was only 6.4 percent of the mean 
thickness. Each mean thickness was based on 25 
measurements. The temperature-rise times were 
plotted against thickness, as shown in figures 5 
and 6. The results derived by least-squares analysis 
represented reasonable fits for power-function rela- 
tionships between time and thickness. The expo- 
nents derived from the analyses were 1.53 for perlite 
plasters, 0.93 for vermiculite plasters, and 1.27 for 
sanded plasters. An exponent of 1.67 had been 
derived previously [4] from data for large-scale 
fire-test specimens over an appreciable range of 
thicknesses and a variety of building materials. 
Because of the fairly good agreement found with 
the value of 1.67, and the much greater thickness 
range upon which it was based, curves based on an 
exponent of 1.67 were fit to the data and are shown 
in figures 5 and -6; that for the perlite plasters is 
repeated in figure 6 to facilitate comparisons among 
the plasters 


4.2. Mixes 


Examination of the distribution of temperature- 
rise time versus thickness data plotted in figure 5 
indicated that there was no grouping along lines 
according to mix within the periite plasters. Statis- 
tical analysis of the data * confirmed that the main 
effect of mix was not significant. 

The results of large-scale fire tests of plaster pro- 
tected building elements [4] had indicated that varia- 
tion of the mix proportions of gypsum to aggregate, 
within the range investigated in the present study, 
lead to variation in the fire endurance among speci- 
mens equivalent in other details. Therefore, it 
appears that the effect previously observed and 
ascribed to mix was either (1) actually not related 
to mix, or (2) an interaction between mix and some 
other factor. Under the latter condition, several 
possibilities arise: (1) structural-thermal behavior 
of the plaster of sufficient magnitude to affect the 
results obtained with large ‘‘standard’’ specimens 
but not to be noticeable with specimens of the size 
used in this study; (2) the same possibility with 


3 The temperature rise of 250 °F used as an end point in this study is the same as 
one limiting criterion in standard fire tests of walls and floors, upon which fire 
ratings are based, but was not measured in this study under the standard condi- 
tions Also, no correlation has been established between small-scale tests, such 
as in this study, and large-scale standard fire tests. Consequently, the tempera- 
ture-rise times from this study do not provide an appropriate basis for ratings or 
approvals 
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regard to the plaster base and support system; or 
(3) restraint imposed on large specimens by the test 
arrangement or by different structural-thermal 
behavior among the components of the assembly. 

It was not within the scope of the present study 
to investigate structural behavior. However, gen- 
eral observations were made of the plaster during 
exposure and, for most specimens, only minor crack- 
ing was noted during the fire exposure. 

The assumption that mix proportions made a 
weak contribution to an interaction effect may be 
further explored by consideration of the data from 
the tests of sanded gypsum plasters. As seen in 
figure 6, these data fall significantly below those for 
the lightweight aggregate plasters. Although this 
is undoubtedly due in part to the effect of aggregate, 
as discussed in the next section, the sanded plasters 
represent a much different mix, in terms of the per- 
cent of dry gypsum by weight, than the lightweight 
plasters. The sanded plasters contained 20 and 
33.3 percent gypsum by weight, whereas the light- 
weight plasters contained from 72.5 to 91.3 percent 
by weight. The results for 100 percent or neat 
gypsum plaster were similar to those of the light- 
weight plasters. From this it appears that the effect 
of mix, if any, is within the experimental scatter 
except for very great changes in the percentage of 
gypsum. It did not seem practical to check this 
possibility with lightweight aggregate plaster. The 
required volumes of perlite to produce the 20 and 
33.5 percent gypsum cement by weight mixes ob- 
tained with sand would be about 40 ft* per bag of 
gypsum cement and 20 ft * per bag, respectively. 


4.3. Aggregate 


The study involved both variation as to type of 
aggregate and variation of density range within one 
type. 


a. Variation as to Type 


Examination of figure 6 indicates that the differ- 
ences of temperature-rise times between the perlite 
and vermiculite plasters were small, whereas that 
between either of them and the sanded plasters was 
comparatively large. Statistical analyses were used 
to evaluate the apparent differences. 

The average time normalized to 1.0-in. thickness 
for the vermiculite plasters aged in intermediate 
relative humidity was 61.1 min, whereas the cor- 
responding time for the equivalent perlite plaster 
specimens was 58.4 min, a difference of 2.7 min. 
However, statistical analysis (t-test) indicated that 
the difference of 2.7 min was not significantly dif- 
ferent from zero, since the 95 percent confidence 
limits on the difference were —1.7 to +7.1 min. 
Therefore, it is possible to say with considerable 
assurance that these experiments showed no real 
difference between plasters of the two lkghtweight 
aggregates. 

The average normalized time for the sanded 
plasters was 39.5 min, 21.6 min less than that for the 
vermiculite plasters. The same statistical test 
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mentioned above was applied to this difference and 
indicated the confidence limits were from 18.8 to 
24.4 min. Therefore the difference in fire endur- 
ances between sanded and lightweight aggregate 
plasters indicated in figure 6 is shown real. 


b. Variation of Perlite Density 


> 


Two-coat specimens of 1:2, 1:3 mix and one-coat 


specimens of 1:2% 
ach of the three densities The average times from 


two specimens of each perlite density for 1:2, 1:3 


mix were 62.2, 63.3, and 66.0 min, respectively, for 


high, intermediate, and low density; for the 1:2% 
mix (five specimens of intermediate and two each 
of high and low density), the corresponding average 
times were 57.7, 63.3, and 59.1 min. The actual 


F rise versus thickness of perlite plaster 


mix were made with perlites of | 


1.05 1.10 


specimens, with curve t 58.2 (Ax,)**. 


cated that they were of doubtful significance. 
Existing standards (ASTM C-35) place limits on the 
allowable density range of perlite plaster aggregate 
at 7.5 to 15 |lb/ft®. Although the perlites used in 
this study had densities between 7.9 and 12.9 lb/ft’, 
it is doubtful that statistically significant differences 
would have been found had the full range been ex- 
plored. 


4.4. Aging Conditions 


The specimens were aged for periods of from 28 to 
103 days and in ambients of various relative hu- 
midities. It had been intended that specimens 
should be aged at relative humidities of about 50 
percent or over 90 percent. However, due to me- 


differences were small and statistical analyses indi- | chanical difficulties, the specimens fell into three 
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Figure 6. Time to 250 °F rise versus thickness of vermiculile and sand aggregate plasters. 


The curve fitted for perlite plasters in figure 5 is shown for comparison. 





groups in terms of humidity conditioning: (1) those 
kept at over 85 percent; (2) 
90 percent but falling to about 70 percent; and (3 
those kept in the range 25 to 60 percent. 

All the specimens were kept at high relative hu- 
midity until removed from the forms at an age of 
one day. The ambient temperature of all the aging 
chambers was kept at 72+4 °F 


a. Effect of Relative Humidity of Ambient 


Due to the variation of the relative humidity of 
the ambient during the aging of some specimens, it 
was decided to compare the results for the specimens 
on the basis of the mean relative humidity during 
the final 7 days of conditioning. For the perlite 
plasters, the average time for 20 specimens at 51 
percent average relative humidity was 58.4 min; for 
8 specimens at 82 percent, 59.2 min; and for 2 
specimens at 95 percent, 71.6 min. This would 
correspond to very slight increase in time with in- 
crease in relative humidity except at very high 
values of the latter, as illustrated in figure 7. It 
would be expected that high relative humidity of 
the ambient would lead to increased moisture con- 
tent for most materials, and increased time to reach 
a temperature-rise end-point in a fire test. It 
appears evident from the results that such increase, 
if actually significant, is small for perlite plasters, 
except at high relative humidities. 

For the vermiculite plasters, the results were only 
slightly different: the average time for 5 specimens 
at an average relative humidity of 45 percent was 
61.1 min; that for 3 specimens at 78 percent was 
62.1 min; and that for 5 specimens at 95 percent 
was 83.3 min. The increase at very high humidity 
was appreciably greater than that for perlite plasters, 





© GYPSUM-PERLITE 
© GYPSUM - VERMICULITE 
© GYPSUM ~ SAND 


~ 
2) 


) 
°o 


> 
o 





TIME TO 250°F RISE,NORMALIZED, min 
o 
o 











50 60 70 80 90 
RELATIVE HUMIDITY OF AMBIENT, % 


Time to 250 °F rise versus relative humidity of 
ambient. 


FIGURE 7. 


those started at over 


as illustrated in figure 7. Since the perlite plasters 
conditioned at 95 percent relative humidity were 
all of the same mix as the vermiculite plasters, 
mix not enter into the comparison. The 
duration of aging, which is discussed later in this 
paper, may have had some influence, but does not 
fully account for the above time differences. 

The times for both perlite and vermiculite plasters 
were not affected appreciably by variation of relative 
humidity of the ambient except above 90 percent. 
There was no evidence of sensitivity of the sanded 
plaster specimens to variation in relative humidity 
of the ambient. 


does 


b. Duration 

The majority of the perlite specimens, about half 
of the vermiculite specimens, and half the sand 
specimens were aged for approximately 12 weeks; 
the remainder were aged approximately 30 days. 
Among the specimens aged at normal humidities, 
the results for those aged 30 days were comparable 
to those aged 12 weeks for each aggregate. How- 
ever, among the specimens aged at high humidities, 
the times for the perlite plasters aged 30 days were 
appreciably longer than for those aged 12 weeks; 
for the vermiculite plasters the times were also 
appreciably longer. The results that would have 
been obtained with specimens aged only a few days 
might have been significantly different, but such 
short durations were not considered of practical 
importance. 

The foregoing interactions between aging periods, 
at various humidities, and times observed in the 
fire tests might be somewhat different for much 
thicker or much thinner plaster applications, or 
for those not exposed to fairly free air circulation 
at both The 1-in. plaster thickness ap- 
proached in the specimens was believed to be a 
representative choice within the range of thicknesses 
normally applied to metal lath. Vapor barriers 
in close proximity to the back surface of plaster 
are not found in the majority of uninsulated walls, 
partit ions, or ceilings. 


faces. 


4.5. Thermal Properties 


The thermal properties that are of terest in any 
given heat-flow situation are influenced by the nature 
of the heat flow, whether transient or steady-state. 
The fire exposure employed in this study, and in 


standard fire-endurance tests, is an example of 
transient heat-flow phenomena that tends toward 
steady State as elapsed time increases. With gvpsum 
plaster, the heat flow during the early stages of a 
test is complicated by latent heat losses to evapora- 
tion of free water in the aggregate, change-of-state 
or calcination of the gypsum, and evaporation of the 
water of crystallization liberated during the calcina- 
tion. These built-in heat sinks account for the 
irregular shapes of the temperature curves in figure 4 
during the first 40 min, and for a large part of the 
value of gypsum as a fire protection material. 
Since both transient heat flow and an approxima- 
tion of steady-state heat flow occurred during the fire 
tests, both thermal conductivity, k, and thermal 
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diffusivity, a, were of interest. The steady-state 
thermal conductivity values of the individual plaster 
specimens, as determined in the NBS pentane hot- 
plate and summarized in table 3, were compared with 
the temperature-rise times from the fire tests. It 
was observed that specimens exhibiting longer times 
in the fire tests corresponded to plasters having lower 
thermal conductivity values. 

Mean values of a, and the values of k, p, and ¢ 
from which @ was computed, are given in table 3. 
Values of & and p, density, were determined from 
the 8-in. square specimens; values of c, specific heat, 
were computed from recognized formulas [3] solved 
using the temperature at which k was measured. 
The temperature-rise times were analyzed, by least- 
squares method, for a relationship with a. It was 
found that the specimens exhibiting longer times 
corresponded to lower values of diffusivity. 

It was of interest to observe the interdependence of 
density and thermal conductivity, from the data on 
both measured on the 8-in. square specimens. A 
plot of k versus p for the plasters of all mixes and 
aggregates is given in figure 8. The data fall into 
groups corresponding to the different mixes. The 
density variations among the mixes of 1:2; 1:2 
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1:3; and 1:2 for lightweight aggregates were so 
small that the data for these three similar mixes 
blended into a single large group. The two lines, one 
for the sanded plasters and the other for lightweight 
aggregate plasters, intersect very near the group of 
data for neat gypsum plaster. This was expected, 
since each line represents a range of solid dilution of 
gypsum: the one, dilution with material of lower 
density and the other, dilution with material of a 
higher density. 

A similar plot of a versus p showed a greater scatter 
than that of & versus p as percentages of the overall 
ranges of k and p. However, the data again fell into 
groupings by mix and reasonable fits were obtained 
with straight lines. The intersection of the lines 
fell moderately near the group for neat gypsum. 
The greater scatter probably reflects uncertainties in 
the chemical compositions of the mixtures and 
therefore in the computed values of specific heat, 
plus the observational and round-off errors in the 
values of & and p. 

For the sanded gypsum plasters, both k and a 
increased with increased density. However, for 
the lightweight aggregate plasters, k increased and 
a decreased with increased density. This 


agrees 
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FIGuRE 8. 


Thermal conductivity and thermal diffusivity versus density, at 117 °F. 
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with previously obtained data for concretes over a 
wide range of densities, which showed k increasing 
over the entire range but a@ decreasing at lower 
densities and then increasing at higher densities. 

The data for the concretes did not fit linear func- 
tions of density. Although the plaster data for k 
versus p were obtained under conditions that give 
considerable confidence in the linear function fit, 
the computed values of ¢ that entered into the values 
of a lead to uncertainty in the latter. Therefore, 
the straight lines drawn in the figure should be 
recognized as merely illustrating the general trend. 
The true function expressing the @ versus p relation- 
ship for these plasters may be other than linear in 
form. 


4.6. Elevated Temperature Estimates of Thermal 
Properties 


Values of density at elevated temperatures were 
calculated from the densities measured at room 
temperature and the known changes in gypsum; 
values of specific heat were calculated by means of 
the same formulas mentioned in the preceding 
section. The accuracy of results obtained by such 
procedures depends on how well the initial compo- 
sition of the plaster is known, and on the assumptions 
made as to the changes that take place under the 
influence of increased temperature. 

Determinations of steady-state thermal conduc- 
tivity at elevated temperatures were made on a 
specimen from the fire brick slab [5]. Similar 
determinations could not be made on the plaster 
specimens. However, the temperature data from 
the fire tests and the experimentally determined 
values of & for the fire brick slab were used to obtain 
estimates of the apparent thermal conductivities of 
the plasters at elevated temperatures. The changes 
of temperature with time became essentially con- 
stant in the later stages of the fire tests; the differ- 
ences in temperature across the plaster, air space, 
and fire brick slab all tended toward constant values 
as shown in figure 9. Although steady-state con- 
ditions were not reached, a quasi-steady state was 
approached. 

In the calculation of apparent conductivity, it 
was assumed that the temperature of the exposed 
plaster surface was equal to that indicated by the 
encased furnace thermocouples. This assumption 
was based on data from two fire tests of plaster 
specimens with fine wire thermocouple junctions on 
the surface. Computations made using these data 
showed that, during quasi-steady state, the average 
error introduced by this assumption was less than 
one percent of the apparent thermal conductivity 
value, and the maximum error was less than two 
percent. 

The derivation of the formula by which the appar- 
ent thermal conductivity values were computed 
began with consideration of steady-state one-dimen- 
sional heat flow through two slabs (see fig. 10), 
which may be expressed by the following heat- 


balance equation: 
T. T 
oT,_, 2 


Q=k, ne 


(1) 


where ( is heat-flow rate, & is thermal conductivity, 
OT/Oxr is temperature gradient, and subscripts p 
and s refer to plaster and fire brick slab, respectively. 
Solving (!) tor /,, we get: 


oT, /oT, s 

Or/ Oz \“ 
Since & varies with temperature, neither it nor 
temperature gradient is exactly constant through 
a slab having different temperatures at the two 
surfaces, even at steady state. Therefore, eq (2) 
implies mean values or knowledge of the values 
of k and the temperature gradients at the surfaces 
z, and x. However, steady-state gradients through 
solid slabs are ordinarily taken as the temperature 
difference between faces divided by slab thickness, 
so eq (2) becomes: 


k=k, 


? 


This equation is based on no heat loss or gain in either 
slab—fire brici. or plaster—whereas some heat was 
being absorbed in each. Therefore, the temperature 
gradients in the slabs were not linear functions of the 
z dimension. However, at some plane within each, 
the gradient must have been equal to the values given 
by the ratios in (3). It was assumed that this was 
the midplane of each, arfd that the heat flow at Z,, the 
midplane of the plaster, was equal to that at 7,, the 
midplane of the fire brick slab, plus the rate at which 
heat was absorbed in heating the shaded portions in 
figure 10. From these assumptions, the following 
was derived: 


oT, /oT, oT,\faz, _ dT, dz, _ dT. 
S| Be +(1 Ox J) 2 Por ay Tg Pssae P 


where all the temperature gradients with respect to 
thickness and with respect to time are at Z, or Z,. 
Substituting approximations for the partial deriva- 
tives gives: 


re) } (4) 


Chee X1—Zo — 23-2, =— 
where ¢ is time, —5 1-2) — Ls 


other symbols are as previously used. 
double prime are used with k, in eqs (3) and (4), as 


ZY and the 


The prime and 
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well as in table 4 and the subsequent discussion, to 
emphasize that the formulas, when used with the 
data from the fire tests, gave values of “apparent”’ 
thermal conductivity, obtained under quasi-steady- 
state conditions. 

It is of interest to note that the portion of eq (4) 
to the left of the first plus sign is identical to eq (3); 
the remainder of eq (4) being a correction for devia- 
tion from steady state. Strictly, the major brackets 
in (4) should include a correction factor for heating 
the air between the slabs. In fact, sample calcula- 
tions indicated the correction for air would be of the 
order of one-fortieth of one percent of k,’. This would 
not show in the numerical values in table 4 unless 
they were given to two additional figures. Also, a 
correction should be made for edge effects. However, 
examination of the temperature data for the bottom 
face of the fire brick slab with that for the bottom 
face of the asbestos board (see fig. 2) showed small 
temperature differences and, therefore, small edge 
effects due to lateral heat flow. 

Values of k were computed from data for neat 
plaster, gypsum-perlite and gypsum-vermiculite plas- 
ters at 1:24 mix and gypsum-sand plaster at 1:2 mix. 
The values were plotted against mean _ plaster 
temperature and curves faired through them, 
except for the sanded plaster. The latter indicated 
k’; of the order 3 to 5 but were too widely scattered 
to be a satisfactory basis for a curve. This result is 
not surprising in view of the comparatively short 


TABLE 4. Apparent and computed values of thermal properties 


at elevated temperatures 
Apparent conductivity ’’ determined from data; specific heat computed from 


formulas; density based on measured density at room temperature; apparent 
diffusivity a’ computed by a’ =k’’/pe 


Temperature 
| 


1,300 °F 


1,400 °F 


Apparent conductivity, k’’, Btu/hr ft? (°F/in.) 


Neat gypsum -. 2 2.00 2.35 
Gypsum-perlite, 1:2'9-. 1. 50 i, 1.60 
Gypsum-vermiculite, 1:24 nieltg 1. 65 N 1.70 
Gypsum-sand, 1:2..... oe (¢ 


Btu/lb °F 


Specific heat, c, 


Neat gypsum 
Gypsum-perlite, 1:24 
Gypsum-vermiculite, 
Gypsum-sand, 1:2 


283 
281 
281 
299 


0. 302 
. 298 
. 298 
317 


ih/ft8 


Neat gypsum...‘ 
Gypsum-perlite, 1:24 
Gypsum-vermiculite, 1:24 
Gypsum-sand, 1:2 


75.0 
42.2 
39.9 
102.0 


ipparent diffusivity, a’, ft?/hr 


Neat gypsum. 
Gypsum-perlite, 1:24 ‘ 0106 
Gypsum-vermiculite, 1:24 " .0122 
Gypsum-sand, 1:2 : ot ee (2) 


0. 0080 0. 0083 
0105 


.0120 


0. 0086 
.0107 
.0118 


0.0090 
.O112 
-O117 


® Data gave k”’ values in range 3 to 5 but widely scattered; therefore they were 
not tabulated, nor used to compute a’, 


| 
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durations of the tests with gypsum-sand plasters, in 
which even quasi-steady state was not approached. 

The curves obtained are shown in figure 11. 
The data for gypsum-perlite and gypsum-vermic- 
ulite were so intimately associated that the latter 
were omitted in order to avoid unduly reducing the 
clarity of the figure. The scatter of the omitted 
gypsum-vermiculite points was similar to that of 
those for gypsum-perlite. Table 4 gives values of 
k’; taken from the curves, along with computed 
values of c, p, and a’. 

As a check on the values obtained by using eq (4), 
two gypsum-perlite plaster specimens of 1:2'5 mix 





MEAN PLASTER TEMPERATURE ,°F 


FIGURE 1] 


were tested with the same instrumentation emploved 
throughout the study, but with the furnace pro- 
gramed for temperatures. Data were 
obtained at four temperatures, under steady-state 
conditions, and are shown in figure 11. The curve 
through these is significantly above, but of the same 
shape as that for the data from eq (4). Two factors 
undoubtedly contribute to this discrepancy between 
ky andk,. First,in eqs (3) and (4) there is the product 
i ‘ d d 
%—z, 1,—T, 
sion for plasters are so small that values of the z’s 
are essentially constant but, as shown in figure 9, 
T,— T; (curve B) approaches steady state from lower 
values and 7)— 7, (curve C) approaches from higher 


constant 


The coefficients of thermal expan- 


values. Therefore, the ratio 


~ = would be ex- 


pected to remain below corresponding steady state 
values throughout fire tests, 


of &.. 


resulting in low values 
Secondly, the values of 


w—(To+T, 
ZAt 
and of 
(T,+T,)... T, 
2At 


were used to approximate values of 07'/O¢ at the mid- 
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planes of the plaster and firebrick slabs. Although it 
was assumed that the midplane gradients were equal 
to the overall gradients, the midplane temperatures 
would be less than the average of the surtace tem- 
peratures. The difference between actual and ap- 
proximated midplane temperatures would diminish as 
time increased, and steady state was more nearly ap- 
proached, so the rate of increase of the midplane 
temperature would be greater than the mean rate of 
increase of the two surface temperatures. 


Therefore, 
the correction terms in eq (4 


for both the plaster 
and the slab should be expected to be low. However, 
the agreement between the two curves for gypsum- 
perlite plasters in figure 11, both in shape and order 
of magnitude, were such as to give some confidence 
in the use of eq (4). 

Several important points must be stressed in con- 
nection with the estimates. First, they are at best 
only estimates. Second, they include the effects of 
cracks and all the other changes in the plaster 
resulting from fire exposure. It is because of the 
tendency of gypsum plaster to crack when heated to 
elevated temperatures that thermal conductivity at 
high temperatures could not be measured by labora- 
tory equipment and techniques used to obtain the 
values in table 3. However, the fact that the values 
in table 4 include the effect of cracks is considered 
an advantage because the values are, therefore, 
more closely related to the behavior of 


gypsum 
plaster as a fire-protection material. 





A third important point is in regard to extrapola- 
tion of the curves in figure 11. No extrapolation to 
lower temperatures should be attempted to obtain 
values for fire-exposure or other non-steady-state 
computations since the effects of moisture evapora- 
tion and migration will be important at lower mean 
plaster temperatures. However, limited extrapola- 
tion to lower temperatures might be acceptable to 
obtain values for steady-state computations. Con- 
sideration of plaster thickness is a fourth point. 
For thicknesses much less than those involved in this 
study, e.g., 4 in., the effects of moisture would be 
eliminated more rapidly, corresponding to lower 
plaster mean temperatures. In this: case, some 
relaxation of the above prohibition against extrap- 
olation to lower temperatures might be permissible. 
Conversely, at thicknesses of 14 in. or greater, the 
effects of moisture might not be eliminated until 
mean plaster temperatures exceed 1,200 °F. 

It is worth noting that, for plasters of the same 
mixes, the values of k, given in table 3, for 117 °F, 
are higher than those given for 1,200 °F in table 4, 
although the latter and figure 11 both indicate higher 
k, at increased temperatures. This points up the 
significant effect of the loss of chemically combined 
and other moisture on k,. This fact is very desirable 
from the fire-protection viewpoint, but complicates 
the problem of predicting fire endurance, or any heat- 
flow solution involving the calcination of gypsum, 
by computational or analytical means. 


5. Summary 


This study was designed to examine the effects of 
mix, aggregate, and conditioning on the fire endur- 
ances, in terms of a limiting temperature rise, of 
gypsum plasters. It was not intended to examine 
the effects of large variation of thickness, nor was it 
felt that such was needed because of the published 
data from many standard fire tests and conclusions 
based thereon. 

Under the conditions of specimen size and _ test 
conditions employed, the study showed that the 
mix ratio of gvpsum to aggregate has little, if any, 
influence on the time to a 250 °F temperature rise, for 
a single aggregate within the range of mixes ordi- 
narily employed. The perlite and vermiculite 
aggregate plasters exhibited significantly longer 
temperature-rise times than did the sanded gypsum 
plasters. This was attributed in part to the differ- 
ences between the aggregates and in part to the fact 
that the conventional mixes for sanded plasters 
correspond to quite different ranges of gvpsulm-to- 
aggregate ratios on a weight basis. Variation of the 
density of perlite aggregate over a range of about 8 
to 13 Ib/ft*had no. significant influence on the 
limiting temperature-rise time. 

The effect of aging or conditioning was pronounced 
for the combination of very high relative humidity 


atmospheres and relatively short periods of aging. 
For perlite and vermiculite aggregate plaster speci- 
mens conditioned 30 days in atmospheres at relative 
humidities near 95 percent and temperatures near 70 
to 75 °F, significantly longer times were observed 
than for equivalent specimens aged about 12 weeks 
under similar temperature and humidity conditions. 
Duration of aging beyond 30 days in atmospheres 
below 85 percent relative humidity had no influence. 

The data from the fire tests of neat gypsum, 
gypsum-perlite, and gypsum-vermiculite plasters 
served as bases for the calculation of estimates of the 
apparent thermal conductivities and thermal dif- 
fusivities of such plasters at temperatures in the 
range 1,200 to 1,500 °F. Values of the same thermal 
properties were obtained for steady-state heat flow 
at about 120 °F mean plaster temperature. No 
attempt was made to calculate values from the fire- 
test data for the range 200 to 1,000 °F because of the 
latent heat effects associated with calcination of 
gypsum early in fire-test exposures. However, the 
values for higher temperatures, although only 
estimates, should be useful in the prediction of the 
fire endurances of plastered structures by analog 
devices or computational procedures. 

Examination of the data in terms of plaster thick- 
ness, over the small range of variation observed, 
showed agreement with the earlier findings from 
large specimen tests, so far as a limiting temperature 
rise is concerned. It is possible that large departures 
from the nominal 1-in. plaster thickness may affect 
the application of the high-temperature values of 
thermal properties, due to the effects resulting from 
the water in the hydrated gypsum. 


The author expresses his appreciation to Garnett 
Robinson, who prepared the plaster bases and did 
the plastering. 
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Calorimetric calibration of an ionization chamber for deter- 
mination of X-ray total beam energy, J. S. Pruitt and S. R. 
Domen, J. Research NBS 66A (Phys. and Chem.) No. 6, 
(Sept.—Oct. 1962). 70 cents. 

The design and use of a high sensitivity lead calorimeter is 
described in this paper. It can measure the energy trans- 
ported by an X-ray beam with an intensity as small as 2 
microwatts/cm?, and has been used to calibrate a duraluminum 
chamber developed in this laboratory with a maximum error 
of +2 percent for all bremsstrahlung energies between 20 and 
170 Mev. The results obtained are in excellent agreement 
with independent calibrations obtained with a scintillation 
spectrometer. The calibration of the copper chamber de- 
signed by Edwards and Kerst has also been determined over 
this energy range, to within +2.6 percent, and agrees well 
with the high energy calibrations reported by these authors. 


Effect of exposure site on weather resistance of porcelain 
enamels exposed for three years, D. G. Moore and A. Potter, 
NBS Mono. 44 (Apr. 10, 1962) 15 cents. 

An exposure test of porcelain-enameled steel and aluminum 
specimens is being conducted jointly by the Porcelain Enamel 
Institute and the National Bureau of Standards. The expo- 
sure sites are Dallas, Tex.; Los Angeles, Calif.; New Orleans, 
La.; Pittsburgh, Pa.; Washington, D.C.; and two sites at 
Kure Beach, N.C 

The present report gives the results of the third-year inspec- 
tion. Changes in gloss and color were taken as criteria of 
weathering. Based on averages for all enamels, the 80-ft site 
at Kure Beach caused the greatest change in gloss and color, 
while the conditions at Dallas, Los Angeles, and New Orleans 
caused the least change. The conditions at Washington, 
Pittsburgh, and the second Kure Beach site (800 ft from the 
ocean) were of intermediate severity. Comparison of the 
gloss and color changes with data obtained by the National 
Air Sampling Network indicated that air pollution by acidic 
contaminants was a factor in site severity. 

A direct relation existed between acid resistance and weather 
resistance. This relation was apparent, however, only when 
averages were considered. There were individual exceptions 
within groups of enamels of the same general type. In 
addition, enamels of different types, such as aluminum and 
steel enamels, having the same acid resistance (citric acid spot 
test) did not necessarily show the same weather resistance. 
Further, some red and yellow enamels with good acid resist- 
ance showed poor color stability. It was found, however, that 
this poor stability could be predicted by a specially developed 
cupric sulfate test. 

As a group, the regular glossy enamels for steel showed the 
best weather resistance among the various types tested. 


Fire tests of precast cellular concrete floors and roofs, J. V. 
15 


Ryan and E. W. 


cents. 


Bender, NBS Mono. 45 (Apr. 12, 1962) 


The results of an investigation of lightweight, precast cellular 
concrete planks are given. Fire tests were made of two floor 
and five roof specimens made up of these planks. Variables 
included density of the cellular concrete, thickness and span of 
the planks, reinforcement, and cover for the latter. A steel 
beam encased in blocks of cellular concrete was included in 
one floor specimen. The flexural strengths of 14 individual 
planks determined. The investigation showed fire 
endurances up to 2 hr for 6-in. thick slabs tested under load 
and up to 4 hr for other slabs not loaded. Estimates were 
made of the probable results to be expected for slabs of thick- 
nesses other than those actually tested, 
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Structure shielding against fallout radiation from nuclear 
weapons, L. V. Spencer, NBS Mono. 42 (June 1, 1962) 75 cents. 
The theory of structure shielding from fallout gamma radiation 
is developed to the point of applications to elementary struc- 
ture types. Examples discussed in the text include the 
density interface, foxhole, shielded foxhole or basement, light 
superstructure, vertical wall, blockhouse, vents, compart- 
mentalization effects, and mazes. A large number of engineer- 
ing charts and graphs are presented for engineering calcula- 
tions, including many obtained from angular distributions of 
the exposure dose. Results are given for a fission spectrum, 
and for Co® and Cs!*7 sources. This information has been 
obtained almost completely by machine calculations utilizing 
basic cross section data. A number of sources of experi- 
mental data are mentioned, but detailed comparisons with 
experiment are not included. 


Basic magnetic quantities and the measurement of the 
magnetic properties of materials, R. L. Sanford and I. L. 
Cooter, NBS Mono. 47 (May 21, 1962) 30 cents. 

This paper gives general information regarding the two basic 
quantities, magnetic induction, B, and magnetizing force, H, 
and also the magnetic constant T,, (often designated by the 
symbols uw, and »w,). Information is also given regarding the 
magnetic properties of various materials and methods and 
apparatus commonly used in the Magnetic Measurements 
Section for measuring these properties by means of reversed 
direct current or alternating currents of low frequency. 
Magnetic measurements peculiar to high frequencies are not 
discussed. In view of the gradual adoption of the rationalized 
mksa system of units, this system is included as well as the 
classical cgs electromagnetic system. 


Determination of total X-ray beam energy with a calibrated 
ionization chamber, J. S. Pruitt and S. R. Domen, NBS 
Mono. 48 (June 5, 1962) 20 cents. 

This report describes the use of an air-filled aluminum-alloy 
ionization chamber to determine the energy transported by 
a bremsstrahlung beam with maximum photon energy in the 
range 6 to 170 Mev. The experimental calibrations of this 
chamber over this energy range are given, as well as the results 
of calibration experiments made with a 250-kv constant- 
potential X-ray tube and with Cs"? and Co” gamma-rays. 
Information is presented about the change in calibration when 
the chamber is used with different experimental conditions, 
and when either its dimensions or its alloy composition are 
changed slightly. This report can be used to construct a 
replica chamber and to determine its absolute calibration 
between 6 and 170 Mev under a variety of experimental 
conditions. 


Some problems of fatigue of bolts and bolted joints in aircraft 
applications, L. Mordfin, NBS Tech. Note 136 (PB161637) 
(1962) $1.25. 

The profuse variety of aircraft bolts which is available has 
made the evaluation and specification of bolts for engine and 
structural use extremely complex, particularly insofar as 
fatigue and hot fatigue environments are concerned. The 
state of knowledge of fatigue of bolts and bolted joints is 
surveyed and critically appraised in terms of aeronautical 
practices. Using this material as a basis, recommendations 
are made regarding the evaluation and specification of air- 
craft bolts for fatigue situations and regarding the growing 
problem of errors in fastener replacement. 


Cryogenic temperature measurement with platinum resistance 
thermometers—Is fixed-point calibration adequate? R. J. 
Corruccini, NBS Tech. Note 147 (PB161648) (1962) 50 cents. 
An analysis of extensive calibration data for strain-free 
‘“capsule’’-type platinum resistance thermometers indicates 
that the above question can be answered affirmatively for 
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Characteristics of point-to-point tropospheric propagation 
and siting considerations, R. S. Kirby, P. L. Rice, and L. J. 
Malonev, NBS Tech. Note PB161596) (1961) $2.50. 
This discussion is primarily intended to provide an under- 
standing of the fundamental characteristics of tropospheric 
propagation with particular application to point-to-point 
telecommunications The concept of service probability is 
introduced and its application to the planning of tropospheric 
circuits. This concept provides an objective means for taking 
into account the variables and uncertainties connected with 
a planned circuit so that a reasonable balance can be made 
the cost of installation and operation as compared 
with the probability of suc Principles of siting based on 
taking advantage of favorable technical characteristics of the 
as well as considerations such as radiation 
hazards discussed. methods of making pre- 
liminary estimates of performance adequate for field use are 
presented. 


OMNIFORM I: A general purpose machine program for the 
calculation of tables of functions given explicitly in terms of 
one variable, J. Hilsenrath and G. M. Galler, NBS Tech 
Vote 125 (PB161626) (1962) $1.00. 
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The thermodynamic properties of nitrogen 
between 0.1 and 200 atmospheres, T. R. 
Tech. Note 129 (PB161630) (1962 
The internal energy, enthalpy 
of molecular nitrogen are 


from 64 to 300 °K 
Strobridge, NBS 
$2.95. 
entropy, and specific volume 
derived and tabulated as functions 
of temperature and pressure. In addition to a mathematical 
model for the pressure-volume-temperature surface, accurate 
functions are given for the representation of the vapor 
pressure, density of saturated liquid, specific heat of saturated 
liquid, and the specific heat at zero pressure. 
Tabular values in British units over the same pressure and 
temperature range available as Supplement A of this 


Technical Note. 


are 


Detailed techniques for preparing and using hard gallium 
alloys, G. G. Harman, NBS Tech. Note 140 PB161641 
1962) 75 cents. 
This report presents expansion and clarification of tech- 
niques for preparing and using dental-amalgam-type gallium 
alloys for industrial and scientific uses that were previously 
published in the Review of Scientific Instruments In addition, 
new material, such uardness of the alloys and discussions 
of the bonding mechanism are presented. A different class 
of bonding alloys based on the same principles but not in- 
corporating gallium are 
details. 


‘ } 
as I 


described, along with application 


Dielectric constant of liquid parahydrogen, R. J. 
NBS Tech. Note 144 (PB161645) (1962) 50 cents. 
It is shown that the available data on the dielectric constant 
of hydrogen conform to the Clausius-Mossotti equation within 
the probable experimental errors. The published data cover 
temperatures from the triple point (~14 °K) to well above 
room temperature and a thousandfold range of densities 
Using an average value of the specific polarization, tables of 
the dielectric constant of the liquid have been computed for 
temperatures from the triple point (13.803 °K) to 32 °K 
(58 °R) and pressures from saturation to 340 atmospheres. 


Corruccini, 


| 


A wire exploder for generating cylindrical shock waves in a 
controlled atmosphere, D. L. Jones and K. B. Earnshaw, 
NBS Tech. Note 148 (PB161649) (1962) 50 cents. 

A design for a rugged exploding wire device is given. This 
device permits the study of strong cylindrical shock waves in 
controlled atmospheres using optical and microwave tech- 
niques Adequate detail and pictures given to allow 
construction of the device. 


are 


Displacement and strain-energy distribution in a longitudi- 
nally vibrating cylindrical rod with a viscoelastic coating, 
r. Hertelendy, J. Appl V ech. Trans. AS VE 29, Series E, 
No. 1, 47-52 (Mar. 1962 

A numerical solution by R. M. Davies of the Pochhammer 
frequency equation is used to determine the displacement and 
strain energy distribution the cross-section of an 
infinite elastic circular cylindrical rod for a number of wave- 
lengths of the first, second, and third modes of symmetric 
iongitudinal vibration With these results the effect of a 
thin uniform layer of damping material on the surface is 
investigated, and it is shown that while first mode high 
frequency waves are damped out very rapidly, a large group 
of second and third mode waves experience virtually 
damping. 


across 


ho 


Properties of silico-phosphate cements, J. N. 
G. C. Paffenbarger, Dental 
1962). 

Seven brands of silico-phosphate cements, 
their manufacturers 


Anderson and 
’rodgress _ No. 2, 72-75 (Jan. 
recommended by 
as dental cementing media, were tested 
by standardized procedures (A.D.A. Specifications No. 8 
and 9 Their use in countries is accepted, but in 
others their value is questioned About half of the cements 
tested did not set 10 minutes All had an excessive 
film thickness if mixed thick enough to have a proper setting 
time. The solubilities of the silico-phosphate cements in 
water is higher than those of the zine phosphate cements. 
However, the silico-phosphate cements are not more soluble 
in the mouth since the solubility test is valid within a type of 
cement and cannot be used for intercomparison of dental 
cement based on phosphoric acid. The opacity of some of 
the cements is in the of the opacities of enamel and 
dentin The silico-phosphate cements therefore, the 
most suitable cementing media for translucent restorations. 


some 


within 


range 


are, 


Surface effect on bond strength of steel beams embedded in 
concrete, J. O. Bryson and R. G. Mathey, J. 
Inst. 59, No. 3, 397-406 (Mar. 1962 

Wide flange structural steel beams with different surface 
conditions were embedded in concrete and subjected to push- 
out tests to determine the effect of the surface condition on 
the bond between the concrete and steel. The surfaces of 
the embedded steel beams were either freshly sandblasted, 
sandblasted and allowed to rust, or left with normal rust and 
mill seale. 

The steel beams with a sandblasted surface, and those with a 
sandblasted surface which was allowed to rust, developed 
considerably higher ultimate bond stresses than beams with 
normal rust and mill seale. However, at a free end slip of 
0.001 in. there was no significant difference in the bord stress 
for all three types of surface conditions. 


Am. Concrete 


Preparation of and electroplating on uranium, D. 
Plating 49, No. 4, 363-367 (Apr. 1962). 

The black film that forms between uranium and 
electrodeposited on it results from the solutions 
aqueous cleaning and plating procedures. It was eliminated 
by using a completely anhydrous plating solution. <A 
porosity test, using hydrogen, was developed for testing the 
porosity of plated coatings on uranium. 


E. Couch, 


coatings 
used in 


Investigation of the spectrophotometric method of measuring 
the ferric ion yield in the ferrous sulfate dosimeter, K. Scharf 
and R. M. Lee, Radiation Research 16, No. 2, 115-124 (Feb. 
1962 

The ferric ion yield produced by ionizing radiation in a 
ferrous sulfate dosimeter is usually determined by measuring 
the absorbance (optical density) of the irradiated solution at 
a wavelength of 304 mu. The lowest absorbance which can 
be measured with sufficient accuracy sets the lower limit of 
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the radiation dose which can be determined by this method. 
!t was found that by measuring the absorbance at an absorp- 
tion peak of 224 muy, instead of at 304 my, the sensitivity of 
this method can be approximately doubled. This is due to 
the higher molar extinction coefficient of ferric ions at 224 
my which is 4565 liter mole~'em~', as compared with a value 
of 2196 liter mole~!em=! at 304 my, both measured at 25 °C. 
In the temperature range from 20 °C to 30 °C, the value of 
the extinction coefficient at 25 °C changes by .13% per degree 
at 224 my, compared with a change of .69% at 304 mu. 
The dependence of the extinction coefficient of ferrie ions on 
the normality of the sulfuric acid (solvent) is less at 224 my 
than at 304 my. Values of the ferric ion yield produced by 
Co® gamma rays, determined by measuring the absorbance 
at 304 my and 224 my, agreed within experimental errors. 
Tensile strength and modulus of elasticity of tooth structure 
and several restorative materials, KR. L. Bowen and M. 8. 
Rodriguez, J. Am. Dental Assoc 64, Vo. 3, 278-287 (Mar. 
1962). 

Values for strength and modulus of elasticity were obtained in 
tension for various dental materials 
stored in distilled water at 99 °F. for seven days. The speci- 
mens were loaded at 0.02 inches per minute head speed. 
Dentin and enamel specimens were prepared by cutting with 
They 
were formed into rectangular cross-section rods with a nar- 
rowed middle portion and notched ends, which were imbedded 
in direct filling resin. The materials had the following average 
tensile strength in pounds per square inch: enamel, 1,500; 
dentin, 7,500; silicate cements, 700; zine phosphate cements, 
500; direct filling resins, 4,200 and an experimental direct 
filling material containing about 70 percent treated fused 
silica, about 4,800 depending on surface treatment of the 
silica powder. The materials had the following average 
moduli of elasticity when measured in tension (times 108 
pounds per square inch): dentin, 2.8; a silicate cement, 3.1; 
direct filling resins, 0.27 and an experimental direct anterior 
filling containing 70 percent vinyl silane 
treated fused silica powder, 1.6. 


The specimens were 


rotating diamond tools using continuous water spray. 


material about 


Internal magnetic fields in 
L. H Jennett and R. L. Streever, Jr., J 
33, Vo. 3, 1093 1094 | Var. 1962). 
The nuclear magnetic resonances of both the Ni®! nuclei and 
the Co” nuclei have been observed in a series of ferromagnetic 
nickel alloy powders containing up to 2.4 atomic percent co- 
balt. The room temperature value of the magnetic field at 
the site of the Co” nucleus in any of these alloys is found to be 
about L111 koe. Thisis considerably smaller in magnitude than 
the magnetic field at the nickel site in either the alloy or in 
pure nickel, which, assuming a Ni®! nuclear moment of 0.30 
nuclear is about 170 koe, or at a cobalt site in 
pure cobalt, which is about 211 koe. The line widths of both 
the nickel and cobalt resonances increase with concentration 
with no detectable shift in frequency. The temperature 
dependence of the Co™® resonance frequency in the nickel 
0.6 percent cobalt alloy has been measured, and found to 
be different from that of Ni®' in pure nickel and in the alloy. 


An aluminum magnet cooled with liquid hydrogen, J. R. 
Purcell, Proc. Intern. Conf. on High Magnetic Fields, Mass. 
Inst. of Tech , Cambridge, Vass... Vov. 1 ,. 1961, pp. 166-169 
( Vass Inst Tech Press, Cambridge, Vass., and John W ile y 
& Sons, Inc Vew York, N.Y., 1962). 

A high field, steady state magnet has been developed that 
utilizes liquid-hydrogen forced-convection cooling. The 
solenoid is constructed of stacked reels wound with high 
purity aluminum foil, and employs unique radial ducts for 
coolant flow This design allows for high stacking factor, 
low current requirements and simplicity of construction, 
along with a minimum of total power needs. The construc- 
tion of the magnet and auxiliary equipment is described, and 
the results of preliminary testing are given. 

Static dielectric constant of rutile (TiO,), 1.6—-1060°K, R. A. 
Parker, Phys. Rev. 124, No. 6, 1719-1722 (Dec. 1961). 

The static dielectric constant e« of rutile has been measured 
as a function of temperature from 1.6 to 1060 °K. Rutile 
does not appear to be ferroelectric or antiferroelectric in this 
temperature range. At low 


nickel-rich nickel-cobalt alloys, 
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magnetons, 


temperatures, « approaches a | 
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limiting value of 257 in the c-direction and 111 in the a-diree- 
tion. eg90xK=170 and 86, respectively, and  e:o0K=97 
and 58. The ionic polarizability of titanium is calculated 
using the previously derived Lorentz correction in the rutile 
structure. The polarizability at all temperatures is very 
close to the critical value for a ferroelectric polarization ca- 
tastrophe. The temperature variation of the ionie polariz- 
ability is explained qualitatively on the basis of other 
properties of the crystal. 


Electronic conduction in rutile (TiO,.), H. P. R. Frederikse, 
W. R. Hosler, and J. S. Becker, Proc. Intern. Conf. Semi- 
conductor Phys., Prague, Czechoslovakia, pp. 868-871 (1960). 
Measurements of conductivity, Hall effect and thermoelectric 
power performed on slightly reduced crystals of rutile can be 
explained with a model involving two conduction mechanisms: 
one of free electrons in a narrow 3d-band of titanium (above 
5 °K , and the other of self-trapped electrons or polarons in 
a polaron band (below 5 °K). The ionization energy of the 
polaron increases from ~0.01 ev to 0.07 with increasing 
temperature. 


Modulated 


photoelectric measurement of vibration, V. A. 
Schmidt, 8. 


Edelman, E. R. Smith, and E. T. Pierce, J. 
Acoust. Soc. Am. 34, No. 4, 455-458 (Apr. 1962). 

This paper describes an improvement in the method of 
setting vibration amplitudes for calibrating vibration pickups 
by the disappearance of an interferometric fringe pattern. 
In the usual method, one plate of an interferometer is station- 
ary, the other vibrates with the pickup, and the fringe pattern 
disappears at zeros of the Bessel function 


Jo (4a &) 


p 


where £ is the amplitude of vibration and \ is the wave length 
of light used In the improved method, the fringes are 
observed by a photomultiplier, the previously stationary 
plate is vibrated at a modulating frequency much lower 
than the calibrating frequency applied to the pickup, and the 
signal from the photomultiplier, filtered at the modulating 


J (4a €) 


frequency, has minima at the zeros of ° Observation 


of the nulls on a meter allows faster calibration with greater 
precision and less observer fatigue. 

A high mode tunable cavity 
A. J. Estin, Rev. 
A high-Q 


or microwave-gas interactions, 
Sct. Instr. 33, 369-371 (Mar. 1962). 

resonant microwave cavity, 
TEoin(n=10, 11) mode, is described. A fused quartz liner, 
with suitable inlet and pumping lines, is used to contain 
gases, so that interactions between the molecular transitions 
and the electromagnetic radiation can be generated. 
Techniques are discussed for suppression of unwanted modes 
of oscillation, for dielectric tuning of the cavity, and for 
elimination of radiation through a large port. 


operating in the 


Study of the effect of large aperture on the performance of an 
Ebert spectrometer, L. R. Megill and L. Doppleman, J. Opt. 
Soc. Am. 52, No. 3, 258-261 (Mar. 1962). 

A description of the effect of large aperture upon the resolution 
limit of an Ebert spectrometer is discussed. The calculations 
were made by ray tracing on a high-speed computer. Varia- 
tions of the aberration pattern at the exit slit are noted for 
several parameter changes which may be useful in design 
considerations. 


Prediction of pressure drop in two-phase single-component 


fluid flow, M. R. Hatch and R. B. Jacobs, A. J. Ch. E. 
Journal 8, No. 1, 18-25 (Mar. 1962). 


Data on pressure drop in two-phase, single-component fluid 
flow, both with and without heat transfer, are presented in 
terms of the Lockhart and Martinelli correlation parameters. 
The fluids used were trichloromonofluoromethane and 
hydrogen. 

The results are compared with the correlation curve recom- 
mended by Martinelli and Nelson and give frictional pressure 
drops that are about 40% lower than the curve. 
for the deviations are discussed in terms of 


The reasons 
the effects of 





friction factors, the 
accuracy of data 
procedures. 

It is concluded that the Martinelli and Nelson correlation and 
a simple momentum pressure drop computation can be 
superposed to predict roughly the total pressure drop in 
tubes containing steady state, two-phase, single-component 
fluid flow with appreciable vaporization. 


existence of 
and 


metastable 


equilibrium, 
instrumentation, and 


calculation 


Potential heat: A method - measuring w; heat release of 
materials in building fires, J. J. Loftus, D. Gross, and A. + 


tobertson, Am. Soc. Testing Materials Proc. 61, 1336-134 
(June 1961 


Modifications have been made to a method, currently a 
standard of the French government, for assessment of 
the combustible characteristics of building materials. This 
method makes use of standard calorimetric techniques in 
which the burning of small quantities of combustible in an 
otherwise inert material is assured by use of a combustion 
promoter which is added prior to test. By performing cal- 
orimetric Measurements both before and after exposure to a 
“standardized fire’’ (two hours in air at 750 °C) the differ- 
ence may be considered as the potential heat. The test 
method has been used for measurement of the potential heat 
of a variety of building materials, including metals and ma- 
terials of low combustibility, with results in general agree- 
ment with fire experience discussion is included on the 
precision of the test method 


Comparison of standard 
rubber hardness, R. D 
Res. & Std 


and microtests for 
Stiehler and F. L. 
ASTM Bu 2, No. 5, 388-391 


international 
loth, Vater. 
Vay 1962). 


The two standard testers for 
tional Rubber Hardness 
D 1415 were compared 
Twenty-two rubber 
about 30 to 90 
(IRHD 


with the 


the determination of 
ASTM 
electronic 
ranging in 


Interna- 
Designation 
microtester. 
hardness from 
tubber Hardness Degrees 
comparison. Results obtained 
agreed within about one IRHD. The 


described in 
with an 
compounds 
International 
used in the 


testers 


were 
three 


standard tester had slightly better precision, but a visual null 


indicator for the microtester 
The microtester was more 
e.g., hardening caused by 
in the 


should 


increase its precision. 
responsive 


to surface conditions 
oxidation) and to non-uniformities 
rubber at or near the surface; this makes the instru- 
ment particularly useful in research investigations as well as 
for the measurement of hardness on small rubber 
This comparison indicated that a microtest would be 
able addition to ASTM D 1415. 


parts. 
a valu- 
Designation: 


Emission flame photometry, M. 


Margoshes, Anal. 
No. 5, 221R-224R (Apr. 1962). 


Chem. 34, 


During the period covered in this review, the beginning of 1960 
through February 1962, there has been some change in the 
emphasis of research in flame photometry. The interest in 
high temperature flames, such as cyanogen-oxygen, has 
lessened, while the application of organic solvents as a means 
of improving sensitivity has become more general. Methods 
have been developed for elimination of some anion inter- 
ferences by means of releasing agents which prevent the 
formation of refractory compounds of the analytes. A num- 
ber of papers on chemiluminescence have been presented at 
meetings, although little of this information has yet appeared 
In print. 
The determination of sodium and potassium is still by far the 
most common use made of flame photometry. Because this 
area of application is so well established, very few of the 
papers reporting on sodium and potassium determinations are 
cited in this review. 


Thermal explosives: 
propellants, D. 
(Intern.) on 


Wilkins 


‘s: adiabatic self-heating of explosive and 
Gross and A. B. Amster, Fighth Symp. 


728-734 (The Williams and 
Co., Baltimore, Md., 1962 


Combustion, pp. 
An explosive may be defined as a substance in a metastable 
state of internal equilibrium which, being subjected to a com- 
paratively small external activation, turns rapidly into a hot 


gas producing a very high local pressure. The sensitivity of 
an explosive is defined in terms of the magnitude of this 
activation, i.e., the greater the required activation, the less 
sensitive the mate rial is said to be. Many different methods 
have been used to measure explosive sensitivity often with 
conflicting results. Any attempt to explain these anomalies 
leads to a requirement for knowledge of the thermal, thermo- 
chemical and kinetic properties of the systems involved. For 
a number of explosives and propellants, this paper reports some 
measurements of these properties and suggests a correlation 
with one type of sensitivity, namely sensitivity to impact. 


Tungsten helical-spring microbalance, BS: Ra 
Vacuum Microbalance Techniques 2, 47 b? 


Inc., New York, N.Y., 1962). 


This paper describes a helical spring microbalance 
mil tungsten wire and enclosed in a Pyrex glass housing. It 
was used in connection with a problem involving measurement 
of rates of thermal degradation of polymers in a vacuum by 
the loss-of-weight method. The balance has a sensitivity of 
about 550 w/mg and is so constructed as to avoid errors that 
might be caused by the thermal expansion of its parts or of the 
supports holding the balance and the cathetometer used in 
connection with it. 


Madorsky’ 


(Plenum Press, 


, made of 30 


Electrical properties of nonstoichiometric semiconductors, 
J. H. Becker and H. P. R. Frederikse, J. Appl. Phys. Suppl. 
33, No. 1962). 


The electrical conductivity ¢, thermoelectric power Q, and 
Hall coefficient R are examined as a function of the ratio of 
hole-to-electron concentrations p/n for a non-degenerate semi- 
conductor at constant temperature. From these relations 
the fundamental parameters of the material (forbidden hand 
gap, mobilities and effective masses) can be derived. This 
approach is particularly applicable to materials whose stoichi- 
ometry varies as a function of temperature and vapor pressure 
of the constituents, P. For any model of this equilibrium 
decomposition, it is easy to transform the calculations in 
terms of p/n into results as a function of P. As p/n increases, 
o passes through a minimum, while Q and F# traverse minimum 
(negative), zero, and maximum (positive) values. These 
extrema are of special interest. In the simple case of one 
kind of imperfection, ¢, Q, and R become independent of P 
in a certain pressure range (i.e., when the intrinsic condition 
n=p has reached It is then possible to derive the 
ratio of mobilities u,/u, and the ratio of the average effective 
masses m*,/m*, from o(P) and Q(P) only. Hence, if wa or 
m*, are known (i.e., from measurements at lower tempera- 
tures), one can calculate these parameters for the other 
charge carrier. 


1, 447-453 (Jan. 
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Other NBS Publications 


Journal of Research 66A (Phys. and Chem.) No. 5 (Sept.-Oct. 
1962) 70 cents. 


Calorimetric calibration of an ionization chamber for deter- 
eo of X-ray total beam energy. J. S. Pruitt and 
S. Domen. (See above abstract. 

Zine a as a standard substance in the solution calorimetry 
of portland cement. E. 8. Newman. 

Nuclear optical model analysis of neutron elastic scattering 
for calcium. R. 8. Caswell 

Pyrolysis of some polyvinyl polymers at temperatures 
1,200 °C. SS. Straus and 8. L. Madorsky 

Lattice frequencies and rotational barriers for 
carbonates and nitrates from low temperature 
spectroscopy. R. A. Schroeder, C. E. Weir, 
Lippincott. 

Foreign gas broadening of the lines of hydrogen chloride and 
carbon monoxide. Kk. K. Plyler and R. J. Thibault. 

Monolayers of adipate polyesters at air-liquid interfaces. 
W. M. Lee, R. R. Stromberg, and J. L. Shereshefsky. 


up to 


inorganic 
infrared 


and E. R. 


Journal of Research 66B 


(Math. and Math. Phys.) 
(July—Sept. 1962), 


75 cents. 
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M. Newman. 


multivariate normal distributions. © I. 


Two theorems on matrices. 
Mill’s ratio 
Savage. 
Angle as a fourth fundamental quantity. 
Invalidity of Meixner’s theorem in 
dynamics. R. E. Nettleton. 

Selected bibliography of statistical literature, 1930 to 1957: 
VI. Theory of estimation and testing of hypotheses, 
sampling distributions, and theory of sample surveys. 
L. S. Deming. 


for R. 
J. KE. Romain. 
irreversible thermo- 


Journal of Research 66D (Radio Prop.) No. 5 (Sept.—Oct. 
1962), 70 cents. 


Theory of magneto-telluric fields. J. R. Wait. 

Propagation characteristics of magneto-ionic plasma columns. 
D. Formato and A. Gilardini. 

Dielectric loading of electric dipole antennas. J. Galejs. 

Possible influence of the ionosphere on the impedance of a 
ground-based antenna. J. R. Wait. 

Some statistical theory for the analysis of radio propagation 
data. M. M. Siddiqui. 

Auroral sporadic-F ionization. 
Owren 

Comparative study of the correlation of seasonai and diurnal 
cycles of transhorizon radio transmission loss and surface 
refractivity. B.R. Bean. 

Enhancement of the lunar tide in the noon critical frequency 
of the F, layer over the magnetic equator. R. G. Rastogi. 

Scattering from a conducting sphere embedded in a semi- 
infinite dissipative medium.. J. Galejs. 

High-frequency scattering from a coated sphere.  V. 
Weston and R. Hemenger. 

Propagation of spherical waves through an _ ionosphere 
containing anisotropic irregularities. K. C. Yeh. 


R. D. Hunsucker and L. 


H. 


Tables of spectral-line intensities. Part I. Arranged by 
elements, W. F. Meggers, C. H. Corliss, and B. F. Scribner, 
NBS Mono. 32, Pt. I (Dec. 29, 1961) $4.00. 

Radiation patterns in the lower ionosphere and Fresnel zones 
for elevated antennas over a spherical earth, R. G. Merrill 
and W. V. Mansfield, NBS Mono. 38 (Apr. 2, 1962) 70 
cents. 

An ultraviolet multiplet table, C. E. Moore, NBS Circ. 488, 
Sections 3, 4, and 5 (Apr. 6, 1962) Section 3, 60 cents; 

Section 4, 45 cents; Section 5, 30 cents. 

Specifications, tolerances, and regulations for commercial 
weighing and measuring devices (corrected through 1961), 
Handb. 44, 2d Ed. (1961) $1.00. 

Index to the reports of the National Conference on Weights 
and Measures. From the first to the forty-fifth, 1905 to 
1960, NBS Mise. Publ. 243 (June 28, 1962) 30 cent 
(Supersedes M203). 

Quarterly radio noise data—December, January, February 
1961-62, W. Q. Crichlow, R. T. Disney, and M. A. Jenkins, 
NBS Tech. Note 18-13 (PB151377-13) (1962) $1.75. 

Mean electron density variarions of the quiet ionosphere 
September W. Wright, L. R. Wescott, and D. J. 
Brown, NBS Tech. Note 40-7 (PB151399-7) (1962) $1.50. 

Mean electron density variations of the quiet ionosphere No. 
13; Summary of one year of data, May 1959—April 1960, 
J. W. Wright, NBS Tech. Note 40-13 (PB151399-13) 
(1962) $1.50. 

Bibliography on auroral radio wave propagation, W. Nupen, 
NBS Tech. Note 128 (PB161629) (1962) $2.75. 

Airborne television coverage in the presence of co-channel 
interference, M. T. Decker, NBS Tech. Note 134 (PB161635) 
(1962) $2.00. 

Ionosonde observations of artificially 
clouds: Firefly 1960, J. W. Wright, 
135 (PB161636) (1962) $2.50. 

A bibliography of the thermophysical properties of oxygen 
at low temperatures, J. C. Hust, L. D. Wallace, J. A. Crim, 
L. A. Hall, and R. B. Stewart, NBS Tech. Note 137 
(PB161638) (1962) $2.25. 

Vertical cross sections of the ionosphere across the geomag- 
netic equator, J. W. Wright, NBS Tech. Note 138 
(PB161639) (1962) $1.00. 

Siting criteria for HF communication centers, W. F. Utlaut, 
NBS Tech. Note 139 (PB161640) (1962) $1.25. 
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Atlas of Fourier coefficients of diurnal variation of foF2, 
W. B. Jones, NBS Tech. Note 142 (PB161643) (1962) 
$2.50. 

Numerical results for the surface impedance of a stratified 
conductor, C. M. Jackson, J. R. Wait, and L. C. Walter, 
NBS Tech. Note 143 (PB161644) (1962) $1.25. 

Analysis of ionospheric vertical soundings for electron density 
profile data. III. Procedures for obtaining monthly 
summary virtual height curves for N(H) analysis (com- 
posite virtual height curves), J. W. Wright, NBS Tech. 
Note 146 (PB161647) (1962) 75 cents. 

Dissociation constant of 2-ammonium-2-methyl-1,3-propane- 
diol in water from 0 to 50° and related thermodynamic 
quantities, H. B. Hetzer and R. G. Bates, J. Phys. Chem. 
66, 308-311 (1962). 

An introduction to flame photometry and a review of recent 
studies, M. Margoshes, Phys. Tech. Biological Research 4, 
215-260 (1962). 

Congruences for the partition function to composite moduli, 
M. Newman, Illinois J. Math. 6, No. 1, 59-63 (Mar. 1962). 

Research and the saving of teeth, G. C. Paffenbarger, J. 
Prostetic Dentistry 12, No. 2, 369-383 (Mar.-Apr. 1962). 

Comments on paper by W. D. Westfall, ‘“‘Prediction of VLF 
diurnal phase changes and solar flare effect.”’ J. R. Wait, 
J. Geophys. Research 67, No. 2, 916-917 (Feb. 1962). 

The shape of the geomagnetic field boundary under uniform 
external pressure, R. J. Slutz, J. Geophys. Research 67, 
No. 2, 505-513 (Feb. 1962). 

Plating standards and specifications, F. Ogburn, Electro- 
plating Eng. Handb. 2d ed., Ed. K. Graham, ch. 7, pp. 
257-262 (Reinhold Publ. Co., New York, N.Y., 1962). 

Inequalities for the permanent function, M. Marcus and M. 
Newman, Ann. Math. 75, No. 1, 47-62 (Jan. 1962). 

Study of electronically excited hydroxyl radicals in 
H-+0O,; atomic flame, H. P. Broida, J. Chem. Phys. 
No. 2, 444-448 (Jan. 1962). ; 

Theory of thermal diffusion in dilute alloys, R. E. Howard 
and J. R. Manning, J. Chem. Phys. 36, No. 4, 910-916 
(Feb. 1962). 

Effect of monomeric reagents on the melting (contraction) 
and recrystallization of fibrous proteins, L. Mandelkern, 
W. T. Meyer, and A. F. Diorio, J. Phys. Chem. 66, 375-376 
(1962). 

A correction to the exospheric electron density estimate using 
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